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Executive  Summary 

The  objectives  of  the  project  were  to  develop  basic  quantum  scientific  foundation  and  technical 
know-how  required  for  major  breakthroughs  and  innovations  in  “nanojet  technology,  which  will 
contribute  to  the  technology  developments  and  the  preservation  of  a  superior  leadership  of  US 
Airforce  in  frontier  quantum  nanojet  technologies. 

Nanojets  are  the  matter  jets  of  nanosize,  which  serve  as  a  unique  method  of  transporting 
particles,  emitted  from  a  single  or  multiple  jet  nozzle  of  nanometer  size  to  a  target  site  at  a 
designed  flow  rate,  with  design  particle  density  and  energy.  The  controlled  transport  and 
manipulation  of  the  streams  of  particles  offer  potential  applications  in  atom-molecular 
assembling,  nanofabrication,  single  electron  transistor,  atom  optics,  atom  lithography,  atomic 
interferometer,  drug  delivery  in  biomedical  application,  nano-propulsion  and  energy  sensors. 

Nanojets  are  classified  into  two  major  categories  depending  on  three  primary  factors: 
characteristic  dimensions  of  the  jets,  the  particle  density  and  injection  energy.  In  general, 
nanojets  of  the  liquids  of  high  molecules  weight,  operating  at  high  injection  energy  exhibit 
classical  jet  like  behavior  which  are  predicted  by  molecular  dynamics  or  classical 
Navier-Stokes  type  equation  with  analytical  formulation  based  on  lubrication  theory,  see  for 
example,  Moseler  and  Landman  and  Eggers.  However,  the  jets  of  electrons  or  particles  with 
sufficiently  low  atomic  weights  operating  at  low  injection  energy  will  exhibit  quantum  wave 
behavior,  which  are  markedly  different  from  the  conventional  liquid  jets.  Hence  the  classical 
models  cease  to  apply.  It  is  worthy  of  mention  that  in  prior  to  the  present  study  there  had  been 
little  or  no  in-depth  investigations  of  the  nanojets,  which  operate  in  the  region  where  the 
quantum  effects  play  the  dominant  dynamic  roles.  Thus  it  was  assessed  that  the  key  to  the 
successful  development  of  nanojet  technology  is  to  establish  the  fundamental  knowledge  of 
the  quantum  mechanical  structures,  dynamics  and  energetic  of  the  nanojets  operating  in  the 
ranges  of  interest  in  special  classes  of  particle  nanojets  applications. 

In  order  to  achieve  the  goal  of  quantum  nanojet  technology  described  above,  we  developed  a 
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fundamental  analytical  and  numerical  simulation  research  program  focused  on  the  jet  fluid 
dynamics  based  on  two  quantum  mechanical  perspectives:  Schrodinger’s  wave  mechanics 
and  quantum  fluid  dynamics  based  on  Hamilton-Jacoby  formalism.  Each  of  the  formalisms 
offers  unique  ontological  and  methodological  advantages,  which  provided  valuable 
complementary  knowledge  for  the  understanding  of  the  structures,  dynamics  and  energetic  of 
quantum  nanojets. 

We  accomplished  the  comprehensive  research  project  in  the  following  three  major  categories: 

I.  Quantum  Theory  of  Nanojets: 

Mathematical  formulation,  analytical  solution  methods  and  numerical  simulations  by 
modern  computational  schemes  for  Schrodinger’s  equation  and  quantum  fluid  dynamics, 
are  developed.  Principles  of  quantum  mechanical  equivalence  between  two  formalisms 
are  established. 

1)  Inverse  Problems, 

a) .  Schrodinger’s  wave  mechanics  are  used  to  construct  the  analytical  solution  of 

wave  fields  from  which  the  nanojet  flow  field  structures  are  predicted. 

References  2,  3,  4,  5,  6,  7,  8 

b)  Two-slit  jet  flow  field  structures  are  calculated  in  a  broad  jet  operating  range  of 
“quantum  Reynolds  number”,  defined  by  the  ratio  of  the  inertia  force  to  the 
quantum  force.  This  is  the  key  parameter  that  affects  the  jet  behaviors  and 
structures  quantum  behavior,  structures,  dynamics  and  energetic. 

References  2,  3,  4,  5,  6,  7,  8 

c)  Multi-slit  electron  jets,  formed  by  periodically  arranged  array  of  nanojets  nozzles, 
are  analyzed  in  broad  operating  range  of  quantum  Reynolds  number. 

References  8 

2) .  Direct  Problems 

a) .  Quantum  fluid  dynamics  formalism  based  on  Hamilton-Jacoby  equation  are 

adapted  for  the  numerical  simulation  and  analysis  of  nanojet  flow  fields. 

Reference,  6. 

b)  Quantum  computational  fluid  dynamics  numerical  codes  are  developed  for  flow  field 
calculation. 

Reference,  6. 

c)  Two-slit  jet  flow  field  structures  are  calculated  for  a  wide  jet  operating  ranges  of 
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quantum  Reynolds  number. 


Reference,  6. 

3).  Quantum  mechanical  Equivalence: 

Theory  of  quamlet,  see  nomenclature,  is  developed  to  examine  the  equivalence 
between  the  linear  quantum  wave  field,  described  by  Schrodinger’s  equation  and  the 
non-linear  quantum  flow  field  predicted  by  quantum  fluid  dynamics  formalism. 

Reference  3 

II.  Quantum  Nanojet’s  Structures,  Dynamics  and  Energetrics 

The  spatial  distributions  of  density,  velocity,  quantum  potentials  of  steady  state 
nanojets  are  predicted  in  the  broad  range  of  quantum  Reynolds  number  to  examine 
the  jets  configurations,  structural,  dynamic  and  energetic  behavior.  Quantum  bubbling, 
branching  and  clustering,  see  Glossary,  are  explained  by  wave  propagation  and  axial 
tunneling.  Interference  and  tunneling  and  their  physical  reasons  are  clarified. 

1 ) .  Sensitivity  Analysis 

a)  Sensitivity  of  two-  and  multi-slit  nanojet  structures  on  the  quantum  Reynolds 
number  Is  assessed.  Quantum  wave  Interference  Is  caused  by  multi-branch 
Interaction,  whereas  the  tunneling  is  due  to  the  axial  decay  in  the  probability 
density  when  the  kinetic  energy  associated  with  particle  motion  In  transverse 
direction  exceeds  Inlet  particle  energy. 

References  2,  4,  5,  6,  7,  8 

b)  Quantum  nanojets  are  classified  into  bubble  jets,  multi-branch  jets  and 
quasi-classical  jets. 

References  2,  4,  5,  6,  7 

c) .Key  nanojet  dynamic  properties:  jet  kinetic  energy  and  quantum  potential  and  their 

dependence  on  quantum  Reynolds  number  are  examined.  The  exchange  of  the 
quantum  potential  energy  and  mean  kinetic  energy  holds  a  key  to  the  local  jet 
acceleration  and  deceleration. 

References:  2,  4,  5,  6,  7,  8. 

2) .  Nanojet  Scaling  laws 

Reference  solutions  are  obtained  for  a  set  of  reference  jet  velocity  and  particle  density 
at  the  Inlets.  The  solutions  of  the  jets  at  other  operating  conditions  are  predicted  from 
the  scaling  laws. 

References  2,  4,  5,  6,  7 
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3).  Comparative  Study 


Results  of  the  jet  structure  sensitivity  analysis  of  two-slit  jets  based  on  two  formalisms 
are  compared  and  found  to  be  in  good  agreements. 

Reference,  6. 

III.  Quantum  Nanojet’s  Application. 

Quantum  nanojets  exhibit  number  of  interesting  quantum  structures  and  behaviors 
Involving  bubbling,  branching,  clustering,  interference  and  tunneling,  which  are  uniquely 
suitable  for  various  applications  in  nano-quantum  transport  of  sub-atoms,  atoms, 
molecules,  interferometer,  microscopes,  nanoscopes,  drug  delivery,  nanoelectronics  and 
nano-micro  scale  optical  tools  in  physical,  chemical,  biological  and  material  research,  and 
nanothrusters.  Proper  selection  of  injection  conditions  and  slit  geometry  is  essential  for 
each  application.  Examples  of  applications  are  described  below. 

Reference  3 

1 ) .  Low-quantum  energetic  nanojet  systems.'  Strong  quantum  effects. 

Applications."  nano-sensors,  nano-imaging,  mass-spectrometer*,  energy  sensor 
interferometer,  cryptography,  nano-enrgetic,  nanowires,  nanoelectronics, 
nanooptoelectronics. 

2) . Intermediate  quantum  energetic  nanojets;  Moderate  quantum  effects 

Applications:  Nanofabricaltion,  soft  lithography,  jet-assisted-nano  assembling, 
spectrometer*,  interferometer,  nanoelectronics 

3) . High  quantum  energetic  nanojets.  Weak  quantum  effects 

Applications:  Nano  &  micro  propulsion,  see  item  a),  below.  Particle  &  molecular  beams, 
hard-lithography  quantum  computation 

In  closing,  we  have  successfully  completed  a  comprehensive  research  project  in  a  profoundly 
important  area  of  nanoscience:  Quantum  Nanojets  Technology.  The  one-year  efforts  led  to  the 
establishment  of  basic  knowledge  and  fundamental  data  of  nanojet  physics  and  also  provided 
remarkable  insights  for  the  potential  future  developments  of  nanoscientific  devices:  mass 
spectrometer,  nanoimaging,  lithography,  nano-propulsions,  nanowires,  and  nanotube 
technology.  Among  the  most  important  findings  are  the  critical  roles  played  by  the  quantum 
potential  and  the  quantum  potential  forces,  which  are  induced  by  the  quantum  diffusion  of  the 
probability  fluid.  It  is  remarkable  to  point  out  that  the  quantum  structures  and  dynamic 
behaviors  of  two-slit  nanojet  are  strongly  dependent  on  the  quantum  Reynolds  number,  which 
is  defined  as  the  ratio  of  the  inertia  force  to  the  quantum  potential  force  acting  on  the  fluid. 
Nanojets  exhibits  two  unique  structural  and  dynamic  behavior:  1 )  quantum  bubbling,  caused 
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by  premature  break-up  of  the  jet  at  low  Reynolds  number,  2)  branching,  Induced  by  multi 
branch-jets  interaction  which  yields  Interference  fringes,  and  3)  quantum  clustering:  formation 
of  the  high  density  cluster  core  along  the  jet  axis,  which  yields  the  nanojets  “signatures”.  The 
clustering  Is  due  to  the  intersection  of  right  and  left  running  waves  from  the  silt  surface  which 
creates  compression  force,  which  facilitates  the  formation  of  the  high  density  cluster  core 
along  the  jet  axis. 

These  unique  dynamics  and  energetic  behavior  of  nanojets  offer  broad  areas  of  applications 
as  nanodevices,  sensors,  imaging,  and  space  thrusters. 

The  present  writer  and  his  associates  wish  to  extend  their  appreciation  for  the  unique 
opportunity  in  developing  this  remarkably  Interesting  branch  of  new  nanoscience  under  the 
sponsorship  of  the  US  Air  Force  Office  of  Scientific  Research.  We  acknowledge  that  part  of  the 
works,  in  particular  Ref.  1,  7  and  8,  have  been  supported  by  National  Science  Council  of 
Taiwan.  Special  thanks  are  due  to  Dr.  F.L.  Madarasz  for  his  valuable  contributions  throughout 
the  research  works.  Finally,  I  should  like  to  warmly  thank  all  my  colleagues  Professor  S.Y.  Lin, 
Professor  K.R.  Chen  and  our  graduate  students  Mr  C.T  Lin,  Mr  Y.S.  Lin  and  Mr  TS  Lin  for  their 
technical  contributions  toward  successful  completion  of  the  research  project.  I  should  also  like 
to  thanks  Dr.H.L.  Tsai  and  Ms  Evensa  Su  for  their  exceptional  administrative  assistance. 

Nomenclature 

Symbol: 

b:  half  of  silt  separation  in  two-slit  jet 
E:  particle  energy  at  the  exit  of  slit  =  (1/2m)U^ 
m;  mass  of  a  particle 

N:  Quantum  Reynolds  number,  N=  (2mE)^^^  b/rth  *=  Ub/D  ( Two-sllt ) 

Ns:  Quantum  Reynolds  number,  N=  (2mE)^^^  R/ttIi  *=  UR/D  (Single-slit) 

R:  half-width  of  a  jet. 

U:  particle  velocity  at  the  slit. 

Sub-script 

cr:  critical  condition,  I.e.,  transition  from  quantum  to  classical  state 

Glossary 


•  Quantum  Reynolds  number: 

Quantum  Reynolds  number  represents  the  ratio  of  the  inertia  force  to  the  quantum 
force  acting  on  the  fluid.  (  Chiu  2002  ) 

•  Quantum  nanoJets 

Quantum  nanojet  is  defined  as  the  streams  or  arrays  of  particles  such  as  electrons. 
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atoms,  molecules,  neutrons.  Ions  or  photons.  Injected  directly  or  through  a  Injector  of 
nanoscale,  to  a  designated  target  space  In  the  presence  or  absence  of  an  external 
potential,  and  thereby  exhibits  quantum  phenomena.  (Chiu  2002) 

•  Quantum  potential 

Quantum  potential  is  defined  as  the  sum  of  the  kinetic  energy  and  the  dilatation 
energy  associated  with  quantum  diffusion  velocity  vector  field  of  quantum  fluid. 
(Chiu  2002,  Reference  1 ) 

•  Quantum  potential  force 

The  force  derived  from  the  gradient  of  the  quantum  potential. 

•  Quamlets 

Quantum  flowlets  associated  with  a  single  eigen-state  or  the  linear  superposition  of 
two-degenerate  states.  (Reference  3) 

•  Quantum  nanojet  phenomena 

Three  important  quantum  phenomena  observed  are  quantum  bubbling,  branching 
and  clustering  in  quantum  nanojets  at  specific  values  of  quantum  Reynolds 
number. 


Selected  Highlights  of  nanojets 

I.  “Nanojet  at  a  given  quantum  Reynolds  number  exhibits  unique  structural  and  dynamic 
behavior” 

•  Structural/dynamic  complexities  are  the  principal  features  of  “Quantum  Nanojet  Energetic 
Hierarchical  Systems” 

a) .  Strong  quantum  jets: 

i)  Bubble  jets,  N<5 

ii )  Multi-branch  fan  jets, 

quantum  force  » inertia  force  N  *«  N 

iii)  Diffusion  dominated  flow 

b) .  Intermediate  quantum  jets: 

i)  Two-branch  quantum  jets 

quantum  force  ~  Inertia  force  N  ~  Ncr 

ii)  Mixed  diffusion-convection  flow 

c) .  Two-primary  quasi  classical  jets: 

i)  Two-primary  classical  jets 

quantum  force  «  Inertia  force  ~  N  »  Ncr 
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ii)  Convection  dominated  flow 

•  Passive  control  of  nanojets: 

The  structure  and  dynamic  behaviors  of  nanojets  can  be  controiied  by  the  selection  of 
appropriate  quantum  Reynolds  number  through  the  proper  choice  of  any  or  ail  of  the 
following  parameters. 

i)  inlet  particle  energy 

ii)  Slit  width 

iii)  Size  of  jet  of  each  slit. 

•  Active  control  of  nanojets: 

The  structural  and  dynamic  behaviors  of  nanojets  can  be  controlled  by  the  proper 
application  of  a  tailored  external  potential,  such  as  an  electrostatic  and/or  a  magnetic 
potential.  The  active  control  plays  significantly  important  roles  in  the  innovation  and 
application  of  new  nanodevices. 

II  Quantum  bubbling, 

“Formation  of  “bubbles”  which  are  characterized  by  bubble-like  structures  appear  when  the 
quantum  number  is  very  low.  The  bubble-like  structures  occur  when  the  jet  momentum  force  is 
smaller  than  quantum  potential  force.  The  jet  breaks  up  in  discrete  bubbles. 

III.  Quantum  branching/confluencing: 

“Phenomena  displaying  the  formation  of  multi  branches  jets  interaction,  which  yield 
interference  fringe  formation”. 

•  Quantum  wave  harmonic  generation,  wave  interactions  and  convective  motion 
inter-coupling  are  the  principal  mechanisms  of  jet  state  transition  and  interference  fringe 
formation 

a) .  Phenomena: 

(i) .  Jet-splitting  results  in  the  formation  of  multi  branches  at  low  to  intermediate 

Reynolds  number. 

(ii) .  Center  branch  creates  bright  fringe  on  the  central  axis  at  low  to  intermediate 

Reynolds  number. 

(iii)  Formation  of  jet-streamers  on  the  nanojet  surface  at  low  to  intermediate 
Reynolds  number. 

(iv) .  Interference  fringes  of  various  patterns:  single,  double  peaks,  Mexican  hat, 

M-hat  and  easter-hat  fringe  patterns. 

b) .  Mechanism: 
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(i).  Harmonic  wave  generation  and  convective  motion  intercoupiing. 


(ii) .  Intersection  of  waves  form  the  upper  and  lower  slit  jets  on  the  plan  of  symmetry 
creates  a  center  branch  at  low  to  intermediate  Reynolds  number. 

(iii) .  Excitation  of  left  and  right  running  waves  along  the  characteristic  surfaces 
c).  Dynamic  and  structural  impacts: 

(i) .  Creation  of  rapidly  varying  quantum  force  in  both  transverse  and  axial  direction  in 
Fresnel  region. 

(ii) .  Complex  fringe  patterns  caused  by  quantum  potential  under  branching  and 

confluencing. 

IV.  Quantum  clustering: 

“Formation  of  high  density  particle  jet  core,  which  serves  as  quantum  nanojets  “signatures”” 
•  Signature  of  nanojet  is  formed  by  quantum  clustering 

(a) .  Phenomenon: 

(i) .  High-density  particle  clusters:  base  and  primary  clusters  appear  in  the  jet  cores. 

The  size,  location  and  shape  of  cluster  are  the  structural  signature  of  a  nanojet. 

(ii) .Size  and  the  location  of  cluster  increases  with  respect  to  the  quantum  Reynolds 
number. 

(b) .  Mechanism: 

(i) .  Intersection  of  right  and  left  running  waves  from  the  slit  surface  compresses  the 

probability  fluid  and  creates  high  density  clusters  on  the  nanojet  axis. 

(ii)  Convective  motion  shifts  the  location  and  size  of  the  clusters  toward  the 
downstream  of  nanojet  flow. 

(iii) .  Clusters  decay  in  axial  direction  by  quantum  tunneling. 

(c) .  Dynamic  and  structural  impacts: 

(i) .  Creation  of  rapidly  varying  quantum  axial  force  due  to  quantum  tunneling 

(ii) .  Cluster  spreads  along  axial  direction  and  ultimately  disperses  in  the  full  nanojet 
core  at  large  value  of  N. 

V.  Invention  of  NanoJet  based  mass  spectrometers  and  energy  sensors" 

“Based  on  the  mechanisms  of  quantum  branching,  confluence  and  clustering,  the 
following  four  types  of  nanojet-based  devices  can  be  developed  for  practical  applications 
in  nanoscience  and  quantum  science”. 
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(a) .  Mass  spectrometer:  Low  Reynolds  number  device  (Center-branch  based  design). 

(b) .  Mass  spectrometer:  Intermediate  Reynolds  number  device  (Cluster  based  design). 

(c) .  Energy  sensors:  Low  Reynolds  number  device  (Center-branch  based  design) 

(d) .  Energy  sensors:  Intermediate  Reynolds  number  device  (  Cluster  based  design) 

1)  Mass  spectrometer  and  energy  sensor: 

“Low  Reynolds  number  devices” 

•  Center  branch-based  mass  spectrometer. 

(a).  Principle: 

The  center  branch  of  a  dual-slit  jet  starts  at  an  axial  location  at  distance,  S,  form  the 
center  of  two  slits.  The  S  is  dependent  on  N,  i.e.,  S  (N)-function,  given  by, 

S(N)  =  So+aN“ 

where  So,  a,  and  a  are  constants  for  a  given  two-slit  geometry. 

(b).  Method: 

When  S  is  obtained  from  the  experimental  measurement,  the  mass  of  a  particle,  m, 
can  be  calculated  from  the  above  equation  for  given  particle  energy,  E,  at  the  slit,  as 
follows: 

m  =  (1/2E)  (nt)/bf  [(S-So)/a] 

•  Center-branch  jet  based  energy  sensor. 

(a) .  Principle: 

By  the  same  principle  described  above,  we  can  predict  the  energy  of  the  particle,  E, 
at  the  slit,  if  S  and  m  are  known, 

E  =  (1/2m)(7ih/b)^  [(S-S)/a] 

(b) .  Method: 

When  S  is  measured,  E  can  be  calculated  from  the  above  equation  for  a  given 
atomic  weight 

2)  Mass  spectrometer  and  energy  sensor: 

“Intermediate  Reynolds  number  devices” 

•  Cluster  signature  based  mass  spectrometer. 

(a).  Principle: 

The  length  of  a  cluster,  L,  of  dual-slit  jet  is  given  as  a  function  of  N,  i.e.,  L(N)-fu notion, 
as  follows, 

L(N)  =  Lo+bNP 
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where,  U,  b  and  p  are  constants  for  a  given  two-slit  jet  and  m  is  the  atomic  weight. 

(b).  Method: 

When  L  is  measured  from  the  cluster  signature,  m  can  be  calculated  from  the  above 
equation  for  given  particle  energy,  E,  at  the  slit. 

m  =  (1/2E)(7i  h/bf  [(L-Lo)/b] 

•  Cluster  signature  based  energy  sensor 
(a).  Principle: 

By  the  same  principle  described  above,  the  energy  of  the  particle  at  the  slit  exit  is 
given  by  the  following  expression, 

E  =  (1/2m)(7t1i/bf  [(L-Lo)/b] 

(b).  Method: 

When  L  is  obtained  experimentally,  E  can  be  calculated  from  the  above  equation  for 
a  given  atomic  weight  m. 
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Quantum  fluid  mechanical  and  electronic  structure  of  a 

hydrogen-like  atom 

By  H.  H.  Chiu 

Space  Science  and  Technology  Center  and  Institute  of  Aeronautics  and  Astronautics 
National  Cheng  Kung  University,  Tainan,  70101,  Taiwan 
(hhchiu@htindl  .iaa.ncku.edu.tw) 

Quantum  diffusive  fluid  dynamics  formalism  is  adapted  to  examine  the  quantum  mechanical 
significance  of  quantized  energy  and  the  mechanisms  of  detailed  dynamic  equilibria  among  the 
Coulomb  field,  vortex  induced  flow  and  the  quantum  potentials  of  elemental  processes  induced  by 
the  quantum  diffusion  of  the  probability  fluid  of  the  electron  in  a  hydrogen-like  atom.  The 
quantized  energy  is  found  to  be  the  negative  of  the  kinetic  energy  of  the  asymptotic  diffusion 
velocity  of  the  probability  fluid  of  the  electron.  The  Coulomb  force  balances  with  the  forces  due  to 
the  total  radial  quantum  potential  induced  by  radial  diffusion  minus  the  centrifugal  barrier  force. 
With  a  finite  magnetic  quantum  number,  a  line  vortex  is  present  in  an  atom  structure.  The  kinetic 
energy  of  the  vortex-induced  flow  is  supplied  by  a  part  of  the  quantum  dilatation  energy  associated 
with  the  diffusion  in  polar  angular  direction,  whereas  the  remainder  of  the  dilatation  energy 
balances  with  the  centrifugal  barrier  potential.  An  atom  with  a  finite  magnetic  quantum  number  is 
structured  with  the  vortex  induced-flow  dominated  inner  core  with  the  line  vortex  extending 
beyond  the  Bohr  radius  along  the  polar  axis,  and  the  electrostatic  potential  dominates  the  exterior 
region  of  the  core.  In  a  state  of  zero  magnetic  quantum  number  the  atom  exhibits  a  pure  electronic 
structure.  The  paper  concludes  with  areas  of  future  research  on  quantum  energetic  and  dynamic 
equilibria  in  other  quantum  systems  and  the  validity  of  the  model. 

Keywords:  quantum  diffusive  fluid  dynamics;  quantum  energetic;  hydrogen  atom  structure 


1.  Introduction 

Three  profound  issues,  whieh  remain  un-addressed  or  imperfectly  understood  in  a  hydrogen  atom, 
are:  the  quantum  mechanical  significance  of  the  quantized  energy,  the  mechanism  of  the  dynamic 
balance  of  an  electron  under  the  Coulomb  potential  field,  and  the  source  of  the  kinetic  energy  of 
the  vortex  induced  flow,  in  a  stationary  state,  defined  by  a  set  of  quantum  numbers.  An  atom  emits 
a  photon  when  it  is  de-excited  from  a  higher  energy  level  to  a  lower  state  and  absorbs  quanta  when 
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it  is  excited  by  radiative  absorption  or  atomic  collision. 

The  questions  of  exactly  what  the  nature  of  the  quantized  energy  is  that  accompanies  the 
differential  changes  of  energy  in  atomic  and  molecular  absorption  and  emission,  and  why  the 
energy  of  a  stationary  eigenstate  takes  on  a  particular  constant  value,  however,  have  curiously 
remained  unanswered  despite  the  successful  prediction  of  the  quantized  energy  levels  themselves 
by  quantum  theory.  Additionally,  the  nature  of  dynamic  equilibrium  for  maintaining  the  electron 
without  falling  upon  nucleus  under  the  Coulomb  attractive  force  at  the  discrete  quantized  energy 
levels  of  an  atom  has  been  inadequately  understood.  Another  dynamic  phenomenon  involving  a 
line  vortex  induced  flow  in  an  atom  under  a  finite  magnetic  quantum  number  provokes  the 
question  of  how  does  the  latter  flow  acquire  the  energy  to  sustain  the  motion.  A  review  of 
literature  reveals  that  the  quantitative  criteria  for  dynamic  equilibria  and  detailed  structure  of  a 
global  atomic  system  in  a  stationary  state  are  not  known  to  a  satisfactory  level  of  detail.  It  is 
envisioned  that  the  quantum  phenomena  related  with  the  above  three  principal  modes — quantized 
energy,  Coulomb  potential,  and  vortex  motion — are  interrelated  issues,  which  need  to  be  treated 
concomitantly  by  what  may  be  described  as  global  modal  balance  in  an  atomic  system.  The 
objectives  of  the  present  study  are  to  develop  a  self-consistent  atomic  modal  balance  theory  based 
on  the  frameworks  of  quantum  diffusive  fluid  dynamics,  and  to  examine  the  nature  of  energetic 
and  dynamic  equilibria  of  a  hydrogen  atom  in  a  stationary  state  so  as  to  establish  a  dynamic 
characterization  of  quantized  energy  and  to  gain  a  basic  understanding  of  the  electronic  and  fluid 
dynamic  structures  of  an  atom  from  the  quantum  fluid  dynamics  perspective. 

In  quantum  theory  the  hydrogen  atom  is  taken  as  a  two-particle  system,  a  proton  and  an  orbital 
electron.  However,  in  the  fluid  dynamics  model  presented  here,  the  electron  is  viewed  as  a 
probability  fluid  whose  density  is  the  product  of  the  wave  function  and  its  complex  conjugate, 
which  both  satisfy  the  Schrodinger  equation.  The  quantum  probability  fluid,  whose  stress  tensor  is 
constructed  from  the  quantum  potential  induced  by  non-uniform  spatial  density  distribution,  obeys 
the  laws  of  conservation  of  quantum  fluid  dynamics.  There  are  two  distinct  developments  in  the 
formalism  of  quantum  fluid  dynamics  or  the  quantum  theory  of  motion:  the  Hamilton-Jacobi 
transformation  applied  to  the  Schrodinger  equation;  and  alternatively,  the  formulation  by 
stochastic  mechanics  of  the  classical  theory  of  Brownian  motion  of  particles.  The  development  of 
two  formalisms  and  historical  highlights  are  described  below. 

Quantum  fluid  dynamics  and  the  quantum  theory  of  motion  are  usually  viewed  as 
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generalizations  of  Newton's  classical  mechanics  in  the  Hamiltonian  formalism.  They  have  also 
been  known  for  sometime  via  earlier  theories:  pilot  wave  of  de  Broglie  (1928);  the  hidden-variable 
theory,  Bohm  (1952a,  1952  b);  the  causal  interpretation  and,  later  the  ontological  interpretation  of 
quantum  mechanics,  Belinfante  (1973),  Bohm  &  Hiley  (1993);  and,  most  recently,  the  theory  of 
quantum  motion.  An  excellent  book  by  Holland  (1993)  presents  a  comprehensive  quantum 
theoretical  account  of  the  de-Broglie-Bohm’s  causal  interpretation  of  quantum  mechanics 
covering  the  broad  exposition  of  the  theory  and  critical  review  of  specific  issues:  particle  trajectory, 
quantum  prediction,  classical  deterministic  model,  non-locality,  mathematical  features,  quantum 
intuition,  and  a  particle’s  reciprocal  action  on  the  associated  wave.  Of  particular  interest  are 
historical  dialogues  and  contingencies  surrounding  the  Einstein-Podolosky-Rosen  paradox  (1935), 
Copenhagen  interpretation  by  Bohr  (1958),  interconnectedness  theorem  of  Bell  (1966),  and 
experiments  on  Bell’s  inequalities  by  Aspects  (1982).  While  we  do  not  elaborate  on  a  detailed 
description  of  each  issue  cited  above,  we  will  focus  our  attention  on  the  unique  concept  of  a 
quantum  potential,  which  arises  through  the  Hamilton-Jacoby  formalism,  and  its  significance  to 
the  nature  of  the  global  atomic  problems  we  intend  to  resolve.  The  rationale  for  this  motivation  is 
supported  by  the  observation  made  by  Holland  (1993),  who  foresaw  that  “corresponding  to  the 
new  theory  of  motion,  i.e.,  generalized  Hamilton-Jacoby  theory,  is  a  new  conception  of  matter.”  In 
following  this  observation  we  first  acknowledge  that  a  salient  feature  of  the  new  theory  is  that 
matter  has  an  intrinsic  field  aspect:  the  mass  points  moving  and  interacting  under  the  influence  of 
the  quantum  potential,  which  is  characterized  as  an  organizational,  self-referential  form  of 
potential  energy  that  brings  about  inner  stresses  in  the  material  system  to  which  the  mass  points 
respond,  Holland  (1993).  Secondly,  we  note  that  the  quantum  potential  has  a  defining  feature  of 
being  a  holistic,  context-dependent  nonlocal  potential  that  depends  only  on  the  "form"  of  the 
potential,  and  not  its  intensity  or  strength  and  the  effects  of  such  potential  do  not  fall  off  with 
distance,  Bohm  and  Hiley  (1993).  These  two  unique  aspects  of  the  quantum  potential  appear  as  a 
fundamental  basis  in  formulating  a  concomitant  global  balance  theory  with  the  quantum  potential 
giving  rise  to  key  dynamic  and  energetic  factors  that  govern  the  structure  of  a  hydrogen-like  atom. 

Other  quantum  fluid  descriptions  developed  by  Nelson  (1966,  1985),  Comisar  (1965),  and 
Kershaw  (1964)  are  based  on  stochastic  processes  via  the  classical  theory  of  Brownian  motion.  An 
extensive  discussion  of  quantum  fluctuations  and  a  list  of  references  can  be  found  in  excellent 
lecture  notes  published  by  Nelson  (1985).  The  stochastic  mechanical  formulation,  or  Nelsonian 
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mechanics,  does  not  share  the  physical  basis  of  a  causal  theory  but  adapts  a  probability  concept  of 
classical  theory,  through  which  atomic  processes  are  described  by  particle  motion  in  space  and 
time.  However,  despite  of  the  profoundly  different  theoretical  conceptual  framework  of  classical 
stochastic  mechanics  and  quantum  mechanics,  the  theory  shows  that  for  observations,  which  may 
be  reduced  to  position  measurements,  the  two  formalisms  give  the  same  quantum  mechanical 
prediction.  Nelson  (1966).  Two  unique  steps  have  been  undertaken  in  the  stochastic  mechanical 
formalism.  They  are  the  introduction  of  quantum  osmotic  velocity.  Nelson  (1966),  which  is  the 
negative  of  the  diffusion  velocity  defined  by  Chiu  (2002),  and  the  formulation  of  the  evolution 
equation  of  the  osmotic  velocity  and  the  momentum  equation.  With  the  use  of  a  Hamilton-Jacobi 
transformation,  the  equation  governing  the  evolution  of  osmotic  velocity  can  also  be  obtained 
from  the  continuity  equation  for  the  probability  current  density,  as  we  shall  show  later.  Thus  the 
laws  of  conservation  derived  from  the  two  formalisms  are  mathematically  equivalent  to  each  other. 

While  we  will  not  elaborate  on  historical  dialogues  and  contingencies  between  the 
Copenhagen  vs.  de-Broglie-Bohm  interpretations  of  quantum  mechanics,  we  will  adopt  a  view 
that  is  based  on  an  examination  of  the  measurement  process,  which  suggests  that,  within  a  limited 
framework,  the  two  theories,  quantum  mechanics  and  quantum  fluid  dynamics,  are  equivalent. 
Nelson  (1966).  Thus  the  latter  approach,  i.e.,  the  quantum  fluid  dynamics,  is  used  as  an  alternative 
methodological  framework  for  the  proposed  investigation  of  atomic  modal  balance.  In  §2  we 
recapitulate  the  quantum  diffusive  fluid  mechanics  formulation  of  the  laws  of  conservation  of 
probability  density  and  energy.  In  §3  we  elucidate  the  role  of  a  quantum  potential  and  how  its 
presence  gives  rise  to  the  dilatation  and  kinetic  energies  associated  with  diffusion  field,  which  are 
in  accord  with  the  general  framework  of  Nelson’s  stochastic  formulation  description.  Principal 
modes  and  the  elemental  processes  are  introduced  to  facilitate  modal  balance  analysis  in  §4.  The 
nature  of  atomic  modal  energy  balance  and  dynamic  equilibria  are  investigated  in  §5.  This 
analysis  includes  the  examinations  of  physical  significance  of  the  quantized  energy,  the  detailed 
balance  among  the  Coulomb  field,  vortex  induced  flow,  quantum  potential,  and  the  intercoupling 
modes  of  an  atom.  In  §6  we  present  numerical  results  including  the  detailed 
electronic-hydrodynamic  structures  for  an  atom  for  selected  quantum  states  and  discuss  the  range 
of  validity  of  the  results.  The  paper  concludes  with  the  suggested  areas  of  future  research  in 
quantum  energetic  of  other  atoms,  molecules  and  many-particle  systems. 
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2.  Quantum  Bernoulli’s  equation 

The  hydrogen  atom,  a  proton  and  an  electron,  is  a  two-particle  system  with  a  central  force  acting 
on  the  line  between  the  two  particles.  Canonically,  the  Schrodinger  equation  describing  this 
system  is  separated  into  two  independent  different  equations:  one  for  the  center-of-  mass  motion, 
and  the  other  for  the  relative  motion  about  the  center-of-mass.  The  center-of-mass  problem  is 
treated  as  a  single  particle  with  a  mass  of  M  =  (me+mp),  where  mp  and  mg  are  the  masses  of  the 
proton  and  electron,  respectively,  translating  through  space.  The  relative  motion  problem  is  also 
treated  as  a  single  particle  problem  but  with  a  reduced  mass  defined  by  p=  (mpnig  )/M.  The  ratio  of 
the  proton  to  electron  mass  is  mp/mg  =1837.  Since  the  proton  is  much  more  massive  than  the 
electron  it  is  located  very  close  to  the  center-of-mass,  while  the  reduced  mass  is  close  to  the 
electron’s  coordinates.  Following  Landau  and  Lifshitz  (1985),  we  approximate  the  center  of  mass 
as  the  origin  of  the  relative  coordinate.  Because  of  the  above  approximation,  the  solution  is,  in 
principle,  invalid  in  a  small  spherical  region  of  radius  (mE/mp)re  centered  at  the  origin  of  the 
coordinate,  as  in  the  similar  treatment  presented  in  standard  texts.  Thus  for  the  hydrogen  atom,  the 
time  dependent  Schrodinger’s  equation  for  relative  motion  in  Coulomb  potential,  Ve  ,  is  given  by, 

rt^-~vV+Vr  (2.1) 

Ot  2m^ 

where  the  reduced  mass,  //,  has  been  approximated  by  electron  mass  nie  and  h  is  the  Planck’s 
constant  divided  by  27i. 

In  the  present  quantum  fluid  dynamics  model,  we  assume  that  the  electron  consists  of  the 
Schrodinger’s  probability  fluid,  attracted  by  the  Coulomb  force  induced  by  the  positively  charged 
nucleus.  The  probability  fluid  obeys  the  laws  of  conservation  derived  by  the  Hamilton-Jacobi 
method.  Following  Bohm  (1952a)  we  first  express  the  wave  function  in  the  polar  form, 
y/  =  / K) ,  where  p  is  the  probability  density,  (j) ,  is  a  scalar  potential.  By  inserting  the 

latter  expression  of  y/,  into  Eq.  (2.1),  and  separating  into  imaginary  and  real  parts,  we  obtain  the 
equations  for  the  flow  variables  p  and  (f) .  The  imaginary  and  real  part  gives  the  conservation 
equation  of  probability  density  and  energy,  respectively. 
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^  +  V*J  =  0, 

dt 


(2.2) 


dt 


f-1 

UJ 


u«u 


(2.3) 


where  /  is  the  current  density  or  Lagrangian  density.  The  expressions  for  p  and  /  are: 


P  =  y/y/\ 


(2.4) 


and 


J  =  /7U  = 


2im„ 


(yy  y/  -\lN ip")  , 


(2.5) 


where,  u,  is  the  mean  convective  velocity,  defined  by  the  ratio  of  the  probability  current  flux  J  with 
the  density. 


u  = 


2im 


¥  ) 


e  J 


=  V(I). 


(2.6) 


The  velocity  field,  Eq.  (2.6),  is  irrotational,  except  at  the  nodal  singularities  where  the  wave 
function  vanishes  and  vorticies  are  present.  It  has  been  shown  by  Carlen  (1984)  Nelson  (1985), 
Bacciagaluppi  (1998),  and  Bemdl,  et  al  (1995),  that  under  certain  regularity  conditions  on  the 
initial  wave  function  and  the  external  potential,  the  particle  has  zero  probability  of  entering  the 
nodal  surface  from  outside.  Since  the  determination  of  the  current  density,  J  ,  Eq.  (2.5)  requires 
the  simultaneous  measurements  of  position  and  velocity  of  particle,  which  is  not  in  compliance 
with  the  uncertainty  principle,  J,  can  not  to  be  treated  as  the  averaged  measured  particle  flux  at  the 
point  in  space  and  time  interpreted  in  classical  fluid  dynamics,  except  in  the  case  when  J  varies 
slightly  or  not  at  all  in  a  spatial  coordinate  so  that  an  accurate  velocity  determination  can  be  made 
without  impairing  the  usefulness  of  the  concept  of  flux. 

The  quantum  potential,  appears  in  the  fourth  term  on  the  left  side  of  Eq.  (2.3),  has  the 
functional  form. 
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(2.7) 


(\\ 

V"(ln/7)+  -  V(ln/7).V(ln/7) 

v2y 


We  will  adapt  the  form  of  quantum  potential  and  for  notational  convenience  will  denote  it  by 

just  n  . 

The  momentum  equation  is  obtained  by  taking  the  gradient  of  Eq.  (2.3),  as  follows,  see  for 
example,  Ghosh,  &  Deb.  (1982), 


0u  VV 

—  +u*Vu= - ^+Vn  (2.8) 

dt  m^ 


This  equation  serves  as  guide  for  the  examination  of  dynamic  criteria  of  fluid  equilibria. 


3.  The  roles  of  quantum  diffusion  in  modal  dynamic  balance 

The  Quantum  potential,  IT,  plays  a  major  role  in  quantum  energetic  involved  in  the  large  classes  of 
quantum  flow  fields  with  non-uniform  density  distributions.  In  this  section  we  will  characterize  it 
within  the  context  of  quantum  diffusion  processes.  We  introduce,  in  analogy  with  classical  mass 
transport  theory,  the  diffusion  of  the  probability  fluid  characterized  by  the  gradient  induced 
diffusion  velocity,  V ,  Chiu  (2002),  as  follows, 


^  n  ^ 


2m 


V  In  =  -D 


e  J 


p 


(3.1) 


where  5  is  the  diffusion  current  flux  of  the  probability  fluid,  and  the  quantity  h  /  2m^ ,  which  has  the 

same  physical  dimension  as  the  mass  diffusivity,  is  termed  as  quantum  diffusivity,  designated  by  D. 
It  is  of  special  interest  to  note  that  Nelson  (1966,  1985)  defines  the  osmotic  velocity  as  the  mean 
value  of  the  difference  between  the  forward  and  backward  velocities  of  Brownian  particle  motion 
with  diffusion  coefficient  D,  and  no  friction.  According  to  Einstein’s  theory  of  Brownian  motion, 
the  osmotic  velocity  is  the  velocity  acquired  by  a  particle  in  equilibrium  with  respect  to  an  external 
force,  to  balance  with  the  osmotic  force.  By  adapting  the  available  data,  for  example,  D 
=  0.116x10'^  m^/sec  and  Bohr  radius,  a,  of  a  hydrogen  atom  of  5.29xl0'“’m,  we  estimate  a 
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characteristic  diffusion  velocity  to  be  approximately  2.3x10*’  m/sec. 


The  physical  significance  of  the  first  term  of  the  quantum  potential,  {h^  j  In  p ,  can  be 

obtained  from  the  definition  of  the  diffusion  velocity,  as  follows.  Taking  the  divergence  of  Eq.  (3.1) 
yields. 


dJl  = 


Ami  j 


V^(ln/7)  = 


[  P  J 


=  -Z)V.V 


(3.2) 


By  comparison,  it  is  seen  that  the  first  term  in  the  quantum  potential  per  unit  mass,  Eq.  (2.7), 
represents  the  specific  dilatation  energy  of  diffusion  field,  dJt  =  -Z)V»V  ,  due  to  the  divergence  of 
the  quantum  diffusion  velocity.  The  specific  dilatation  energy,  Eq.  (3.2),  agrees  with  the  negative 
of  the  dilatation  of  the  osmotic  velocity  that  appears  in  Nelson’s  stochastic  mechanical  formulation. 
Nelson  (1966).  The  dilatation  of  the  diffusion  velocity  field  has  a  different  physical  significance 
from  that  of  the  mean  convective  velocity,  discussed  below. 


The  specific  dilatation  of  a  mean  convective  velocity  is  given  by  equation,  (2.1)  as. 


-V.„.^+„.VInp-^-5^  =  -^  (3,3) 

dt  dt  D  dt 

which  is  equal  to  substantive  derivative,  d/dt ,  describing  the  co-moving  time  rate  of  change  of 
expansion  or  contraction  of  the  fluid  seen  from  the  instantaneous  rest  frame  of  the  actual  fluid 
parcel.  According  to  Eq.  (3.3)  we  conclude  that  a  steady  state  flow  is  a  “convective  dilatation  free”, 
when  u  and  V  are  normal  to  each  other  throughout  the  flow  field.  We  shall  see  later  that  the 
probability  fluid  flow  field  in  a  hydrogen  atom  at  a  stationary  state  described  by  a  set  of  quantum 
numbers  is  convective  dilatation  free. 


Similarly,  the  specific  dilatation  of  the  diffusion  velocity  is  obtained  by  taking  the  divergence 
operator  on  (3.1)  as  follows, 


-V.V  =V.Vln/7+(— )  V  V 
P 


(3.4) 


The  second  term  appears  in  the  right  hand  side  of  above  equation,  (Z)  /  /?)  V^/7 ,  originates  from  the 
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Fokker-Planck  equation  governing  the  diffusion  of  objects,  Nelson  (1966).  By  using  an 
appropriate  algebraic  manipulation  of  Eq.  (2.3)  we  find  that  p!  p  is  equal  to 
/)“'[(! /2)(w^  +  V^)  + V^/  m^-E !  m^'\.  Upon  substituting  the  latter  expression  in  (3.4),  we  have, 

-  V.  V  =V.  Vln/7  +  [-  (m  "  +  V" )  +  ^  -  — ]  (3.5) 

2 

Thus,  physically,  the  specific  dilatation  -V»V  is  equal  to  the  rate  of  the  change  of  expansion  or 
contraction  along  the  direction  of  the  diffusion  velocity  plus  the  total  kinetic  energy,  including  the 
mean  convective  and  diffusive  fields,  the  specific  external  potential  energy  minus  the  specific  total 
energy,  E,  divided  by  the  diffiisivity. 

We  define  the  diffusion  field  as  a  “diffusive  dilatation  free”  when  -V»V  =0 .  In  contrast  to  the 
convective  flow  field,  the  dilatation  of  the  diffusion  velocity  field  in  a  hydrogen  atom  in  a 
stationary  state,  as  we  show  later,  is  not  zero,  even  in  atomic  bound  states.  It  is  the  dilatation  energy 
of  the  latter  “compressible  diffusive  field”  and  the  kinetic  energy  of  the  corresponding  diffusion 
process,  both  of  which  synergistically  acting  together  to  play  a  major  role  in  preserving  dynamic 
equilibria  of  probability  fluid  in  an  atom. 

Typical  dilatation  energy  of  a  diffusion  field  in  a  hydrogen  atom,  estimated  by  the  approximate 
expression  of  mfiV  la  is  10  ev,  is  of  the  same  order  of  magnitude  of  the  ground  state  energy  of  a 
hydrogen  atom. 

Next,  according  to  the  definition  of  diffusion  velocity,  (3.1),  the  term  defined  by 
(^^^/4m/)^j/^jv(ln/?)«V(ln/?) ,  represents  the  kinetic  energy  associated  with  the  quantum 
diffusive  velocity. 


Q  = 


r  ^ 

rn 

UJ 

V(ln/?).V(ln/7)  = 


r-i  ^ 


(3.6) 


The  characteristic  specific  kinetic  energy  of  diffusion  in  a  hydrogen  atom  is  estimated  to  be  10  ev., 
which  is  also  of  the  same  order  of  magnitude  of  the  ground  state  energy  of  hydrogen  atom. 
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Thus,  by  introducing  03T  and  Q,  we  formally  express  the  quantum  potential  O  as  the  sum  of 


the  specific  dilatation  energy  and  diffusion  kinetic  energy. 


n=an+Q= 


rn 

-Z)V*V  + 

v*v 

(3.7) 


We  conclude  that  the  quantum  potential,  given  by  Eq.  (3.7),  is  in  agreement  with  the  corresponding 
expression  given  by  Nelson  (1966). 

It  is  to  be  noted  that  the  equation  governing  the  diffusion  velocity  and  the  momentum  equation 
constitute  the  complete  set  of  conservation  equations  given  by  Nelson  (1966).  However,  the 
derivation  of  the  Nelson’s  diffusion  equation  can  be  achieved  through  the  use  of  the  continuity 
equation  (2.2)  and  diffusion  velocity  equation,  (3.1),  as  follows.  By  taking  a  partial  derivative  with 
respect  to  time  of  the  diffusion  velocity  Eq.  (3.1)  and  using  the  continuity  equation  to  replace 
dp/ dt  by  -V»( pu)  ,yNQ  arrive  at, 

^  =  (^)VV.u-V(u.V)  (3.8) 

dt  2m^ 

Thus,  Nelson’s  equation  governing  the  diffusion  velocity  is  identical  to  the  continuity  equation  of 
the  probability  density  derived  from  the  Hamilton-Jacobi  formalism,  together  with  the  definition 
of  the  diffusion  velocity  given  by  (3.1). 


4.  Principal  modes  and  elemental  processes 

By  adopting  the  definitions  of  quantum  diffusion  and  related  energies,  the  Bernoulli’s  equation,  Eq. 
(2.3),  for  the  probability  fluid  of  a  hydrogen-like  atom  can  be  expressed  as, 

^  +  -(u.u-V.V)  +  ^  +  /)V.V  =  —  (4.1) 

dt  2 

Bernoulli’s  equation  is  a  statement  of  the  balance  of  the  mechanical  energy  among  various  modes 
represented  by  two  types  of  flow  variables.  The  first  are  those  belonging  to  the  principal  modes: 
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quantized  energy,  Coulomb  potential,  the  kinetic  energy  of  the  mean  convective  motion,  and  the 
energy  associated  with  the  temporal  variation  of  the  velocity  potential.  The  next  are  those  which 
belong  to  the  quantum  potential  modes:  dilatation,  and  kinetic  energies  of  diffusion  induced 
potentials.  In  the  modal  analysis  method  to  be  presented  here,  we  decompose  a  quantum  potential 
mode  into  a  group  of  elemental  processes  each  of  which  has  either  a  part  of  the  dilatation  energy  or 
kinetic  energy  or  a  combination  of  the  both.  The  purpose  of  the  modal  analysis  is  to  determine  how 
each  principal  mode  balances,  energetically  or  dynamically,  with  an  appropriate  set  of  elemental 
processes.  The  successful  determination  of  such  sets  of  elemental  processes  provides  the  answer  to 
three  unsolved  issues  raised  at  the  outset  of  this  paper. 

(a)  Quantum  modal  balance  theory 


When  the  atom  is  in  a  stationary  state  Eq.  (4.1)  can  be  recast  in  spherical  polar  coordinates 
resulting  in. 


1 


-u 


Ze" 


M  0r 


2  mr  4m 

e  fc  V 

1  d\  1 

+  2  •  2q  2  +  T 

r  sm  6  ocp  2 


2  5 

r  — Inp 
dr 


1 


dlnp 

dr 


r  sin  0  00 
1  f  dlnp\ 


^  ^sin0— 

50 


+  - 


1 


r  sin  0 


J 

din  p^ 

0(p 


m„ 


(4.2) 


where  Z  (=1)  is  the  atomic  number,  e  the  charge  on  an  electron.  For  algebraic  convenience  the 
terms  of  Eq.  (4.2)  can  be  regrouped  and  written  as  follows 

J?(n,/,w)  +  S5-9Jt(n,/,m)-Q(n,/,m)  =  (£(«) .  (4.3) 

Here,  n,  I,  and  m  are  the  principle,  orbital  angular  momentum  and  magnetic  quantum  numbers, 
respectively,  for  an  atom  in  a  stationary  state.  The  terms  of  Eq.  (4.3),  which  are  all  per  unit  mass, 
are  sequentially:  J?[n,l,m)  is  the  kinetic  energy  of  mean  flow  induced  by  a  line  vortex;  X>  is  the 

Coulomb  mode,  which  is  independent  of  the  quantum  state;  is  the  dilatation  mode; 

Q(n,/,m)  is  the  diffusion  kinetic  energy  mode;  and  (£(n)  is  the  quantized  state  energy  mode. 

Formulation  of  the  quantum  modal  analysis  begins  with  the  definition  of  elemental  processes, 
described  below.  First,  the  dilatation  and  the  kinetic  energy  associated  with  the  diffusion  of  the 
probability  fluid  are  decomposed  into  the  radial,  polar  and  azimuthal  and  components,  as  follows. 
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where  the  subscripts  R,  0,  and  (p  designate  the  components  in  radial,  polar,  and  azimuthal 
directions,  respectively. 

Secondly,  each  component  of  dilatation  mode  («,/,m)  and  diffusion  kinetic  mode 

Q.Rg^[n,l,m)  is  further  divided  into  a  number  of  elemental  processes  («,/,m)  and 
as  follows. 


(n,  l,m)  =  Yj  .  («.  1,  m) , 

j 

(4.7) 

(n,  /,  m)  =  ^  {n,  /,  m) , 

(4.8) 

where  the  subscript  a  refers  to  the  R,  6,  or  (p ,  and  j  represents  the  j  th  elemental  process.  These 
processes,  and  ,  are  the  unit  processes,  each  of  which  contributes  in  such  a  manner  that 

the  sum  of  a  set  of  selected  elemental  processes  of  DJI  or  Q  modes,  or  combination  of  the  two, 
balances  with  a  specific  principal  mode:  X>  ,  or  (£(n)  .  Additionally,  there  are 

cross-intercoupling  elemental  processes,  which  make  no  contribution  in  energy  balance  with  a 
principal  mode.  They  do  serve,  however,  to  maintain  equilibrium  among  the  remaining  elemental 
processes,  which  are  not  participating  in  the  principal  modal  balance. 

(b)  Quantum  modal  balance  equations  and  scaling  laws 

Quantum  modal  balance  analysis  seeks  to  determine  appropriate  sets  of  elemental  processes 
such  that  each  set  satisfies  the  conditions  of  modal  balance,  given  below. 
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Coulomb  mode: 


-Zg" 

m^r 


Vortex  mode: 


2m/  sin^  9 


Quantized  state  energy  mode: 


(«,/,m)  +  2;Q®  (n,/,m) 

J  J 


-€(«) 


Z^g" 

Im^an^ 


in^tr' 


Cross  inter-coupling  mode: 

+  {n,l,m)  =  0 

j  j 


(4.9) 


(4.10) 


(4.11) 


(4.12) 


The  kinetic  energy,  J?(^n,l,m),  of  vortex  induced  flow  appears  in  the  right  hand  side  of  Eq.  (4.11) 
is  discussed  in  a  later  part  of  this  sub-section. 

In  order  to  carry  out  the  modal  analysis  it  is  convenient  to  first  expand  the  probability  density, 
p,  in  terms  of  hydrogen  eigenfunctions  expressed  in  spherical  polar  coordinates  as  follows, 
(Landau  &  Lifshitz  1985), 

Pnim  {r,  d,  (j))  =  ¥nlm¥l„n  =  K  (^’  K  (^’  ^ 

where  R„i  (r)  is  given  by, 

(4.14) 


in  which  P  =  2Z/na,  a  is  the  Bohr  radius,  dL  =  Ml  ,  and  ri  =  Pr  .  The  function  F„i,  is 

proportional  to  the  associated  Laguerre  polynomials,  and  is  governed  by  the  following  differential 
equation: 

9F 

;7-^  +  (2/  +  2-;7)-^  +  (^-/-1)F„,=0,  (4.15) 

arj  arj 
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where,  the  allowed  values  of  /  are  /  =  0,  1 ,  2,  3, ...  n- 1 . 

The  spherical  harmonics,  [9,(1))  can  be  separated  into  polar  and  azimuthal  angular 
components, 

=  (4.16) 

where  0,^  [9)  are  the  solutions  of  the  Legendre  equation. 


_l _ ^ 

sin6>  dO 


+ 

d0  ) 

/(/  +  1) 


m 


sin^  8 


0;«=O, 


(4.17) 


and  have  the  form. 


0.(^)=(-ir 


1(21 +  !)(/- 

m|)! 

I  2  (/  + 

m 

1)!' 

.p\^' 


(cos6>) , 


(4.18) 


where  m  takes  on  the  values  of  m  =  0,  ±1,  ±2,  ...  ±/,  and  if'''  are  the  associated  Legendre 
polynomials.  The  are  normalized  with  respect  to  integration  over  the  entire  solid  angle. 


In  order  to  calculate  the  kinetic  energy,  M[n,l,m),  of  vortex  induced  flow,  Eq.  (4.10),  we 
first  determine  the  velocity  of  the  mean  flow  from  Eq.  (2.4), 


u 


nlm 


Vln 

^  mH  '' 

y2mi  J 

V  ^ n/m  ) 

sin  8  ^ 

(4.19) 


The  above  velocity  distribution  corresponds  to  the  flow  induced  by  a  line  vortex  extending  in  the 
direction  of  polar  axis  coinciding  with  the  z-axis,  which  passes  through  the  center  of  the  mass  of  an 
atom.  The  kinetic  energy,  per  unit  mass,  of  the  flow  induced  by  the  vortex  mode  is  then. 


J?[n,l,m) 


rn  2  ^ 


(4.20) 


By  definition,  the  circulation  strength  of  the  vortex  flow  is  given  by. 


^  (•  ,  mh  r?.ind  ,  iTimh 

r  =  u»ar  =  —  — ^^—d(p  = - , 

•'  rsin^ 

e  e 


(4.21) 
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where  again  m  takes  on  the  values  of  0,  ±1 ,  ±2,  ...  ±1.  Thus,  the  vortex  in  the  atom  is  quantized 
with  the  magnetic  quantum  number,  m. 

We  will  now  show  that  the  Bernoulli’s  equation  (4.3)  can  be  cast  in  the  form  of  a  scaling  law, 
£  +  [n,l,m)-s\  =  0 ,  (4.22) 

(  a  j  a  j  j 


where  is  a  quantum  scaling  parameter,  and  /?,  S),  371,  and  s  are  the  non-dimensional 
quantities  corresponding  to  J?,  2) ,  071,  Q,  and  ^  of  Eq.  (4.3).  The  quantum  scaling  parameter  (E^ 


is  found  to  be  equal  to  four  times  of  the  negative  of  quantized  energy  per  unit  mass, 

=1 


4m 


mn  a 


(4.23) 


(c)  Decomposition  of  quantum  diffusion 

The  density  distribution,  Eq.  (4.13),  can  be  expressed  as  the  product  of  four  density  functions. 


where 


2Z 

p{n,l,m)  =  ( — f  pip^pppg 
na 


Pi='if\  Pa=^ 


(4.24) 


(4.25a,b) 


pF=FnM),  p,=&]„{0) 


(4.26a,b) 


where,  represents  the  power  type  distribution  associated  with  the  orbital  angular  momentum 
quantum  state,  is  an  asymptotic  distribution  for  large  r,  Pp,  is  the  density  distribution 

defined  by  Lagueere  polynomials,  and  Pg  is  the  density  distribution  due  to  Legendre  polynomial 

function.  By  using  the  definition  Eq.  (3.1),  the  diffusion  velocity  is  expressed  by  the  sum  of  four 
components. 


V  =  -DVln/7  =  V,  +Vj+V^+V^ 
=  (V,+V,+VJe,+V,e, 


(4.27) 
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where,  and  are  unit  vectors  and  the  four  diffusion  components,  %  ,  and  are 

sequentially:  power  type,  Lagueere,  asymptotic,  and  polar  diffusion  velocity,  respectively.  These 
are  given  by 


V,  =  -Z)Vln/7^  =  -DVIylt]^'  =  -{2/3D-) 

7 

(4.28) 

=  -DVln/?^  =  -Z)VlnFj  =  -2pD{—^)  e,. 

F„,  dt] 

(4.29) 

=  -DV\a.p^  =  -FVlne""  =  /3Dt^ 

(4.30) 

V,  =  -FVlnp,  =  -Z)V  In0^  =  -2/?D(^^)e, 

77©  dr] 

(4.31) 

where,  /3D=  Ze^ /nfi. 


The  dilatation  energy  of  diffusion  field  is  also  decomposed  into  the  corresponding  four  parts. 


m  =  -Dv.v  =  -Dv.v,  -  Dv^Vj,  -  -  Dv*y^ 


(4.32) 


where. 


=  -DV-\,  =  V,)  =  2/3^D\-^) 

r  or  Tj 


(4.33) 


=  -DV.V^  =  -~{r"  V^)  = 
r  or 


1 

d"F„, 

f  1 

Fn, 

djf 

Vf„,  dr]  ) 

riFni  dri 


(4.34) 


^,=-Dy^\,=-^^{r^W,)  =  2/3^D\k 

r  dr  Ti 


(4.35) 


m,= 


rsin^  dO 


2p^D^ 

1  d'^%  1 

(  1  d&^, 

2 

1  ^  coiO  d@ 

if 

©  d9'^  1 

t©  dO ) 

1 

© 

h 

2/3^D^ 


Tl 


/(/  +  !)- 


m 


sin^  0 


(4.36) 


Next,  the  kinetic  energy  of  diffusion  £2,  is  calculated  from  Eq.  (4.27  )  as  follows. 
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(4.37) 


a  =  iv-=i(v,  +  v,+v,).(v,  +  v,+v,)+i(v,.v,) 

2  2  2 

=  )<V,=  + v;  + V])  + V,  V,  + V,  V.  + V.  v,  +  iv,' 

2  2 


where,  are  radial  diffusion  components  and  \  is  the  polar  component  given  by  Eqs. 

(4.28)-(4.31),  respectively. 

(d)  Elemental  processes  of  dilatation  mode 

The  dilatation  mode  has  three  components,  radial,  azimuthal,  and  polar  modes,  each  of  which 
is  calculated  below.  The  radial  component  of  the  dilatation  mode  is  decomposed  into  three 

elemental  processes,  dJta.  {n,l,m),j=l,2,3  ,  as  follows. 


=  -z)v.v^  =  -z)v.v,  -z)v.v^  -z)v.v^ 

where,  dJt^  ,  971^  and  971^  are. 


(4.38) 


1  d^F„,  2  dF„,  1 


-  +  - 


4m/  ^  drj^  F„,  dp  p 


- )  =  -Z)V.V^-Z)V.V^+-V;  (4.39) 


4m^  F„,  dp  2 


(  4.40) 


= - y 


(2f^) 


^  I  ^ 


\n  j 


=  -Z)V.V, 


(4.41) 


The  dilatation  mode  in  -  direction,  971^  ,  is  decomposed  into  three  elemental  processes. 


971^  =971^  +971^  +971^  =-Z)V.V^ 


( 4.42) 


where,  971^^  for  i  =  1,  2  and  3  ,  is  given  by 


971,  = - r 

'  4m/ 


^  m^  ^ 

p^  sin^  6 


—  n,  +  n  , 


(4.43) 


971.  =■ 


4m„ 


1 

r  1 

d®,.] 

2 

7" 

1®/™ 

dO  j 

=  --V" 
2^’ 


( 4.44) 
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(4.45) 


where,  n„  =9Jl„+Q^. 

The  component  of  the  dilatation  mode  in  (fh  angular  direction  vanishes  and  makes  no 
contribution  in  the  energy  balance, 


9Jl,=0. 


( 4.46) 


(e)  Elemental  processes  of  diffusion  kinetic  mode 

The  diffusion  kinetic  mode  has  the  radial,  azimuthal,  and  polar  angular  components.  The  radial 
component  of  diffusion  kinetic  mode,  («,/,m),  can  be  broken  down  into  four  elemental 

processes,  Qri  for  /  =  1,  2,  3  and  4  respectively. 


+^fi3  + 


( 4.47) 


n  =^(2/3^)  -  1-^  -1  =(YV,  +  V,V,  +  V,V,) 

4m^  I  77j[F„,  dp  I  p 


( 4.48) 


/)2  I  1  jp  \  1 

Q  =Q  =^—(2^2^  J - =iv^ 

4m/'  ^  dp]  2  " 


(4.49) 


f  \  \ 


(4.50) 


Q  =Q  = -  {2f)-  =  -Wl 

"^4  ^  4^2V  ^^42^ 


(4.51) 


The  diffusion  kinetic  mode  in  direction  0  is  represented  by  a  single  component,  ’  given  below. 


Q  ^J!—{2/3f\\  ^ 

4mf'  ^  d0 


--v;. 

2  ^ 


(4.52) 


The  azimuthal  component  vanishes  and  thus  does  not  contribute  in  the  modal  balance. 
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a = 0 . 


(4.53) 


(f)  Kinematic  and  dynamic  characteristics  of probability  fluid 

The  probability  fluid  flow  in  the  hydrogen  atom  exhibits  interesting  flow  characteristics, 
discussed  below.  The  dilatation  of  the  convective  velocity  is  expressed  in  term  of  the  inner  product 
of  u  and  V,  defined  by  Eq.  (3.4),  and  is  given  by 

V.u  =  =  ue,.[(V,  + V^+ V,)e,.  +  V,  ej  =  0  (  4.54) 

Thus,  the  convective  velocity  fields  of  all  the  states  defined  by  a  set  of  n,  1,  and  m  of  a  hydrogen 
atom  are  convective  dilatation  free.  Alternatively,  we  write  the  following  identity, 

V.u  =  =  ^[(V  yrf  -iVipy]  (  4.55) 

D  2im^ 

If  the  state  of  an  atom  is  described  by  of  a  superposition  of  eigenfunctions  each  of  which  has  the 
different  values  of,  n,  and,  I,  but  the  same  value  of  m,  then  we  can  show  that  the  above  equation 
vanishes.  Accordingly  the  flow  is  convective  dilatation  free. 

In  a  convective  dilatation  free  flow,  the  velocity  scalar  potential  satisfies  Laplace  equation,  i.e., 

V.u  =  =  VV  =  0  (4.56) 

This  is  a  uniquely  interesting  property,  which  permits  the  construction  of  the  convective  field  by 
the  potential  flow  theory  of  classical  hydrodynamics,  whereas  the  density  is  determined  by 
Bernoulli’s  equation  in  conjunction  with  the  convective  velocity  field. 


5.  Atomic  modal  balance 

Modal  analysis  aims  to  determine  the  sets  of  elemental  processes  which  satisfy  the  quantum  modal 
balance  equations,  (4.9)  to  (4. 12),  and  to  ensure  that  the  sum  of  all  the  sets  of  elemental  processes 
is  equal  to  the  total  quantum  potential.  The  procedure,  in  general,  involves  the  trial  and  error 
method  to  properly  choose  prospective  elemental  processes  according  to  the  nature  of  the  radial 
decay,  i.e.,  f^,  for  s  =  0,  1  and  2,  as  described  in  the  modal  balance  equations  (4.9)  to  (4.12). 
However,  the  elemental  processes,  and  ,  involve  functions  such  as  ,  (F^,)”' 
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{d^F^i  / dr}^)  ,  (F^i)  \dF^,  / drj) ,  and  (F^,)  ^ (dF^,  / drjf  .  Since  these  latter  terms  reveal  no 

explicit  functional  form  of  radial  decay,  i.e.,  r\  the  selection  of  appropriate  elemental  processes 
requires  knowledge  of  the  functional  characteristics  of  these  special  functions  in  order  to  facilitate 
the  correct  identification  of  the  prospective  elemental  processes.  While  the  trial  and  error  method 
is  in  general  carried  out  by  inspection,  the  use  of  the  momentum  equation  in  conjunction  it  offers  a 
physical  guide  to  the  selection  of  elemental  processes. 

Detailed  energetic  balance  of  each  principal  mode  and  cross  inter-coupling  modes  are  discussed  in 
the  following  sub-sections.  The  results  of  the  algebraic  analysis  of  the  modal  balance  are 
summarized  in  Table  1 . 


(a)  Coulomb  modal  balance : M  ri, 


The  Coulomb  potential,  W^/m^,  which  is  a  central  force  field,  plays  a  major  role  in  radial 
modal  balance.  We  select  the  following  two  elemental  processes  and  for  the  Coulomb 
modal  balance. 


>,2  r  2yj2A 


4m^ 


\Pnldj 


T7 


d^F  dF 

+(2/  +  2-7)^^-(/  +  1)F„ 


drj 


dij 


nl 


=  .OV-V,  -  DV. V,  +  (V,  V,  +  V,  V, )  +  ^ 


2 


(5.1) 


=  n,-n,--vj 


where,  ^D  =  Ze^ !  nh  ■,  andfl^  is  the  centrifugal  potential  expressed  by 

(5.2) 

2m  r 

In  order  to  verify  that  the  sum  of  the  above  two  terms,  and  ,  is  equal  to  the  Coulomb 
potential,  we  multiply  Eq.  (4.15)  by  {h^ !  Am !  F^,ri)  and  subtract  the  resulting  equation 
from  Eq.  (5.1)  to  prove  that  the  sum  of  and  indeed  add  up  to  Coulomb  potential,  as 
shown  below. 
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\KiVj 

^2^ 

\Pninj 

^2^ 

kKiTIj 


J! 


(fF  dF 

- ^  +  (2/  +  2-7)^-(/  +  1)F„, 


dri^ 

^d^F,  ^ 


drj 


K/)= 


dtl  j 

-h 


+  {21  +  2-t]) 

2  f  ^  n2  \ 


ydrt  j 


+  («-/-l)F„, 


(5.3) 


Am 


.  \  d  j 


Ze  V 


n  =  — 


mr  m 


where  t]  =  jdr,  and  /3  =  2Z  !  na,  were  defined  in  Eq.  (4.14).  The  above  energy  balance  also  gives 
the  following  picture  of  momentum  equilibrium.  Observing  that  =  Ze^  /  nfi  is  constant,  taking 

the  partial  derivative  of  Eq.  (  5.3  ),  with  respect  to  r  gives  the  dynamic  equilibrium  of  Coulomb 
with  the  quantum  force, 

8(!W,,+n.,)  g(n.-n,-(v;)  a(n,-n,)  ^  i 

dr  dr  dr  dr 


The  above  equation  states  that  the  radial  quantum  potential  force,  0n^  /  dr  minus  the  centrifugal 
force,  dHj  /  dr ,  balances  with  the  Coulomb  force  for  an  atom  in  a  stationary  state.  Significantly,  it 
must  be  pointed  out  that  contrary  to  previous  claims,  Bohm  and  Hiley,  (1999),  for  example,  which 
state  that  the  total  quantum  force,  (5n^+Tlg)/0r  ,  is  what  balances  the  Coulomb  force, 

(1  /  m^)(d  Vj  dr) .  The  latter  is  inconsistent  with  the  present  criterion  of  radial  dynamic  equilibrium. 


From  the  view  point  of  dynamic  balance,  we  are  interested  to  know  the  relative  importance  of 
the  dilatation,  971^  and  the  kinetic  energy  of  diffusion  in  the  Coulomb  modal  balance.  The 

fractions  of  the  quantum  potential  contributing  to  the  Coulomb  modal  balance  may  remain 
constant  or  exhibit  spatial  variations  for  the  selected  quantum  numbers  n,  I,  and  m  as  shown  in 
Table  2.  The  relative  contribution  of  and  971^^  on  the  Coulomb  modal  balance  as  function  of 

radial  distance  is  given  by. 


r 


971. 


V 


J_^ 

Fni  drj 


)- 


I 

n 


J_^ 

Ki  dyf 


2^ 

1 

dFA 

1 

UJ 

\^nl 

dt]  J 

d 

(5.5) 
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The  asymptotic  values  of  and  its  derivatives,  which  appear  in  (5.5)  are, 


1  dF^, 


r  1  A 


lim- 

F^i  drj 


+  0 


11 


J 


1  d^F^  ^  ^ 


lim 

F^i  drj^ 


iti  _ 


O 


y 


(5.6) 

(5.7) 


Using  these  asymptotic  forms  in  Eq.  (5.5)  y  becomes, 

limy^n-1.  (5.8) 

rj^cc 

Thus,  the  kinetic  energy  of  diffusion  dominates  over  the  dilation  energy  at  higher  principal 
quantum  number  in  the  Coulomb  modal  balance.  Numerical  analysis,  detailing  the  Coulomb 
modal  balance  and  its  dependence  on  the  quantum  numbers  n,  /and  m,  is  presented  in  §6. 


The  ratio  of  the  contribution  by  the  dilatation  energy  in  Coulomb  modal  balance  to  the 
total  radial  dilatation  energy,  937^ ,  at  a  non-dimensional  radial  distance  is  given  by 


\  I 


F„,  dif 


-  + 


\n) 


f 


IdF^ 

\F„i  dri 


f-l 

UJ 


1  d"F, 

K,  dif 


-  + 


'7^ 


f 


1  dF, 

_ n]_ 

\K,  dt] 


\Ki  dt] 


1  / 
+  — 
11  11 


(5.9) 


By  substituting  asymptotic  forms  given  by  Eqs.  (  5.6)  and  (5.7)  in  Eq.  (5.9)  we  find  the  limiting 
value  of  to  be. 


limcTg^^l.  (5.10) 

77^00 

This  result  implies  that  the  major  fraction  of  dilatation  energy  contributes  to  the  Coulomb  potential 
balance  at  large  radial  distances. 

Similarly,  the  fraction  of  diffusion  kinetic  energy  contributing  in  modal  balance  is 
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(5.11) 


L  V  Kl  dt]  q\ 

(  1 

2 

1 

fl 

ri  ) 

1  + 

fl 

ri  ) 

2 

1 - 

[Ki  dri  j 

The  asymptotic  value  of  cr^  vanishes,  i.e.,  lim  ^  0  .  Hence  the  fraction  of  the  diffusion  kinetic 

r/^cc 

energy  contributing  in  the  Coulomb  potential  diminishes  at  large  radial  distances. 

(b)  Vortex  modal  balance:  Hg  and  11^ 


The  vortex  induced  flow,  in  the  azimuthal  angular  direction  for  a  state  with  magnetic 
quantum  number  m,  is  driven  by  quantum  potentials  fl^  and  H  (=  fl^  ) ,  as  follows, 


1  %  m 

rig+n  =971^  +  +  971^3  +  0^3  =  ^(n,l,m)  =  -ul,^  =- — y  2  ■  in 

'  "  2  2m^  r  sm  0 

=  -DV. ^ \ 


(5.12) 


The  appearance  of  the  radial  potentials  n,(=  Tl^  )  in  the  kinetic  energy  of  the  azimuthal  motion, 
in  Eq.  (5.12),  suggests  the  presence  of  cross  inter-coupling  in  the  modal  balance  of 
n^(^,r)-l-n;(r) ,  and  thus  requires  special  attention.  Since  cross  inter-coupling  also  occurs  in 
radial  direction,  we  shall  discuss  these  aspects  together  in  subsection  (d). 


(c)  State  quantized  energy  modal  balance:  Q 


The  elemental  process  that  contributes  to  the  quantized  energy  is  equal  to  the  negative  of  a  part 

1  2 

of  radial  diffusion  kinetic  energy,  ,  as  shown  by  the  following  equation. 


= - 7 

“  4m  ^ 
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8m„ 
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z  e 


2m  an  2fi  n 


2_2 
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m 


(5.13) 


The  asymptotic  diffusion  velocity  is  a  constant,  =  Ze^  /  nfi ,  and  thus  the  quantized  energy  has  a 

constant  value.  We  shall  discuss  the  physical  explanation  of  the  term  =l/2Vj  as  the 
quantized  energy  in  subsection  (f). 


23 


Another  element  of  the  state  quantized  energy  is  the  non-dimensional  energy  scaling  parameter 
£.  Its  value  can  be  deduced  from  Eq.s  (4.22)  and  (4.23)  as  follows, 


rs  ryT.  2 

2Z  e 


AE  £ 


= - 

m  Am„  m  n  a  m„ 


(5.14) 


The  above  relation  yields  ^ 


The  fraction  of  diffusion  kinetic  energy  contributing  to  the  quantized  energy  is, 


p  _ 


Q, 


'■R 


J_^ 
\Pnl  dV 
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\Fni  dri 


ri ) 


+ 


v4y 


rt ) 


(5.15) 


Using  the  asymptotic  expressions  given  by  Eqs.  (5.6)  and  (5.7),  the  above  equation  reduces  to, 

\\m£^^=\.  (5.16) 

Jj^cc 

Thus,  the  majority  of  diffusion  kinetic  energy  at  large  radial  distance  contributes  to  the  quantized 
energy.  The  spatial  variation  of  the  contribution  by  diffusion  kinetic  energy  on  quantized  energy  at 
quantum  numbers,  n,  I,  and  m  is  presented  in  §6, 

(d)  Cross  inter-coupling  modal  balance 


Three  cross  inter-coupling  modes  occur  among  the  components  of  quantum  potentials  in  the 
hydrogen  atom  are  identified  below. 


The  expression  for  dilatation  energy,  given  by  Eq.  (3.5),  can  be  recast  in  the  Bernoulli’s 
equation, 


-Dv.v + - V  =  n. + n, + n,  =  (»i, + n,)+ (ffli, + a) + +  q.) 


(5.17) 


The  above  equation  can  be  regrouped  as  the  sum  of  the  following  six  sets  of  potential/energy 
terms, 
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(5.18) 


(971^,  +  - -)  +  (97lg,  -K)  +  (0^4  +  — )  +  (971^2  +  ^?2 )  + 

(971^2  +^?2)  +  (^(?3  +^^3)  =  0 


Since  the  first  three  sets  describe  principal  modal  balance,  as  shown  by  Eqs.,  (5.3)  (5.12)  and 
(5.13),  hence  each  set  vanishes.  The  remaining  three  sets  do  not  contribute  to  the  principal  modal 
balance,  hence  they  represent  the  cross  inter-coupling  modal  balance,  as  detailed  below. 

(1)  Cross  inter-coupling  modes  in  vortex  mode 

We  identify  two  sets  of  cross  inter-coupling  modes  in  Eq.  (5.18),  as  follows. 


®>„,+£\=-tv„>+ivJ  =  o, 


(5.19) 


^0^  +  n«3  =  971^^  +  971^3  +  Q^3  =  97lg^  -  DV»\,  +  -  V/ 
4m„  r^ 


(5.20) 


By  substituting  Eqs.  (5.19)  and  (5.20)  into  (5.12),  we  find  the  kinetic  energy  consists  of  the  first 
part  of  the  dilatation  in  polar  angular  direction,  971^  ,  as 


^{n,l,m)  =  dJlg  = 


fd 


(2^^) 


m 


rf  sin  9 


(5.21) 


Observe  that  since  971^  »  the  two  modes  balance  each  other  in  Eq.  (5.19),  hence  they 

constitute  cross  inter-coupling  modes.  Equation  (5.20)  expresses  the  cross-intercoupling  between 
the  centrifugal  potential,  n  =  (fl^j),  with  dilatation  energy,  971^  ,  in  the  polar  angular  direction. 

It  is  interesting  to  observe  that  the  ratio  of  the  radial  diffusion  kinetic  energy,  ,  to  the  radial 
dilatation  971,  in  the  centrifugal  potential  n  is  equal  to  the  orbital  angular  momentum  quantum 

number,  /,  hence  the  effects  of  the  kinetic  energy  of  diffusion  become  more  pronounced  than  those 
of  the  dilatation  energy  in  their  prospective  contribution  to  the  centrifugal  potential  at  large  value 
of  /. 
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(2)  Cross  inter-coupling  modal  balance  in  radial  mode 


We  note  that  the  dilatational  energy,  ,  Eq.  (4.40)  and  the  kinetic  energy  of  diffusion,  , 

(4.49)  are  equal  but  have  an  opposite  sign,  hence  both  processes  are  cross  inter-coupling  modes, 

M«.+^>.=-|V+5V/=0  (5.22) 


(e)  Global  atomic  modal  balance 

Global  modal  balance  among  all  the  principal  modes  and  the  quantum  potentials  can  be 
validated  by  demonstrating  that  the  sum  of  all  the  elemental  processes  contributing  to  the  three 
principal  modes,  Eq.  (5.3)  (5.12)  and  (5.13),  is  equal  to  the  total  quantum  potential.  This  is  shown 
below. 


^+^-^  =  (n,-n„-ivj)+(n,+n„)+iv;- 

mm  z  z 

e  e 

= 

Here  we  included,  fl^  =  0 ,  in  the  quantum  potential  for  the  completeness. 


We  summarize  the  global  atomic  modal  balance  of  three  major  types  of  atomic  states  in  Table  3. 
Table  3.  Global  atomic  modal  balance  in  hydrogen  atom 


States 

V 

e 

m 

e 

m 

e 

m  m 

e  e 

n^0,l  =  0,m  =  0 

0 

-V" 

2 

total 

n  ^  OJ  ^  O^m  —  0 

n.-n„-tvj 

+n^3(=  0) 

-V" 

2 

total 

n  ^  OJ  ^  O^m  ^  0 

+n^3(?i  0) 

-V" 

2 

total 

Table  3  clarifies  three  issues,  which  are  raised  at  the  outset  of  this  paper:  (1)  Quantized  energy  is 
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equal  to  the  negative  of  the  kinetic  energy  of  the  asymptotic  diffusion  field  in  an  atom  at  a 
stationary  state;  (2)  Coulomb  potential  force  on  electron  probability  density  fluid  is  balanced  by 
the  total  radial  potential  force  minus  the  centrifugal  force.  (3)  Kinetic  energy  of  vortex  induced 
flow  is  provided  by  the  sum  of  the  components  of  quantum  potentials  in  the  polar  angular  direction 
and  the  centrifugal  barrier  potential  fl,  . 

(f)  Physical  significance  of  the  quantized  energy 

A  rather  remarkable  result  is  that  the  negative  of  kinetic  energy  of  asymptotic  diffusion, 
=  1/2  Vj ,  is  equal  to  the  quantized  energy,  as  given  by  Eq.  (5.13).  Moreover,  its  value  remains 

constant  with  respect  to  radial  coordinate.  These  two  facts  merit  further  elucidation.  First  note  that 
the  asymptotic  density  distribution  p  a  has  a  log  value  of, 

=  (5.24) 

Since  this  is  a  linear  function  of  the  radial  coordinate,  r,  multiplied  by  a  constant  coefficient 
f  ,  hence  the  asymptotic  diffusion  velocity,  ,  is  a  constant  equal  to  D/5  -Z] !  nh  .  The 

physical  significance  of  the  asymptotic  diffusion  velocity  is  elucidated  by  the  observation  that  the 
density  distribution  of  probability  fluid  is  governed  by  the  Coulomb  force,  which  attracts  the  fluid 
inward  on  one  hand,  and  the  quantum  potential  force  that  disperses  the  fluid  outward  on  the  other, 
see  figure  1.  Indeed,  we  show  that  the  asymptotic  density  distribution  depends  on  the  relative 
preponderance  of  the  absolute  value  of  the  characteristic  Coulomb  potential  Vj  and  quantum 

potential,  fl*  ,both  measured  at  the  Bohr  radius  Tl* ,  as  follows, 

8  4/i^  r  8  V*  r 

Pa  =  exp(-/?r)  =  exp( — )( - j){-)  =  exp( — )(— ^)(-)  (5.25) 

n  m^a  a  n  mJA  a 

where  V*  =  {ze^  /  a)  and  fl*  =  4h^  /  respectively.  The  ratio  of  the  two  characteristic 
potentials  W*J  mfl*  is  equal  to  z/4  .  For  a  hydrogen  atom  z  =  l,  thus  the  ratio  of  two 
characteristic  potentials  is  1  to  4. 

It  is  noteworthy  to  mention  that  the  linear  dependence  of  the  logarithm  of  the  asymptotic 
density  with  respect  to  the  radial  coordinate,  as  shown  above,  gives  a  constant  diffusion  velocity, 
Va,  which  in  turn  yields  a  constant  kinetic  energy  of  diffusion,  given  by. 
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(5.26) 


We  have  thus  established  the  following  dynamic  energetic  relationship  between  the  quantized 
energy  with  the  kinetic  energy  of  the  asymptotic  diffusion  process: 

^  =  -Q  =--V/  (5.27) 

*4  2 

The  fact  that  the  kinetic  energy  of  gradient  induced  asymptotic  diffusion  velocity  field  is  equal  to 
the  negative  of  the  quantized  energy  bears  some  generality  in  quantum  systems  with  a  central  force 
field  with  the  asymptotic  potential  for  large  radial  distances,  which  behaves  like  r  ^  For  example, 
our  recent  studies  show  that  in  a  diatomic  molecule  with  a  inter-atomic  potential,  U,  described  by 
the  central  field,  A/r^  -  B/r ,  the  negative  of  kinetic  energy  of  the  asymptotic  diffusion  is  also  equal 

to  the  quantized  energy.  On  the  other  hand,  however,  we  have  found  that  the  zero-point  energy  in  a 

three  dimensional  harmonic  oscillator,  which  has  a  potential  kr^ ,  is  equal  to  the  dilatation  energy 
of  the  asymptotic  diffusion  velocity  field. 


6.  Fluid  mechanical  and  electronic  structure  of  a  hydrogen-like  atom 

In  this  section  we  present  the  results  of  numerical  analysis  of  the  quantum  fluid  mechanical 
characteristics  of  the  modal  balance,  quantized  energy  and  the  structures  of  a  hydrogen  atom  at 
selected  states.  The  relative  contribution  of  each  elemental  process  on  modal  balance  is  also 
examined. 

(a)  Coulomb  modal  balance,  fluid  mechanical  and  electronic  structures 

The  general  energetic  characteristics  of  Coulomb  modal  balance  for  selected  quantum 
numbers  are  presented  in  figures  2~6.  The  Coulomb  modal  balance  of  Eq.  (5.3)  can  be  recast  into 
the  form,  in  which  the  Coulomb  potential  is  normalized, 

27  —  ~  1 

— (fm«.+QO  =  -T  ’  (6-1) 

n  ^ 
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r  77  —  ^ 

where  ^  =  —  = - - — .  The  analytical  expressions  of  the  spatial  distributions  of  OJIr,  and  Qs, , 

a  i2Zln) 

and  the  relative  contribution  of  two  processes  / OJIrj  for  selected  states  are  summarized  in 

Table  2. 

(1)  Fractional  and  complementary  modal  balances. 

Close  examination  of  Table  2  and  figures  2-6  suggests  that  there  are,  in  general,  two  types  of 
Coulomb  modal  balance:  fractional  and  complementary  balances.  In  the  former  case,  each 

elemental  process,  i.e.,  OJIr,  and  Qs, ,  assumes  a  fraction  of  the  Coulomb  potential.  Remarkably, 

the  latter  fraction  turns  out  to  be  j'l)  =  \l  n  ,  Q«,  l'D  =  {n-\)l  n  ,  and 

y  =  =n-\  throughout  the  space  of  range  0  <r  <co.  The  fractional  balance  occurs  in 

those  states  wherein  the  orbital  quantum  number  meets  the  following  condition, 

/  =  n  -  1  for  all  values  of  n  and  m  (6.2) 

For  example,  the  states,  \nlm)=  |100),  |210),  |21±l),  |320),  |32±l),and  |32±2),  exhibit  the 

fractional  balance  as  shown  in  Table  2  and  selected  numerical  examples  are  given  in  figures  2,  4 
and  6. 

From  the  expression,  y  =  n-l,  we  conclude  that  the  dilatation  plays  a  major  role  in  the  Coulomb 
modal  balance  for  the  ground  state,  whereas  the  diffusion  kinetic  energy  becomes  the  dominating 
balancing  mechanism  for  n>2,  see  figure  7.  The  above  results,  y  =  Qs,  /OTr^  =n-\,  and  the  ratio, 

Q,  /  dJt,  =  Qr^  =  I ,  discussed  in  section  5(d),  suggest  that  at  higher  quantum  numbers  n  and  /, 

the  diffusion  kinetic  energy  plays  a  greater  role  in  dynamic  and  energetic  atomic  modal  balance 
than  dilatation  mode. 

In  a  complementary  balance  mode,  the  energy  of  each  elemental  process,  i.e.,  Mrj  and  Qs, ,  is 
not  a  constant  fraction  of  the  Coulomb  potential  but  varies  in  space  as  shown  in  Table  2  and  figures 
3,  and  5.  The  spatial  variation  in  yalso  has  singular  behavior  at  the  location  where  OJIri  vanishes, 
as  shown  in  figures  8.  However,  at  a  large  radial  distances  where  r/a  — >■  oo ,  OJIr^  I'D  —>■  \/n,  and 
Q.rJ'D  —>  {n-\)ln.  These  ratios  are  the  same  as  those  of  the  fractional  balance  mode  described 
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above.  The  complementary  balance  pattern  occurs  in  the  states  where  the  condition  of  Eq.  (6.2)  is 
not  met,  for  example,  in  the  states  |200),  |300),  and  |310)  and  two  examples  are  shown  in 
figures  3  and  5. 

(2)  Singularities  in  complementary  balance  distributions. 

The  spatial  distributions  of  the  energies  of  an  elemental  process  exhibit  interesting  behaviors  at 
the  nodal  singularities  where  the  dilatation  and  diffusion  kinetic  energies  diverge.  Since  the 
Laguerre  polynomial  F„i  has  the  maximum  of  n-\  singularities,  there  are  corresponding  numbers 
of  singularities  where  the  divergence  phenomena  occur.  For  example,  the  state  1 200) ,  n=2,  has  a 

singularity  at  2  as  shown  in  figure  3,  whereas  the  state  |300)  has  two  singularities  at  ^~2, 

and  ^~7,  shown  in  figure  5.  The  quantities  and  Qr,  in  the  1 200)  state  exhibit  singularities  at 

4^=  2,  which  are  shown  in  figure  3,  and  will  be  discussed.  The  probability  fluid  diffuses  outward 
from  the  interior  region,  i.e.,  <4  ~<  2,  where  R20  >  0  .  The  diffusion  velocity,  defined  in  Eq.  (3.1), 
increases  as  the  probability  fluid  approaches  the  singularity.  The  expanding  flow  is  accompanied 
by  a  large  value  of  negative  dilatation  energy.  The  density  reduction  continues  toward  the  singular 
point  where  the  diffusion  velocity  becomes  infinite  and  the  density  reaches  zero.  The  diffusion 
kinetic  energy  diverges  to  positive  infinity,  whereas  the  dilatation  approaches  negative  infinity  at 
(4=2.  However,  as  ^  slightly  exceeds  2,  say  ^2'^,  the  density  of  probability  fluid  increases  in 
outward  direction,  and  hence  the  fluid  diffuses  inward  to  the  nodal  point.  At  <4  =  2^  the  diffusion 
velocity  approaches  negative  infinity,  and  thus  the  kinetic  energy  diverges.  This  abrupt  change  in 
the  direction  and  magnitude  of  diffusion  velocity  from  positive  infinity  to  negative  infinity  is 
anticipated  because  of  the  vanishing  inertia  of  the  fluid.  Likewise,  a  sudden  change  in  the 
dilatation  from  negative  infinity  to  a  finite  negative  value  is  achieved  at  ^=2^.  A  large 
deceleration,  experienced  at  4  =  2  yields  a  positive  dilatation  energy.  Hence,  upon  crossing  the 
4^2  singular  point,  the  dilatation  energy,  OJIrj  ,  increases  rapidly  to  a  finite  but  negative  value.  The 

same  divergence  phenomena  are  observed  at  other  singularity  points  for  the  states  1 300) ,  shown  in 
figure  5. 

It  is  interesting  to  note  that  although  the  specific  kinetic  energy  of  diffusion  and  the  dilatation 
energy  per  unit  mass  possess  divergent  singularities  their  corresponding  energies  per  unit  volume. 
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and  ,  do  not  since  the  probability  density  vanishes  more  rapidly  than  the  energies 
diverge  at  the  singularity.  That  is, 


),„g  - 

^1 


nl 


dr] 


2l\ 

1-  — 

7  j 


dr]  7] 

(6.3) 

- — F  ^  =  0 

2  ^nl  ^  ’ 

(6.4) 

where  the  subscript  “sing”  implies  an  evaluation  at  the  singularity  point. 

(3)  Fractional  contribution  of  the  energy  of  elemental  processes. 

The  fractions  of  the  dilatation  energy  agjj=  DJts, and  diffusion  kinetic  energy 

<7^  =Q.rJq.r  ,  for  selected  states,  |100),  |210),  and  |300),  are  plotted  in  figures  2,  4,  and  9, 
respectively.  In  general,  the  fractional  contributions  are  spatially  dependent,  with  the  exception  of 
the  special  case,  1 200)  ,  in  which  OJls,  / OJts  and  Qs;  / remain  constant  throughout  the  space. 

The  distributions  of  =Q.r^Iq.r  for  the  states  |100),  |210)  and  |300)  are  shown  in  figures 
2,  4,  and  10,  respectively.  Note  that  since  the  asymptotic  value  of  =  Qs,  at  a  large  radial 
distances  is  unity,  see  Eq.  (5.16),  all  the  dilatation  energy  contributes  to  the  quantized  energy. 

(b)  Vortex  modal  balance  and  core  structure 

When  the  magnetic  quantum  number  m  is  associated  with  a  bound  state  a  line  vortex  is  present 
in  the  core  of  the  atom  and  extends  from  -  co  to  co  along  the  z-axis  as  shown  in  figures  11,12  and 
13.  The  flow  field  induced  by  a  vortex,  which  is  superimposed  on  the  Coulomb  field,  makes  the 
atomic  core  structure  rather  complex  in  the  vicinity  of  the  nucleus.  A  comparison  of  the  relative 
magnitude  of,  and  2) ,  suggests  that  there  are  two  different  types  of  core  structures  depending  on 

the  magnetic  quantum  number.  The  first  type  is  the  Coulomb  potential  dominated  electronic 
structure  that  occurs  under  any  principal  quantum  number  n ,  and  orbital  angular  momentum 
quantum  number  /,  but  only  when  the  magnetic  quantum  number  is  zero,  w  =  0 .  The  second  type 
is  the  composite  core  structure,  which  consists  of  a  region  dominated  by  the  vortex  induced  core 
flow  and  the  Coulomb  potential  dominated  outer  region  as  shown  in  figures  11,  12,  and  13.  The 
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latter  structure  occurs  at  bound  state  values  of  I,  n,  and  m.  Two  regions  of  the  composite  structure 
are  separated  by  a  critical  radius,  ,  which  is  defined  by  the  equality  of  ^  with  D .  The  critical 

radius,  shown  in  figure  1 1 ,  is  given  explicitly  by 


a  1 

2Z  sin"  0  ■ 


(6.5) 


The  kinetic  energy  of  the  vortex-induced  flow  is  greater  than  the  Coulomb  potential  in  the  interior 
region  r  <rc,  whereas  the  opposite  is  true  in  the  exterior  region  r  >rc  .  The  critical  radius  depends 
on  the  magnetic  quantum  number  m  and  is  independent  of  n  and  I .  We  compare  the  critical  radius 
with  the  expected  value  of  the  radial  distance  of  the  electron  from  the  nucleus,  given  by  , 

('■>„/  =  jfl  sin  edOdcpdr  =  +  (6.6) 


In  order  to  examine  the  effects  of  the  vortex-induced  flow  on  the  value  of  we  consider  the 
limiting  state,  characterized  by  the  largest  magnetic  quantum  number,  at  a  principal  quantum 
number  n.  Let  /  =  /max  =  n-\,  and  \m\  =  wtmax  =  /max  =  n-\ .  The  ratio  of  to  is  , 

C  _  sin"  6  _  /g  yx 

{r)^j  [3n" -«(«-!)]  (2«"+«)sin"^ 

This  ratio  is  plotted  in  figure  11  for  the  state  |433).  The  region  marked  by  a  darker  color  is  the 
vortex  dominated  inner  core  whereas  the  exterior  region,  in  a  lighter  color,  is  Coulomb  potential 
dominated.  The  size  of  the  inner  vortex  core  is  0.234  (r^  43.  The  core  size  increases  with  an 
increase  in  the  principal  quantum  number,  which  is  seen  in  figure  12.  It  reaches  a  maximum  value 
of  O.Sir^  n,n-i  when  co,  at  the  state  of  ionization.  Thus,  the  inner  core  of  high  magnetic 

quantum  number  is  circulation  flow  dominated  in  contrast  to  a  stagnant  inner  Coulomb  core  for  the 
case  of  m  =  0.  However,  it  is  noted  that  the  structure  described  above  is  not  valid  within  the 
spherical  region  of  the  radius  of  (me/mp)re,  centered  at  the  origin,  figure  13,  because  of  the 
approximation  made  that  the  nucleus  is  located  at  the  center-of-mass,  which  coincides  with  the 
origin  of  the  coordinate,  as  discussed  previously. 
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7.  Conclusions 


Based  on  the  quantum  diffusive  fluid  dynamics  formalism,  we  examined  the  quantum  mechanical 
significance  of  quantized  energy,  and  dynamic  equilibria  of  an  electron  among  the  Coulomb  force 
field,  the  vortex  induced  flow  and  quantum  potential  in  a  hydrogen  atom.  In  our  model,  a  hydrogen 
atom  in  a  stationary  state  is  envisioned  as  an  electron  probability  fluid  attracted  by  the  Coulomb 
potential  field  induced  by  the  positively  charged  nucleus  located  at  the  origin  coinciding  with  the 
center  of  the  mass  of  two  particles.  Fundamental  characteristics  of  quantum  potential,  the  key 
internal  potential  for  the  energetic  and  dynamic  equilibria  in  an  atom,  are  examined  through  the 
dependence  of  the  potential  on  the  diffusion  velocity  field  created  by  the  non-uniform  spatial 
distribution  of  the  probability  density  in  a  stationary  atomic  state.  The  quantum  potential  field  is 
classified  on  the  basis  of  the  four  different  types  of  density  functions:  the  asymptotic  diffusion 

due  to  the  asymptotic  density  distribution;  Pp  ,  the  Laguerre  diffusion  induced  by  the  Laguerre 
function  type  distribution;  p^ ,  the  power  law  diffusion  due  to  the  orbital  angular  distribution 
function  ;  and  Pg ,  the  polar  angular  diffusion  induced  by  density  distribution  in  polar  angular 

direction.  Each  diffusion  field  exhibits  a  non-uniform  spatial  distribution  that  creates  the  quantum 
potential.  The  latter  potential  consists  of  the  dilatation  and  kinetic  energies,  whose  spatial  gradient 
yields  a  quantum  force.  The  global  energetic  bearing  diffusion  process  of  the  electron  probability 
fluid  is  comprised  of  the  sum  of  elemental  processes.  The  quantum  energetics  of  these  elemental 
processes  plays  a  major  role  in  the  modal  balance  between  the  quantum  potential  with  each 
principal  mode:  quantized  energy.  Coulomb  potential  mode,  and  vortex  induced  flow  mode. 

Modal  balance  analysis  predicted  that  the  quantized  energy  is  equal  to  the  negative  of  the 
kinetic  energy  of  asymptotic  diffusion  process  of  the  probability  fluid  and  is  given  by  En  /me,  = 
-  (1/2)  \ K  .  The  analysis  also  gives  a  fundamental  physical  explanation  of  why  the  energy  of 
eigenstate  has  a  constant  value,  as  summarized  below.  Since  the  logarithm  of  the  asymptotic 
density  InpAn  =  -  ^  is  a  linear  function  of  the  radial  coordinate,  the  asymptotic  diffusion  velocity , 

Va  =  -DdlnpAn/dr  =  (ze^/n  h) ,  takes  on  a  constant  value  for  each  quantum  state  n.  Thus,  based  on 
the  result  of  modal  analysis,  E„/me  =  -  (1/2)  Va^,  the  quantized  energy  has  a  constant  value,  i.e.,  E„ 
/me  =  -  (1/2)  (ze^/n  hf  .  It  is  noteworthy  to  point  out  that  the  classical  law  of  equipartition  of 
energy  ceases  to  apply  in  an  atomic  energy  balance.  We  found  that  not  all  the  kinetic  energies  of 
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diffusion  contribute  equally  into  quantized  energy.  For  example,  the  kinetic  energy  of  diffusion  in 
polar  angular  direction  does  not  add  up  to  the  quantized  energy  because  the  former  energy  is 
balanced  by  an  equal  part  of  the  dilatation  energy.  Furthermore,  the  kinetic  energy  of  diffusion  in 
azimuthal  angular  direction  is  zero,  and  thus  yields  no  contribution  to  the  quantized  energy.  Modal 
analysis  reveals  that  the  Coulomb  attractive  force,  Fc,  is  balanced  with  an  equal  but  opposite 
quantum  force,  Fq  =  F^-  F^  =  5(Ft^  -n,)/0r ,  to  preserve  radial  equilibrium,  and  thus  prevents 

the  electron  from  falling  on  the  nucleus.  The  latter  result  suggests  that  some  of  the  previously  held 
concepts,  claiming  the  Coulomb  force  being  balanced  by  the  total  quantum  potential  force 

Total  /rir ,  or  by  the  centrifugal  force  alone,  are  invalid. 

It  was  also  found  that  the  diffusion  kinetic  energy  dominates  over  the  dilatation  energy  for 
relatively  large  principal  and  angular  momentum  numbers,  except  in  the  ground  state.  As  a  result, 
an  atom  in  an  excited  state  is  portrayed  as  a  diffusive  probability  fluid  with  high  diffusion  kinetic 
energy  in  comparison  of  the  dilatation  energy. 

A  uniquely  interesting  fluid  dynamics  feature  of  a  hydrogen  atom  with  a  finite  magnetic 
quantum  number  is  the  presence  of  a  line  vortex  along  the  polar  axis.  The  vortex-induced  flow  is 
found  to  be  driven  by  the  sum  of  the  quantum  potential  due  to  the  polar  diffusion  fig  and 

centrifugal  potential,  fl, .  However,  only  a  part  of  polar  dilatation,  911^  ,  provides  the  kinetic 

energy  of  the  vortex  induced  flow.  The  remaining  energies  are  balanced  with  cross  inter-coupling 
modes,  and  thus  do  not  contribute  in  the  principal  modal  balance.  Specifically,  the  part  of  the  polar 
dilatation,  971^  canceled  with  the  centrifugal  barrier  potential  11, . 

By  combining  the  electrostatic  and  fluid  dynamic  features,  as  described  above,  a  hydrogen 
atom  with  a  finite  magnetic  quantum  number  is  envisioned  as  a  fluid  structure  having  a  line 
vortex-flow  with  a  dominant  core  that  extends  beyond  the  Bohr  radius  along  the  polar  axis.  The 
electrostatic  potential  becomes  dominant  in  the  region  exterior  of  the  vortex  core,  as  shown  in 
figure.  13.  For  an  atom  with  m  =  0,  the  vortex  structure  is  absent,  and  the  atom  has  a  Coulomb 
potential  dominated  core.  Our  model  approximated  the  nucleus  to  be  located  at  the  center-of-mass 
between  the  proton  and  electron.  The  latter  approximation,  however,  places  a  limit  on  the  validity 
of  model:  the  solution  is  invalid  within  the  spherical  region  of  the  radius  r  =  (me  /mp)re,  where  re  is 
the  electron  coordinate  measured  from  the  center-of-mass,  as  shown  in  figure  13. 


34 


Modal  balance  analysis  shows  a  number  of  unique  features  of  quantum  potential  in  a  hydrogen 
atom.  First  of  all,  Bohm’s  description  of  non-locality  of  the  quantum  potential  is  indeed  observed 
but  in  a  rather  discriminate  manner.  Of  all  the  elemental  processes  of  the  quantum  potential  in  a 
hydrogen  atom,  the  kinetic  energy  due  to  the  asymptotic  diffusion,  Qa=  ‘/z  Va  ,  is  constant 
throughout  the  full  spatial  domain  and  henceforth  falls  in  the  Bohm’s  characterization  of  non-local 
potential  that  does  not  fall-off  with  distance. 

It  is  interesting  to  mention  that  in  our  recent  studies,  we  also  arrived  at  a  conclusion  that  the 
Bohmian  non-local  dilatation  energy,  ,  for  a  harmonic  oscillator  with  a  spring  force  potential, 

which  is  proportional  to  the  square  of  a  position  coordinate,  is  equal  to  the  specific  zero-point 
energy,  whereas  the  sum  of  both  non-local  kinetic  energy  and  dilatation  energy  gives  the  n-th 
quantized  energy.  Secondly,  there  are  elemental  processes  of  quantum  potential,  which  do  not 
exhibit  non-local  behaviors.  For  example,  all  the  dilatation  energy  and  the  kinetic  energy  of 
diffusion,  with  the  exception  of  asymptotic  diffusion,  in  a  hydrogen  atom  do  not  exhibit  non-local 
behavior,  i.e.,  the  effects  of  the  quantum  potential  decay  with  distance.  These  latter  parts  of 
quantum  potential  contribute  in  the  local  dynamic  equilibria  associated  with  the  Coulomb  force, 
centrifugal  force,  and  vortex  induced  flow. 

It  is  relevant  to  mention  that  the  results  obtained  from  the  present  study  preserve  the  same  level 
of  knowledge  as  rendered  by  standard  quantum  mechanics.  This  is  so  because  the  conservation 
equations  of  quantum  fluid  dynamics,  which  are  used  in  the  present  modal  balance  equations,  are 
formulated  from  the  Schrodinger’s  equation  through  Hamilton-Jacobi  transformation. 
Additionally,  the  solution  of  Schrodinger’s  equation  for  a  hydrogen  atom  has  been  used  to 
construct  the  solution  of  the  present  modal  balance  problem.  The  merit  of  the  present  approach  and 
its  results  are,  however,  in  the  provision  of  the  physical  interpretation  of  quantum  mechanical 
energetic  and  dynamic  equilibria  of  an  atom  in  the  framework  of  the  quantum  theory  of  motion. 

It  is  remarkable  to  conclude  that  within  the  framework  of  a  quantum  fluid  description,  results 
from  a  spectrum  of  interesting  problems  recently  investigated,  suggest  that  the  quantum  energetic 
balance,  dynamic  equilibria,  and  the  organized  structures  of  the  atoms,  molecules  and 
many-particle  systems,  are  manifestly  governed  by  the  quantum  potentials  of  the  sets  of  elemental 
processes  associated  with  quantum  diffusion  provoked  by  non-uniform  probability  density 
distribution  in  quantum  systems  with  an  external  force  field  and  interaction  potentials.  The 
propositions,  described  above,  are  the  potential  areas  of  future  research  for  their  validation.  The 
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conclusion,  drawn  above,  appears  to  be  in  full  accordance  with  Holland’s  profound  view  that  “this 
new  (quantum)  theory  of  motion  introduces  a  new  conception  of  matter,  which  describing  that  the 
internal  potential  is  an  organizational  or  self-referential  form  of  energy  which  brings  about  an 
inner  tension  in  a  the  material  system  to  which  the  mass  points  respond”,  Holland  (1993). 
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Table  1 .  Dilatation,  Tl ,  and  diffusion  kinetic  energy,  Q ,  of  quantum  elemental  processes  in  hydrogen  atom 
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Table  2.  Modal  energy  balance  of  hydrogen  atom  at  selected  states 
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Figure  1.  Quantized  energy  is  the  diffusion  kinetie  energy  of  a  negatively 
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Figure  2.  Fractional  balance  of  Coulomb  and  quantum  potentials  at  ground  state  |l00) 

=  Q  =  —s  and  =  — 1/<^  . 


H.H.  Chiu 


(^) 

(0 


^=r/a 


Figure  3.  Complementary  modal  balance  of  Coulomb  and  quantum  potential  at  state 
1 200)  ,v,=-l/^ 
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Figure  4.  Fractional  modal  balance  of  Coulomb  and  quantum  potentials  at  state 
|210)  ^  and  v^=-\l^. 
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Figure  5.  Complementary  balanee  of  Coulomb  and  quantum  potentials  at  state  |300) 
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Figure  6.  Fraetional  modal  balance  of  the  Coulomb  potential  with  dilatation  and  diffusion 
kinetic  energy  for  the  selected  stated  :  1 1 OO)  ,|210),|21  +  l),|320),|32  +  2)  ,and 
CO  with  I  =  n-\ 
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Figure  7.  Variation  in  the  relative  importanee  of  diffusion  kinetic  energy  and  dilation  energy 
in  Coulomb  fractional  balance  mode.  In  complemental  mode  represents  of 

the  figure  also,  the  asymptotic  value  of  the  relative  magnitude  of  Qj,  and 
911  „  at  large  radial  distance. 
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Figure  8.  The  ratio  of  (^)  and  (^)  for  Coulomb  mode  balance  at  state  1 300) 
The  asymptotic  value  of  the  ratio  approaches  to  n-1  =2  at  >oo  . 


Figure  9.  Fraction  of  dilatation  energy  and  diffusion  kinetic  energy  contributing  in 
Coulomb  mode  balance  at  state  |300  >  .The  asymptotic  values  of 

and  approach  to  1  and  0  respectively. 
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Figure  10.  Distributions  of  diffusion  kinetic  energy  and  the  ratio  of  andQ^  at  1 300) 
The  asymptotic  value  of  approaches  to  1  at  ^co  . 
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Figure  11.  Electronic  and  hydrodynamic  structure  featured  by  vortex  dominated  inner  core  and 
Coulomb  potential  dominated  exterior  region  for  |433  >  . 
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Figure  12.  Dependence  of  the  vortex  dominated  core  size  on  the  principal  quantum, 

orbital  angular  momentum  quantum  number  and  magnetic  quantum  number 
for  |n,n-l,n-l )  state. 
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Figure  13.  Schematic  of  electronic  and  fluid  dynamic  structure  of  hydrogen 
atom  at  the  state  |32±2  >  ,  exhibiting  the  modal  balance  between  the  principal 
modes  ,  quantum  dilatation  and  diffusion  kinetic  energy  modes. 
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Abstract 


A  theory  of  dual  electron  beam  nanojets  is  developed  from  the  perspective  of  quantum  fluid 
dynamics.  The  structural  and  dynamic  behavior  of  the  dual  beams  are  investigated  over  a  large 
range  of  quantum  Reynolds  numbers,  N.  For  low  to  intermediate  N  <  100,  the  dual  beams 
exhibit  a  strong  wave-like  character.  This  wave-like  character  manifests  itself  as  a  merging  of 
the  two  separate  jet  streams  into  a  multiple  branched  beam  with  each  branch  fanning  out  at 
different  angular  displacements  relative  to  the  jet  axis,  overlapping  a  creating  and  interference 
pattern.  As  N  increases  the  merging  and  branching  occur  further  away  from  the  jet  nozzles. 
At  higher  N  >  400,  branching  occurs  at  a  considerable  distance  from  the  jet  nozzles;  the  fanning 
collapses  into  a  forward  propagating  beam  with  little  angular  spread.  Eventually,  at  very  high 
N  >  800  the  beam  splits  into  two  distinct  forward  propagating  beams  very  much  like  that 
produced  by  classical  jets.  Quantum  clustering  is  another  feature  that  arises  when  N  >  20. 
This  is  the  formation  of  regions  in  the  jet  with  abnormally  high  probability  density  and  is  the 
combined  result  of  self-interference  and  tunneling  within  each  branch.  The  clusters  shift  down 
stream  with  higher  N.  The  size  and  density  of  the  cluster  also  increases  with  an  increase  in  N. 
The  complexities  of  these  properties  are  studied  in  detail.  The  unique  properties  of  a  dual  beam 
quantum  nanojet  are  applicable  to  a  variety  of  nano-scale  applications,  including:  mass 
spectroscopy,  imaging,  sensors,  and  lithography. 


Key  Words:  quantum  nanojets,  quantum  Reynolds  number,  branching,  clustering 
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I.  INTRODUCTION 


Electron  double-slit  experiments,  resulting  in  the  interference  fringes  of  their  associated 
matter  waves,  have  long  been  at  the  heart  of  quantum  theory.  These  dual  beam  experiments 
have  provoked  a  historical  controversy  as  to  the  interpretation  of  their  exact  nature,  e.g., 
Copenhagen  interpretation,  pilot-wave  theory,  hidden  variable  theory,  and  causal  effects  which 
more  recently  have  been  folded  into  the  quantum  theory  of  motion'.  There  have  been  extensive 
and  impressive  experimental  studies  focused  on  the  interference  phenomenon  and  its  possible 
applications  (See,  for  example.  Refs.  2,  3,  and  4).  Nevertheless,  little  in  the  way  of  analytical 
studies,  aimed  at  a  detailed  explanation  or  clarification  of  the  dual  beam  interference 
characteristics  over  a  broad  range  of  operating  conditions,  has  been  reported.  Phillipidis,  Bohm, 
and  Kayes^  were  among  the  first  to  examine  analytically  a  two-slit  jet  operating  at  high  particle 
energy.  Their  work  suggests  that  the  quantum  theory  of  motion,  based  on  the  causal  theory  of 
de-Broglie-Bohm,  explains  how  a  “quantum  potential”  plays  a  dominant  role  on  the  electron 
behavior  in  the  dual  beam  interference  phenomenon.  According  to  the  causal  theory,  the 
quantum  potential  is  an  organizational  or  self-referential  form  of  energy,  which  brings  about  an 
inner  tension  in  the  material  system  to  which  the  mass  points  respond  (See  Ref  1).  The  present 
work  builds  upon  our  previous  efforts^’’  that  used  the  theory  of  de-Broglie-Bohm  as  a  starting 
point. 

The  present  work  contrasts  well  with  the  seminal  work  of  Moseler  and  Landman*  (ML),  and 
that  of  Eggers.^  ML  studied  the  formation,  stability,  and  break  up  of  a  liquid  jet  a  few 
nanometers  in  diameter  through  the  use  of  a  molecular  dynamics  simulation.  They  compared 
these  results  with  those  obtained  from  a  classical  continuum  hydrodynamics  model.  The 
atomistic  description  was  related  to  the  hydrodynamics  model  through  their  derivation  of  a 


stochastic  lubrication  equation,  which  allow  for  the  description  of  free  surface  flow  problems. 
Eggars  extended  that  work  by  employing  a  path  integral  formulation,  which  describes  the  most 
probable  configuration  for  jet  break  up. 

In  both  the  ML  and  Eggars’s  calculations  classical  Newtonian  laws  of  motion  govern  the 
particle  motion  and  no  quantum  effects  are  present.  The  present  work,  however,  approaches  the 
problem  from  the  viewpoint  of  quantum  fluid  dynamics.  This  allows  us  to  study  the  wave-like 
features  of  the  jet  particles  with  low  to  intermediate  quantum  Reynolds  numbers.  At  higher 
quantum  Reynolds  numbers,  viz.,  higher  energies,  our  theory  automatically  approaches  the 
classical  regime  of  the  particle  motion.  That  is,  since  we  have  limited  our  theory  to  low  particle 
densities,  the  theory  at  high  particle  energies  produces  results  that  would  be  obtained  from  the 
use  of  the  conservation  equations  of  classical  fluid  dynamics. 

In  summary,  the  work  of  ML  and  Eggars  focuses  on  the  classical  regime  of  the  jet  particle 
motion  whereas  the  present  work  focuses  on  the  quantum  regime;  the  two  approaches 
complement  each  other. 

In  section  II  we  use  the  independent  single  particle  model  to  introduce  definitions  and 
concepts  necessary  for  the  transition  to  a  quantum  fluids  model.  Based  on  a  quantum 
mechanical  formalism,  specific  elements  of  the  mathematical  analysis,  which  relate  to  quantum 
branching  and  clustering  are  detailed  in  section  III.  In  section  IV  we  discuss  the  structural  and 
dynamical  aspects  of  the  results  of  our  dual  nanojet  beam  calculation.  Lastly,  in  section  V,  we 
summarize  the  salient  features  of  our  results,  discuss  modifications  of  the  theory  for  our  next 
publication,  and  review  possible  nanojet  technologies  for  practical  use. 


II.  DEFINITIONS  AND  CONCEPTS 


From  elementary  quantum  theory  the  total  wave  function  of  a  dual  slit  is  given  as 
V  =  V  siiti  +  W  siit2  •  Each  of  the  slit  wave  functions  satisfies  the  Schrodinger  equation  for  an 
independent,  single  particle.  Boundary  conditions  are  built  into  each  of  the  slit  wave  functions. 
Squaring  the  above  expression  results  in  a  dual  jet  beam  probability  density,  which  includes  the 
necessary  interaction  terms.  With  this,  a  parametric  study  can  now  easily  be  carried  out  by 
varying  the  slit  width,  the  slit  separation  distance,  and  the  energy  of  the  emerging  particle  plotted 
as  a  function  of  the  axial  distance  from  the  slit  screen.  This  is  a  very  direct  approach. 
However,  the  approach  we  take  in  the  present  paper  is  to  view  the  system  as  a  quantum  fluid. 
In  this  model  we  do  not  separate  out  the  wave  functions  for  each  individual  slit  but  construct  a 
single  wave  function,  representing  the  dual  slit  system,  directly  from  the  slit  boundary  conditions. 
While  slightly  more  complicated  this  approach  results  in  a  physical  interpretation  of  quantum 
interference  and  tunneling  which  leads  to  the  processes  of  quantum  branching  and  quantum 
clustering.  A  schematic  of  the  dual  slit  system  we  will  treat  in  this  paper  is  shown  in  Fig.  1 . 

We  start  by  reviewing  the  properties  of  an  independent,  single  particle  solution  to  the 
Schrodinger  equation,  from  which  the  solutions  of  a  multiple  particle  dual  jet  stream  are 
constructed.  Since  the  electrons  are  taken  as  free  particles  of  energy  E,  we  assume  their  wave 
functions  are  stationary  state  plane  waves  and  are  of  the  form  'F(x,y,t)  =  \|/(x,y)  exp(-iEt  /  h) . 
The  time-dependent  Schrodinger  equation  of  motion  governing  the  electrons  is, 


ih 


d'F 


(1) 


where  h  is  the  Planck’s  constant  divided  by  2k,  and  nig  is  the  rest  mass  of  an  electron.  The 
probability  density,  p,  of  the  electron  is  then, 
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P=W  , 


(2.a) 


* 

where  V|/  is  the  complex  conjugate  of  the  wave  function  V|/  .  A  probability  current  density,  J,  is 
defined  as, 


J  = 


=  pa, 


(2.b) 


(2.C) 


where  the  particle  velocity  u  is  defined  as  the  ratio  of  the  probability  current  flux  J  with  the 
density.  These  two  quantities  are  related  through  a  continuity  equation,  which  is  analogous  to 
what  one  finds  in  classical  fluid  dynamics.  The  transition  from  a  wave  mechanical 
interpretation  to  a  fluid  mechanical  one  begins  with  a  focus  on  these  two  quantities. 

A  time  independent  equation  is  obtained  by  the  substitution  of  the  stationary  state  T(x,y,t) 
into  Equation  (1).  Casting  the  time  independent  Schrodinger  equation  into  a  non-dimensional 
form  allows  us  to  define  a  quantum  Reynolds  number,  N:  the  first  step  toward  our  fluids  model. 
Thus,  with  the  appropriate  algebraic  manipulation  and  redefining  of  independent  variables  such 
that  X  — >  Tix  /  b  and  y  ^  iry  /  b  ,  where  2b  is  the  slit  separation  in  the  dual  slit  problem  shown 
in  Fig.  1,  the  time  independent  non-dimensional  Schrodinger  equation  becomes, 

vVx,y)+NVx,y)=0,  (3) 


where 


N"  = 


2m  E  b^ 


(4.a) 
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which  may  also  be  written  in  terms  of  the  deBroglie  wavelength,  /  2m^E , 


N' 


(4.b) 


where  =  b  /  7t . 

Alternatively,  N  can  be  viewed  as  a  quantum  Reynolds  number  that  is  defined  by  the  ratio  of 
typical  inertial  force  to  quantum  viscous  force.  Thus  N  =  u^b  /  (h;:  /  m^)  =  u^b  /  Vq  ,  where  me 


is  the  rest  mass  of  electron,  u^  the  initial  velocity  of  a  particle  corresponding  to  its  energy  E, 
and  Vq  =  h7r  /  m^  is  the  quantum  viscosity. 

The  dual  slit  wave  function  in  the  fluids  model  is  written  as  the  sum  of  several  Fourier  transform 
integrals.  Similar  to  classical  optics,  in  which  there  are  near  and  far  regions  of  diffraction,  these 
integrals  represent  physical  distinct  regions  as  well.  The  definition  of  these  regions  can  be 
obtained  by  taking  the  Fourier  transform  of  the  Schrodinger  equation  given  in  Equation  (3).  The 
resulting  equation  is, 


^  +  (N^-7=)(|>  =  0, 


(5) 


where  wave  function  ^  is  given  by  the  Fourier  transform, 


(|)(X,77) 


e'''^  dy . 


(6) 


In  transforming  Equation  (3),  the  boundary  conditions  for  v|/(x,y)  and  its  spatial  derivative  with 
respect  to  y  are, 

v(x,y)  ly^iarge^  0  D  D  D  y)/5y ^  0 ,  (7.a,  b) 
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for  all  values  of  x.  In  Appendix  A  we  derive  a  detailed  expression  for  the  total  wave  function  in 
terms  of  a  sum  of  Fourier  integrals.  In  Appendices  A  and  B  we  apply  boundary  conditions  to 

the  initial  particle  velocity,  u^ ,  and  initial  number  of  particles,  ,  emanating  from  the  slit,  as 

well  as  to  the  wave  function  and  its  spatial  derivatives,  which  results  in  an  expression  for  the 
wave  function  of  a  multiple  particle  dual  jet  beam.  This  wave  function  is  succinctly  written  as: 

v|/(x,y)  =  v|/^(x,y)  +  vi/g(x,y)  +  v|/^(x,y)  .  This  component  notation  will  be  used  in  the 

following  sections  in  which  the  sub  and  superscript  notation  are  defined. 

Equation  (5)  possesses  two  regions  of  validity  for  its  solution.  The  first  is  in  the  region 
where  1 77 1  <  N .  We  call  this  region  the  sub-critical  region.  In  this  region  the  general  solution 

is  made  up  of  two  propagating  waves  of  the  form  exp(±iax) ,  where  a  =  +-^N^  -  r\  .  Similarly 
the  second  region  is  where  |  7 1  >  N ,  is  call  the  super-critical  region.  The  general  solution  in 

this  region  has  the  form  exp(-  p  x ) ,  where  p  =  .  This  has  some  elements  of  analogy 

to  particles  incident  on  a  step  potential  in  one  dimension.  The  two  propagating  waves  in  the 
sub-critical  region  give  rise  to  interference,  whereas  the  solution  in  the  super-critical  region 
decays  exponentially  is  viewed  as  tunneling.  This  analogy  will  be  exploited  later  sections  of 
the  paper. 

For  computational  purposes  to  be  used  later  in  the  paper,  we  define  the  transformed  functions 
F(rj)  and  G(  rj) : 


F(r])  =  <^(0,rj)  = 


G(?j)  = 


dx 


1 


\|/(0,y)  e'''My, 

—00 

-#(0,y) 

-«>  dx 


(8.a) 


(8.b) 
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Equation  (8)  corresponds  to  the  boundary  conditions  on  v|/(x, y)  and  its  derivative  with  respect 


to  X  denoted  by  V|/(x,y)  lx=o,y=siit=  [5¥(x,y)/dx]^^^  =  5v(O,y)/0x  .  That  is, 

they  are  evaluated  at  x  =  0  and  over  width  of  the  slit. 

We  have  now  laid  the  groundwork  for  the  dual  slit  problem  viewed  as  non-interacting  multi¬ 
particle  system,  which  will  be  developed  in  the  next  section.  The  physical  parameters  of  that 
system  are  as  follows.  We  consider  dilute  jets  of  non-interacting  electrons  emitted  from  a  dual¬ 
slit/nozzle  injector  in  the  x-direction  with  an  energy  E  into  a  vacuum  of  two-dimensional  space, 
in  the  absence  of  an  external  potential.  The  slits  are  separated  by  a  distance  of  2b,  and  their 
widths  are  taken  to  be  0.4b.  Again,  a  schematic  diagram  of  this  dual-slit  configuration  is  shown 
in  Fig.  1. 

III.  WAVE  PROPAGATION,  INTERFERENCE  AND  TUNNELING 
A.  Sub-critical  wave  field:  |  ;;  |  <  N 

Since  1 77 1  <  N  it  is  advantageous  to  write  77  =  N  -  4"  so  that  the  sub-critical  wave  field  can  be 
expressed  as  the  sum  of  four  integrals  with  respect  to  4"  as  shown  in  Equation  (9).  The 
separation  of  the  wave  function  as  such  is  explained  in  Appendix  A.  The  subscript  notation  of 
F  and  G  on  the  wave  functions,  \\i  and  'F  ,  is  related  to  the  Fourier  transforms  of  the  boundary 
conditions  given  in  Equation  (8).  The  superscript,  <,  is  just  a  reminder  that  these  functions  are 
associated  with  the  sub-critical  region  1 77 1  <  N .  Thus, 
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„  G,,  ,  ,Jsin(N-f)-™,[0.6(N-f)]l 

/I  ^  ^  ^ 

471  [  V2N4'-4'"(N-4') 

(lO.b) 

(lO.c) 

III 

II 

1 

(lO.d) 

(lO.e) 

Corresponding  expressions  related  to  the  coefficients  and  are  derived  in  Appendix  A.  The 

coefficients  and  are  scaling  factors’,  which  are  derived  in  Appendix  B. 

Their  functional 

forms  are: 

II 

'  hJ  ' 

o 

(11. a) 

m  mb 

hji 

(ll.b) 
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Equation  (9)  may  be  rewritten  in  a  more  succinct  form  in  terms  of  what  we  have  designated 
as  ^  waves,  each  of  which  has  a  wave  characteristic  function  defined  in  the  exponential  terms 
in  Equation  (9), 

'P;(x,y,t,N)  +  4';l(x,y,t,N)  =  ^;+^+<;+4;.  (12) 

Each  term  corresponds  respectively  with  those  in  Equation  (9).  Here  again,  wave  propagation 
is  in  both  the  “forward”  and  “backward”  directions,  which  coincide  with  the  directions  normal  to 
the  corresponding  wave  characteristic  surfaces,  thus  the  subscripts  f  and  b  on  the  4"  waves. 
Moreover,  these  waves  propagate  at  a  radial  angle,  0,  with  respect  to  the  x-axis.  In  fact,  they 
propagate  in  a  symmetric  fashion  above  and  below  the  x-axis,  thus  the  superscripts  plus  (+)  and 
minus  (-),  respectively. 

The  angle  of  propagation  is  given  by, 

0  =  +tan“^  =  +tan“' 

K  . 

where  the  plus  (+)  sign  applies  to  the  Ci  C  waves,  and  minus  (-.)sign  to  Q  and 

^  waves,  respectively.  A  review  of  Equation  (13)  reveals  the  following. 

For  small  values  of  N,  and  fixed  i^,  0  is  large.  On  the  other  hand,  when  N  takes  on  larger 
values  0  decreases  as 


lim0:±|— I  ^0,  (14) 

Essentially,  the  beam  splits  into  two  distinct  forward  propagating  beams  parallel  to  the  x-axis 
with  a  structure  very  much  like  that  produced  by  classical  jets.  This  suggests  that  at  some  point. 
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a  transition  occurs  between  quantum  behavior  and  that  of  a  classical  sand/grain  like  jet.  Our 
study  shows  that  this  transition  occurs  at  a  critical  value  for  the  Reynolds  number  of  Nc  =  N  ~ 
100.  As  we  will  see  below,  low  Reynolds  numbers  are  associated  with  a  predominately  wave¬ 
like  behavior  consistent  with  quantum  behavior  and  high  Reynolds  numbers  give  rise  to  a  purely 
particle  picture  consistent  with  a  classical  behavior. 

Another  feature  of  note  is  that  for  a  fixed  value  of  N  the  Fourier  components  with  a  larger 
values  of  ^  (equivalently  smaller  values  of  has  a  large  radial  angle.  This  basic  rule  helps  us 

to  understand  the  nature  of  wave  structure  of  the  nanojets  as  we  shall  see  in  the  section  4. 

B.  Super-critical  wave  field:  | ;/ 1  >  N 

As  with  the  sub-critical  region  the  super-critical  wave  field  can  be  expressed  as  a  Fourier 
transform  with  respect  to  4",  which  is  also  shown  in  Appendix  A. 


'F;(x,y,t,N)  =  \i/Ux,y)e“ 


(15) 


j*  exp(-px)  cos  (N  -I-  exp(-iEt  /  h)  d^ 


The  subscript  F  on  the  psi  functions  is  again  related  to  the  Fourier  transforms  of  the  boundary 
conditions  given  in  Equation  (8).  The  superscript,  >,  is  just  a  reminder  that  these  functions  are 

associated  with  the  super-critical  region  |7|>N.  The  parameter  P  =  (2N(^-l-(^^)^^^.  As  with 

and  B^,  the  coefficient  has  a  related  expression  derived  in  Appendix  A.  Then,  C^is 

identified  as, 


C  ^0.^  ,  Jsin(N  +  40-sin[o.6(N  +  40] 
^  71  (N  +  0 


(16) 
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So  in  this  super-critical  region  we  have  a  standing  wave,  which  decays  exponentially  in  an 
axial  direction  as  indicated  by  the  exp(-px)cos(N  +  Qy  exp(-iEt  /  h)  factors.  The  axial  decay 

is  caused  by  quantum  tunneling.  We  note  that  when  x<(2Ni^+4^)  i.e.,  there  is  negligible 

amplitude  decay  resulting  in  a  high  probability  density.  However,  when  x>(2N(^+(f)  '^^the 
wave  decreases  rapidly.  Later  in  the  paper  we  will  show  that  the  formation  of  a  cluster  down 

stream  is  in  part  due  to  a  region  in  which  x<(2Ni^+4"^)  is  followed  by  a  rapid  exponential 
decay. 

As  with  the  sub-critical  region  it  is  advantageous  for  further  discussion  in  later  sections  to 
divide  this  region  into  its  own  set  of  4"  waves  as  follows. 

ft  00 

'T^(x,y,t,N)  =  C^exp(-px)  exp(+i2y^y)exp(-iEt  /  h)d4' 

ft  00 

-h  C^exp(-px)exp(-i/l  y)exp(-iEt  /  h)  (17) 

=  (l+Q 

where  we  have  defined  2^^  =  N  -i-  Here  again,  we  have  identified  two  terms  as  ^waves, 

Q  and  Q ,  with  plus  (+)  and  minus  (-)  superscripts  referring  to  their  propagation  in  the  direction 

of  increasing  and  decreasing  values  of  the  y-coordinate,  respectively,  and  a  subscript  S  that 
indicates  these  waves  are  standing  waves. 

IV.  STRUCTURAL  AND  DYNAMIC  COMPLEXITIES  OF  ELECTRON  DUAL-SLIT 
NANOJETS 

In  the  present  work  our  numerical  studies  focus  on  the  variations  in  the  structure  and 
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dynamics  of  a  two-slit  electron  jet  having  quantum  Reynolds  numbers  ranging  from  10  to  1000. 
This  corresponds  to  the  initial  velocities  (kinetic  energies)  ranging  from  3.6x10  cm/s  to  3.6x10 
cm/s  (0.004  eV  to  37  eV),  respectively.  The  geometric  parameters  used  in  the  calculations  were: 
slit  separation  of  2b  =  200  nm,  and  a  slit  width  of  0.4b  =  40  nm.  The  initial  probability  density 

at  the  slit  inlet  was  normalized,  iio=  1.  Finally,  given  the  conditions  above,  the  current  density 

per  unit  width  of  slit  in  y-direction  ranges  from  9.0  x  lO"  particle  sec''  to  9.0  x  lO'^  particle  sec  ’, 
respectively. 

The  analysis  of  our  numerical  results  is  broken  up  into  to  primary  regions  of  quantum 
Reynolds  numbers:  10  <  N  <  100,  and  100  <  N  <  1000.  When  N  is  approximately  100  the 
systems  exhibits  a  transition  from  a  quasi-quantum  regime,  just  below  N  =  100,  to  a  quasi- 
classical  region,  just  above  N  =  100.  This  critical  value  of  N  is  designated  as  Nc  =  100. 

We  further  divide  the  analysis  into  two  ranges  of  the  coordinate  x.  The  range  of  0  <  x  <  5  is 
a  near  field  region,  or,  in  analogy  with  classical  optics,  the  Fresnel  region.  The  range  of  5  <  x  < 
30  is  the  far  field  region,  or  the  Franuhofer  region.  The  Fresnel  and  Fraunhofer  regions  are 
illustrated  by  our  results  depicted  in  Figs.  2a  -  2f,  and  Figs.  3a  -  3f,  respectively. 

A.  Quantum  branching  and  structural  transition:  quantum  to  classical 

In  general,  quantum  jets  corresponding  to  lower  Reynolds  numbers,  10  <  N  <  100,  exhibit 
“multi-branching”  effects.  Branching  is  clearly  observed  in  both  the  Fresnel  and  Fraunhofer 
regions  depicted  in  Figs.  2a  -  2c,  and  Figs.  3a  -  3c.  In  these  figures  the  magnitude  of  the 
electron  probability  density  parallels  the  visible  spectrum  with  the  lower  densities  on  the  blue 
end  of  the  spectrum  and  the  higher  densities  toward  the  red  end.  Note  that  each  branch  has 
associated  with  it  a  width  and  branching  angle  0  given  by  Equation  (13).  The  quantum 
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branches  are  formed  by  two  distinct  mechanisms. 


The  first  mechanism  is  the  generation  of  high  harmonic  ^  waves  downstream  from  the  slit 
where  the  electron  jet  splits  into  two  branches,  clearly  observed  in  Figs.  2a,  2b,  3a,  and  3b.  The 
second  mechanism  contributing  to  the  jet  branching  is  the  intersection  of  in  phase  and 
waves  on  the  x-axis.  Figures.  2a  and  2b  illustrate  the  formation  of  a  central  branch  parallel  to 
the  x-axis,  by  the  intersection  of  a  wave  emanating  from  the  upper  slit  with  a  wave 
emanating  from  the  lower  slit.  Moreover,  there  is  a  reflection  of  and  waves  back  into 
the  regions  from  which  they  were  generated,  where  upon  interference  with  the  primary  jets, 
which  emanate  from  the  upper  and  lower  slits,  takes  place.  The  strength  of  the  interference  is  a 
function  of  the  phase  between  the  reflected  wave  and  the  primary  one.  The  effect  of 
interference  can  clearly  be  observed  in  Figs.  2a  -  2c  for  the  lower  Reynolds  numbers  where  the 
branching  angles  are  the  largest. 

Figures  2a  -  2d  and  3a  -  3d  correspond  to  Reynolds  numbers  in  the  range  of  10  <  N  <  100. 
Following  the  evolution  of  the  jet  structure  from  N  =  10  to  N  =  100  in  these  figures  clearly 
shows:  1.  a  reduction  in  the  branching  angles  resulting  in  a  forward  propagation  of  the  jets  from 
each  emanating  from  each  slit,  2.  a  reduction  in  the  widths  of  the  branches,  and,  3.  the  formation 
of  the  central  branch  taking  place  further  downstream  from  the  jet  slits.  At  approximately  N  = 
100,  the  two  jets  are  for  the  most  part  disengaged  and  function  independently  of  each  other.  For 
this  reason  we  have  designated  N  =  100  as  Nc  the  critical  Reynolds  number.  For  N  <  Nc  the 
figures  demonstrate  basic  quantum  behavior  including:  spatial  density  fall-off,  complex 
interference  fringes  and  tunneling. 

For  values  of  N  >  Nc  the  system  continues  to  evolve  with  increasing  N  into  a  more  classical- 
like  structure.  Clearly,  in  the  Fresnel  region  when  N  «  ,  Fig.  2d,  the  branching  angles  have 
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pretty  much  collapsed  so  that  the  respective  branches  are  coincident  resulting  in  two  primary, 
forward  propagating,  sand/grain-like  jets  emanating  from  the  upper  and  lower  jet  slits. 
However,  this  is  not  the  case  in  the  Fraunhofer  region  as  can  be  seen  in  Fig.  3d;  multiple 
branching  continues  to  be  present  for  values  of  x  >  10.  As  N  is  increased  beyond  Nc  the 
formation  of  multiple  branching  occurs  further  and  further  downstream  from  the  jet  slits 
extending  the  range  of  the  primary  but  separate  jet  streams  trailing  the  branching  point  (See  Figs. 
3e  -  3f).  When  N»400  the  complete  range  of  0  <  X  <  30  covering  both  Fresnel  and 
Fraunhofer  regions  is  devoid  of  any  significant  branching  and  branch  interactions.  Essentially 
what  are  left  are  the  two  primary  forward  propagating  jets  of  electrons.  Fig.  2e,  N  =  500, 
clearly  demonstrates  this  effect  over  the  range  of  0  <  x  <  50  for  N  =  500. 

Finally,  when  N  »  Nc,  Figs.  2e,  2f,  and  3e,  3f,  the  two  primary  jets  become  quite  focused 
into  pencil  sharp  beams  of  particles  over  the  entire  range  of  x  propagating  in  the  forward 
direction  with  no  interaction  between  them,  which  is  characteristic  of  classical  particle  beams. 
Clearly  the  evolution  of  the  dual  slit  jet  structure  with  increasing  Reynolds  numbers  exhibits  the 
classic  wave-particle  duality.  At  the  lower  Reynolds  numbers  there  is  a  clear  indication  of 
wave  like  properties  with  branching  and  interference.  At  the  higher  Reynolds  numbers  we  have 
observed  a  significant  reduction  in  branching  resulting  in  little  of  no  interference  and  two 
forward  propagating  beams  of  classical  particles  propelled  through  their  respective  slits. 

In  the  following  two  sections,  B  and  C,  we  take  a  closer  look  at  the  phenomena  of 
quantum  clustering  and  interference. 

B.  Quantum  clustering  and  axial  structure 

In  Fig.  4  we  have  plotted  the  electron  probability  density  as  a  function  of  axial  distance  for 
various  values  of  quantum  Reynolds  numbers  over  the  entire  range  of  10  <  N  <  1000.  The 


-14- 


apparent  broadband  outlining  of  these  plots  is  in  fact  due  to  the  rapid  oscillations  of  density 
within  a  short  distance  of  the  wave  propagation.  Most  noticeable  is  the  growth  of  a  large 
maximum  whose  initial  presence  can  be  seen  in  Fig.  4a,  where  N  =  10,  in  the  very  near  field 
regime,  x«0.1.  For  values  of  N  >  20,  this  peak  quickly  grows  to  a  value  of  1.8,  which  is  1.8 
times  greater  than  the  initial  density  at  the  slit  surface  show  in  Fig.  4a.  The  formation  of  an 
abnormally  high-density  region  suggests  a  “quantum  clustering”  of  the  probability  density, 
initially  introduced  and  described  in  section  3.  As  the  Reynolds  number  increases  the  cluster 
shifts  downstream  away  from  the  origin  at  the  slit;  in  addition,  it  broadens  out  covering  a  larger 
range  of  the  propagation  distance.  This  progression  can  be  clearly  tracked  by  sequential 

viewing  figures  4a  -  4e.  Quantum  clustering  can  be  attributed  to  self-intersection  of  the  and 
^  waves  in  each  jet  stream  (see  Equations  (12)  and  (17)).  Since  the  wave  angle  Oof  the 

and  waves  becomes  smaller  at  large  value  of  N,  the  clustering  occurs  at  larger  x  for  higher 

value  of  N.  The  interference  of  these  two  waves  can  clearly  be  seen  in  the  density  plots.  Fig.  3c, 
3d,  and  3f,  where  they  intersect  each  other  on  the  axis  of  each  jet  stream.  The  standing  waves 

due  to  also  contribute  to  the  higher  density  in  the  region  of  smaller  x,  and  reduced  density  at 
larger  x  due  to  axial  tunneling. 

Finally,  in  Fig.  5  we  plot  a  series  of  probability  density  curves  for  different  values  of  N.  By 
placing  different  but  increasing  N  value  density  curves  next  to  each  other  we  are  able  to  display 
the  evolutional  progression  of  cluster  growth  in  one  figure.  In  this  figure  La,  Lb,  and  Lc  are  the 
displacement  values,  on  the  jet  axis  of  propagation,  where  the  cluster  initially  starts  to  grow,  for 
the  cluster  maximum,  and  for  the  effective  end  of  the  cluster,  respectively.  In  passing  we  note 
that  it  is  this  correlation  between  the  peak,  or  maximum,  values  of  each  N  valued  cluster  that 
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forms  the  foundation  for  using  nanojets  as  a  mass  spectrometer.^ 

C.  Quantum  interference  fringes 

Figure  6  is  a  parametric  study  of  the  effects  of  different  slit  separation  width  of  2b  on  the 
probability  density  and  its  corresponding  interference  patent,  where  b  takes  on  the  values  of  b  = 
80nm,  lOOnm,  150nm,  200nm,  250nm,  and  300nm,  corresponding  respectively  to  the  figure 
panels  (a  -  f).  It  shows  the  cross  sectional  contours  probability  densities  and  their 
corresponding  interference  patents  of  the  advancing  jet  streams  at  the  axial  displacement  of  x  = 
10  for  different  values  of  slit  separations,  and  for  a  value  of  N  =  10.  Calculations  have  been 
carried  out  for  several  values  of  x  spanning  the  near  (Fresnel)  and  far  (Franuhofer)  regions.  We 
have  chosen  to  report  the  results  at  x  =  10,  which  is  an  overlap  of  the  two  regions  and  typifies  the 
comparison  of  the  probability  density  and  its  corresponding  interference  pattern  in  areas  of 
meaningful  interest. 

The  effects  of  a  change  in  the  Reynolds  number  on  interference  fringes  are  displayed  in  Fig. 
7  as  a  function  of  axial  displacement.  We  note  that  for  the  lower  values  of  N,  approximately  N 
<  100,  the  familiar  multiple  fringe  patterns  due  to  interference  are  present.  For  values  of  N  > 
500,  the  interference  fringes  increasingly  fade  with  increasing  N  until  there  are  only  two  vertical 
lines  indicating  the  almost  complete  eradication  of  quantum  interference  with  two  distinct 
streams  of  classical  particles  emanating  from  the  slits.  Note  that  the  central  branch  is  present 
for  all  values  of  x  when  N  <  50.  It  does,  however,  disappear  when  N  =  100  at  x  =  5,  but 
reappears  at  x  =  15  and  30.  This  is  consistent  with  the  discussion  given  in  Section  IV  A  of  Figs. 
2  and  3  in  which  the  formation  of  the  central  branch  progressively  shifts  down  stream  away  the 
slits  with  increasing  value  of  N.  Eventually,  as  illustrated  in  Fig.  3i,  when  N  >  500,  the 
formation  of  the  central  branch  has  essentially  shifted  to  infinity  and  thus  only  two  branches,  one 
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for  each  slit,  are  displayed. 


V.  CONCLUSION 

The  theory  of  dual-slit  dilute  electron  nanojets  was  formulated  in  terms  of  a  quantum 
Reynolds  number,  N.  A  critical  Reynolds  number,  Nc  =  100,  was  identified.  For  values  of  N 
<  Nc  the  dual-slit  system  exhibited  the  quantum  features  of  branching  and  clustering.  And  for 
values  of  N  >  Nc  the  system  tended  toward  a  more  classical  picture  in  which  two,  apparently 
independent,  beams  of  electrons  propagate  in  straight  line  trajectories  from  each  slit. 

The  dual-slit  wave  function  was  constructed  from  the  solutions  of  the  non-interacting,  multi¬ 
particle  Schrodinger  equation  and  the  appropriate  boundary  conditions  on  the  jet  slits.  The  field 
of  propagation  was  resolved  into  two  regions:  the  sub-critical  region  of  propagating  waves  and 
the  super-critical  region  in  which  the  waves  decay  exponentially.  A  Fourier  decomposition  of 
the  wave  function  in  these  regions  showed  that  the  wave  function  could  be  written  in  terms  of  six 

basis  functions:  +  47  in  the  sub-critical  region,  and  ^-i-^in  the  super-critical 

region. 

Quantum  branching,  the  appearance  of  multiple  branches  radiating  at  different  angles  in  a  fan 
shaped  design,  is  induced  by  the  interaction  of  the  two  families  of  and  waves  emanation 
from  both  slits.  Quantum  clustering,  spatial  regions  of  high  probability  densities,  on  the  other 
hand,  arise  though  the  interaction  of  the  and  waves  emanating  from  the  same  slit. 
Moreover,  quantum  and  classical  behavior  are  manifest  in  the  wave  angle,  0.  When  the 
wave  angle  is  large  the  jets  exhibit  strong  quantum  structure,  whereas  a  small  angle  leads  to 
sand/grain  like  jet  structure  indicative  of  a  classical  jet  stream  of  particles.  Thus  the  waves 
are  key  to  the  interpretation  of  wave-particle  phenomena. 
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In  contrast  to  Copenhagen  interpretation,  the  causal  interpretation  by  de-Broglie  and  Bohm 
also  gives  the  nanojet  structure  based  on  the  quantum  fluid  dynamics  model.  The  two  rival 
interpretations,  however,  give  the  same  results.^ 

As  a  follow  on  to  this  paper  we  are  in  the  process  of  recasting  the  present  work  and 

O 

interpreting  its  results  in  terms  of  a  full  quantum  “diffusive”  fluid  dynamics  theory  .  Recasting 
the  theory  in  these  terms  will  allow  us  to  identify  a  quantum  potential.  With  this  potential  a 
more  fundamental  understanding  of  the  driving  mechanisms  of  the  problem  can  be  had.  The 
present  calculation  is  limited  to  relatively  dilute  densities  of  non-interacting  particles  within  the 
jet  stream.  Recasting  the  problem  in  terms  of  a  quantum  diffusive  fluid  dynamical  theory  will 
set  the  background  needed  to  address  moderate  to  high  densities  of  jet  particles  and  their  mutual 
interactions.  Quantum  diffusive  fluid  dynamics  has  been  successfully  used  by  one  of  the 
authors,  H.H.  Chiu,  to  identify  the  dynamic  origin  of  quantized  energy  and  the  mechanisms  of 
detailed  dynamic  equilibrium  within  a  hydrogen  atom  induced  by  the  diffusion  of  the  probability 
fluid  of  the  electron. The  new  aspect  of  the  theory  is  the  consideration  of  quantum  diffusion 
hereto  not  formally  addressed  to  any  extent  in  the  literature. 

Once  the  full  theory,  addressing  a  broad  range  of  particle  densities  in  the  jet(s),  is  completed 
it  is  our  intention  to  extend  the  problem  to  multi-beam  components,  such  as  different  ions.  At 
this  stage  we  will  have  a  theoretical  tool  that  may  be  used  for  parametric  studies  of  design 
parameters  to  guide  and  optimize  the  experimental  development  of  specific  device  applications. 

Finally,  nanojets  have  a  myriad  of  important,  modem  day  applications  for  which  our  model 
can  be  used:  nano-  circuitry  and  electronics;  nano-  propulsion  /  thmster  systems;  fuel  injectors 
for  microscopic  engines;  atomic  optics  and  interferometry;  etching  and  lithography;  drag 
delivery  systems;  injecting  genes  into  cells;  mass  spectrometer;  and  quantum  sensors.  Our  goal 
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with  the  present  publication  and  research  in  progress  is  to  provide  a  theoretical  foundation  to:  1. 
understand  the  fundamental  physics  which  governs  nanojet  operation  under  various 
circumstances;  and  2.  act  as  tool  in  guiding  the  design  and  optimizing  the  operational  functions 
for  the  various  applications. 
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APPENDIX  A:  DUAL  BEAM  WAVE  FUNCTION  AND  BOUNDARY  CONDITIONS 

The  time  independent  wave  function  of  two-slit  jet,  v|/(x,y)  =  v|/^(x,y)-l-v|/^(x,y)-l-vi/^(x,y), 

given  in  integrands  of  Equations  (9)  and  (15)  of  the  text,  is  obtained  by  the  method  of  Fourier 
transform  as  follows.  The  superscripts  “<”  and  “>”  refer  to  the  sub  and  super  critical  regions. 

When  the  single  particle  Schrodinger  equation.  Equation  (3),  is  Fourier  transformed  the 
resulting  equation.  Equation  (7),  is 

^  +  (N'-;;-)(li  =  0  (A.1) 

dx 

As  discussed  at  the  end  of  Sec.  II  of  the  text,  the  solution  of  Equation  (A.l)  possesses  two 
regions  of  validity.  The  first  is  in  the  region  where  |  7|<N:  the  sub  critical  region.  In  this 
region  the  general  solution  is  made  up  of  two  propagating  waves  of  the  form  exp(±iax),  where 

a  =  -  ri^  ,  or 

(|)<(x,ri)  =  qe"“’‘+C2e-'“’‘  (A.5) 

The  second  region  is  where  |  7  |  >  N  :  the  super  critical  region.  The  general  solution  in  this 

region  has  the  form  exp(- p  x) ,  where  p  = ,  or 
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(|)^(x,ri)  =  C3e-P’‘  (A.6) 

We  note  that  =  F^(r/)  and  d(|)^(0,/7)/dx  =  G^(/7) ,  where  F^(rf)  and  G^(rf)  are 

the  Fourier  transforms  of  \t/^(0,y)  and  d\|/^(0,y)/dx  given  in  Equations  (8. a)  and  (8.b), 
respectively.  With  this  we  can  determine  the  coefficients  CpC2,andC3in  terms  of 

F^(rf)  and  G^(rf) .  Then  (|)  ^(x,  rj)  may  be  written  in  the  following  general  form: 


where 


^^(x,r])  =  H^(x,r/)F^(r])  +  K^(x,r/)G^(r])  for  |;7|<N,  (A.7) 

ir(x,jj)  =  =  cos(ax) ,  (A.8) 

K^(x,j])  =  ,  (A.9) 

2\a^  a 

(A.10) 

K^{x,rj)  =  0.  (A.  11) 


The  real  space  solution  \|/^(x,y)of  the  Schrodinger  equation,  Equation  (3),  is  obtained  by 
taking  the  appropriate  inverse  Fourier  transforms: 


H^(x,7j)F  dy  +  i — 

yjln 

In  addition,  because  this  is  a  two-dimensional  problem  and  we  assume  reflection  symmetry  in 
the  vertical  plane  about  the  line  of  propagation  of  the  particle  beam  Equation  (A.  12)  may  be 
rewritten  as. 


[”  is:^(x,7)G^(;7)e  ‘”dy  (A.12) 

J  —00 


■\j2n  ^ 
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\(/>(x,y)  = //>(x,;7)F>(7)  +  Z>(x,7)G^(;7))  cos(7y)  d;/  (A.13) 

'J2n 


The  total  real  space  wave  function  is  a  linear  combination  of  the  wave  functions  of  the  sub 
and  super  critical  regions.  Using  Equations  (A. 8)  -  (A.13)  in  the  appropriate  regions  results  in 
the  total  wave  function: 

V(x,y)  =  H'p(x,y)  +  «(x,y)  +  = 


2  fN 


jE’^(N-i0cos  ^2N4'-4’^x  +  ,  =sin  lcos[(N -i0y]c/(^(A.14) 


e-V2Nf.rxcos 


[(N  +  Oy]d^ 


where  77  =  N-^  when|  ?7|<  N  ,  and  ?7  =  N  +  i^  when  |;7|>N. 

Equation  (A.  14)  can  be  recast  into  a  form  consistent  with  the  boundary  conditions  over  the 
dual  slits,  which  are  defined  in  Equations  (B2)  -  (B5).  The  final  result  of  this  dual  slit  wave 


function  is: 


1^  (x,  y)  =  1^  P  (x,  y)  +  1^  <  (x,  y)  +  1^  >  (x,  y) 


^  F  0<  (X,  y)  +  5  ©<  (X,  y)  +  ^  F  0>  (x,  y) 


where  ^  =  1  +  i,  d--\  +  i, 


(A.  15) 


Q^(x,y)  =  -^£^Y^(N-Ocos  +72N^-^'x  cos[(N-Oy]d<^,  (A.16.a) 


Q>(x,y)  = 


Y^(N  +  i0  exp  — y/2N^+^^x  cos[(N  +  i0y]di^,  (A.16.b) 


rTfsin(NmO-sin[0.6(NmO] 

Y>(N  mO  =  .  —  \ - — ^ ^ 

\  271  N 


(A.17.a,  b) 
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0^(x,y)  =  -^|^^Z^(N-Osin  cos[(N - Oy] d^,  (A.  18) 

V271  *-  -* 


sin(N  -  O  -  sin  [0.6  (N  -  O] 

V2N^^(N-0 


(A.  19) 


and  -fp  and  are  coefficients,  determined  by  a  scaling  method  described  in  Appendix  B,  that 
depend  on  the  density  and  velocity  of  the  particles  on  the  slit  surfaces. 

The  three  functions  Y^,Y^,  and  are  determined  consistent  with  the  boundary  conditions 
on  the  probability  density,  J,  and  the  particle  velocity,  u,  of  the  jet  beam  given  by  Equations 
(B.2),  (B.3),  and  (B.4)  in  Appendix  B.  The  probability  density  and  the  particle  velocity  of  the 
jet  are  calculated  by  a  substitution  of  the  dual  slit  wave  function,  Equation  (A.  15),  into  Equation 
(2).  The  methods  for  the  determination  of  Y  and  involve  Fourier  transforms,  whereas 

Z^(N— 4)  involves  the  solution  of  Fredholm  type  integral  equation.  In  the  determination  of 

the  solution  ofZ^(N— 4),  we  used  an  approximation  method.  In  all  cases  the  fractional 
difference — the  difference  between  the  calculated  mean  value  of  the  particle  velocity  and  its 
initial  boundary  value  velocity  divided  by  the  initial  boundary  value  velocity — was  on  the  order 
of  or  less  than  10'  .  For  very  large  N,  the  error  decreases  to  smaller  values.  For  example,  the 
latter  ratio  is  3.8  x  10'^  for  N  =  10,  and  2.6  x  10'"^  at  N  =  100,  and  is  1.5  x  10'"^  for  N  =  1000. 
The  integral  expressions  given  in  Equation  (9)  and  (15)  for  the  time  dependent  solutions, 

^^(x,y,t,N)  +  ^g(x,y,t,N)  and  ^^(x,y,t,N)  ,  respectively,  along  with  their  respective 
coefficients,  and  C^,  are  obtained  with  the  use  of  Equations  (A.  15)  -  (A.  19)  and  some 

tedious  but  straight-forward  algebra  and  trigonometric  identities. 
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Here  we  follow  Phillippidis  et  al.^  who  assumed  that  the  number  of  particles  impinging  on 
the  entrance  of  the  two  slits  is  one  at  any  given  instant.  Thus,  with  no  external  influence  the 
probability  of  finding  the  particle  emanating  from  either  slit  is  0.5.  When  the  number  of 

particles  impinging  on  the  slits  is  more  than  one,  i.e.,  n^as  defined  below.  Equation  (A.  15) 
represent  the  un-normalize  solution  for  both  slits, 

\|f  =  ^FpQ+5GQ0.  (B.l) 

We  no  proceed  to  show  that  this  wave  function  can  be  written  as  the  product  of  a  normalized 
wave  function,  X,  times  a  scaling  factor,  ,  to  establish  a  scaling  law  for  the  wave  function. 

The  boundary  conditions  at  the  surface  of  the  jet  slits  for  initial  particle  velocity,  ,  which 
is  in  the  axial  x-direction,  and  initial  number  of  particles,  ,  are  as  follows: 


yO(0,y)  =  0,  u(0,y)  =  0  , 

0  <  y  <  0.6 

and  0  >  y  >  -0.6 

(B.2) 

P(0,y)  =  /7,,  u(0,y)  =  u^, 

0.6  <  y  <  1.0 

and  -0.6>y>-1.0 

(B.3) 

/j(0,y)  =  0,  u(0,y) 

A 

o 

II 

and  y<-1.0 

(B.4) 

Since  the  Heisenberg  Uncertainty  Principle  limits  the  precision  with  which  one  can 
simultaneously  determine  the  position  and  velocity  of  a  particle  the  velocity  =  u(0,y)  must  be 

interpreted  as  an  average  velocity  equal  to  1,^(0, y)/ f^,y)  =  (2E/ . 

In  order  to  determine  we  follow  Phillippidis  et  al.^  The  expectation  value  of  the 
number  of  particles  per  unit  length  present  in  each  slit  is  given  by, 

= lo  y)  = lo  y)  ^y = ^y = ■  (b-5) 
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Here  the  boundary  conditions  on  n(x,y),  given  in  Appendix  A,  have  been  used  to  obtain  the 
final  result.  In  this  expression  the  scaling  factor  as  been  defined;  it  is 

=  (B.6) 


It  is  of  interest  to  note  that  when  —  or  ^  =  the  corresponding  wave  function  is 
essentially  the  normalized  wave  function  of  Phillippidis  et  al.^ 

An  expression  for  the  scaling  factor  G  is  obtained  directly  from  probability  current  density, 


h  ^  »  »\ 

J=p«=|^— J(y  Vv-^vvi/  j. 


(B.7) 


Since  the  particles  are  assumed  to  emanate  from  the  slits  parallel  to  the  x-axis  (B.7),  the  mean 
particle  velocity  in  axial  direction  is  given  by 


J  h  TT  -0  f2  +  Q  0 

_ ^  _ _ _ X _ X 

P  ^0  b[Q2+(G^/Fj0^ 


(B.8) 


where  the  subscript  x  indicates  the  partial  derivative  with  respect  to  x  as  introduced  in  Appendix 
A.  By  using  the  boundary  conditions  for  the  initial  mean  particle  velocity,  u^,  we  can 

determine  the  scaling  factor  G^  in  terms  of  7^  and  .  Substituting  the  dual  slit  wave  function. 

Equation  (A.  15),  and  its  component  parts.  Equations  (A.  16)  -  (A.  19),  into  (B8)  we  find, 

u(0,  y)  =  0,  0  <  y  <  0.6  and  0  >  y  >  -0.6,  (B.9.a) 

u(0,y)  =  u  =— —  —  (1  +  A),  0.6<y<1.0  and  -0.6>y>-1.0,  (B.9.b) 

“  e;  b 

u(0,  y)  =  0,  y>l  and  y<-1.0,  (B.9.c) 
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In  Equation  (B.9.b),  the  symbol  A  is  the  error  associated  with  the  approximation  occurring  in  the 


evaluation  of  the  Fredholm  integral  related  to  given  in  Appendix  A  for  selected 

values  of  N  used  in  the  calculation. 

Within  this  approximation,  we  write 


m  u  b 

_ ^  _  e  o 

F,~ 


(B.IO) 


However,  the  definition  of  the  Reynolds  number  given  in  the  text  is  N  =  u^b  /  (hji  /  m^) .  Thus 


it  should  be  now  obvious  that  we  can  recast  N  as  ratio  of  the  scaling  factors. 


G.  m  u„b 

N  =  —2.  = _ e  0 

El  h7i 


(B.ll) 


or  alternatively  by  Equation  (4. a),  is  equal  to  the  ratio  =  b  /  7i  with  the  deBroglie  wavelength 


Using  Equation  (B.l  1)  along  with  Eq,  (B.6)  in  Equation  (B.l)  yields  a  dual  slit  wave  function 
as  a  function  of  the  Reynolds  number, 

V(x,y,N)=  +  0  =  [2nj‘''(^Q  +  6  N  ©)  (B.12) 

Equation  (B.  12)  is  normalized  when  no =1/2.  That  is,  the  normalized  wave  function  is 

X(x,y,N)  =  (^Q  +  5N0).  (B.13) 

Substituting  (B.13)  back  into  (B.12)  results  in  the  scaling  law  for  the  wave  function  consistent 
with  Fitts,*’ 

V(x,y,N)  =  [2nj’''x(x,y,N).  (B.14) 
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FIGURE  CAPTIONS 


Figure  I. 
Figure  2. 


Figure  3. 


Figure  4. 
Figure  5. 


Figure  6. 


Figure  7. 


Dual  jet  slit  schematic  diagram. 

Probability  density  of  near  field  (Fresnel)  jet  streams  density  versus  axial  distance  as 
a  function  of  quantum  Reynolds  number.  The  magnitude  of  the  probability  density 
is  color  coded  in  the  key  next  to  each  figure. 

Probability  density  of  near  field  (Fresnel)  and  far  field  (Fraunhofer)  jet  streams 
density  versus  axial  distance  as  a  function  of  quantum  Reynolds  number.  The 
magnitude  of  the  probability  density  is  color  coded  in  the  key  next  to  each  figure. 
Probability  density  versus  axial  distance  as  a  function  of  quantum  Reynolds  number. 
Probability  density  versus  axial  distance  and  quantum  Reynolds  number.  La,  Lb,  and 
Lc  are  the  displacement  values,  on  the  jet  axis  of  propagation,  where  the  cluster 
initially  starts  to  grow,  for  the  cluster  maximum,  and  for  the  effective  end  of  the 
cluster,  respectively 

Parametric  study  of  the  effects  of  different  slit  separation  width  of  2b  on  the  cross- 
sectional  contours,  at  x  =  10,  of  probability  density  and  its  corresponding  interference 
pattern,  b  takes  on  the  values  of  b  =  0.8nm,  l.Onm,  L5nm,  2.0nm,  2.5nm,  and  3.0nm, 
corresponding  respectively  to  the  figure  panels  (a  -  f).  The  value  of  quantum 
Reynolds  number  is  N  =  10. 

Effects  of  a  change  of  quantum  Reynolds  number,  N,  on  interference  fringes  as  a 
function  of  axial  displacement,  x. 


-27- 


-28- 


Y-Vertical  Distance 


Fig.  2a  Chiu  et  al. 
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Fig.  2b  Chiu  et  al. 
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Fig.  2c  Chiu  et  al. 
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Fig.  2d  Chiu  et  al. 
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Fig.  2e  Chiu  et  al. 
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Fig.  2f  Chiu  et  al. 
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Fig.  3a  Chiu  et  al. 
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Fig.  3b  Chiu  et  al. 
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Fig.  3c  Chiu  et  al. 


-37- 


Y-Vertical  Distance 


Density 

Scale 


Fig.  3d  Chiu  et  al. 
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Fig.  3e  Chiu  et  al. 
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Fig.  3f  Chiu  et  al. 
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Fig.  4a  Chiu  et  al. 
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Fig.  4b  Chiu  et  al. 
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Fig.  4c  Chiu  et  al. 
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Fig.  4d  Chiu  et  al. 
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Fig.  4e  Chiu  et  al. 
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Fig.  5  Chiu  et  al. 
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Fig.  6a  Chiu  et  al. 
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Fig.  6b  Chiu  et  al. 
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Fig.  6c  Chiu  et  al. 


-49- 


N=10 


N=50 


N=100 


N=500 


X=5. 


X=15. 


X=30. 


J  ! 


■■ 


-^4 - ^ - D - S - 4 -  --_2  Q - ^ - - -g - - g - 


II 


,2  0  2  4 


-4'  *  "”-2  0  2;  4 

~r~ 


t 


Density 
I  Scale 


"^4  ^  D  2  4~ 


111 


1 


^2  0  2  4  ^  -2  0  2  4  ^  ^2  £P  "2  4 


■ 

1.00 

0.40 

0.16 

0.06 

0.03 

0.01 

Q  ^  4  -  _4  _2  0  '2  4  -:4  ^  £S  2  4' 


N=1000 


■ 


\2-  0  2 - A -  _4----2-  0  2  4  - -2'~-0“~‘2 - A 

A - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1- 

Y-Vertical 


Fig.  7  Chiu  et  al. 
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Theory  of  quamlets  of  quantum  systems 

H.H.  Chiu 

Space  science  and  Technology  Center 
National  Cheng  Kung  University 
Tainan,  70101  Taiwan 

Two  leading  quantum  mechanical  formalisms:  Schrodinger’s  wave  mechanics  and  the 
quantum  fluid  dynamics,  which  is  based  on  Hamilton- Jacoby  equation,  are  compared  to 
assess  their  quantum  mechanical  equivalence  between  the  linear  and  non-linear  fields.  A 
quantum  system  in  a  stationary  state,  described  by  the  linear  superposition  of 
eigenfunctions  corresponding  to  the  same  degeneracy  energy  level  is  found  to  be 
equivalent  to  the  mixture  of  the  “elemental  and  interactive  quantum  flowlets”,  termed 
“quamlets”,  each  of  which  is  constructed  by  the  bilinear  product  of  eigenfunctions  and 
their  complex  conjugates..  We  show  that  the  equivalence  between  two  formalisms  is  met 
if  each  quamlet  obeys  the  non-linear  equations  of  conservation  of  the  probability  density 
and  energy  when  the  global  quantum  flow  field,  which  consists  of  the  quamlet  mixture, 
satisfies  the  non-linear  equations  of  the  conservation  as  well  as  the  quantization  condition, 
as  in  the  quantum  theory.  An  example  of  the  quamlets  in  two-dimensional  flow  is  given 
to  illustrate  the  properties  and  the  compliance  of  the  conservations  of  the  elemental  and 
interactive  quamlets.  The  equivalence  principle  also  provides  the  guides  and  criteria 
required  for  the  analytical  and  numerical  simulation  of  complex  quantum  systems  by 
fluid  dynamic  approaches,  and  quantum  computational  fluid  dynamics. 

I.  INTRODUCTION 

Quantum  mechanics  has  been  successfully  applied  in  various  areas  of  quantum  physics: 
atoms,  molecules,  nuclear  physics,  superconductivity,  and  chemistry,  nevertheless  there 
have  been  number  of  attempts  proposed  to  offer  alternative  interpretations  based  on  the 
fluid  mechanical  description,  notably  by  Madelung  [1],  de  Broglie[2]  since  the  early 
developmental  stage  of  quantum  theory  and  later,  followed  by  Bohm’s  hidden  variables 
[4,5,6],  and  the  theory  of  stochastic  mechanical  approach  by  Nelson  and  others[6-16]. 
Various  aspects  of  fluid  dynamics  including  basic  formulation,  Takabayashi  [17-18],  Kan 
and  Griffin  [19],  vortex  motion,  Bialynicki-Birula  and  Sliwa  [20],  quantization  condition, 
[17,18],  Wallstrom[21,22].  The  motivations  and  the  physical  significance  underlying  the 
reformulation  of  quantum  mechanics  to  fluid  dynamic,  which  will  be  referred  to  as  QM- 
QFD  reformulation,  are  presented  by  Holland  [23].  One  of  the  basic  motivations  for  the 
QM-QFD  reformulation  is  prompted  by  the  criticism  on  the  lack  of  the  ontological 
interpretation  by  quantum  mechanical  theory,  which  renders  limited  physical  insight  for 
the  mechanisms  of  quantum  phenomena.  Unique  advantage  of  the  reformulation  is  the 
conversion  of  abstract  wave  physics  to  the  familiar  physics  of  the  fluid  flow  of  the 
equivalent  quantum  fluid,  characterized  by  the  density  and  velocity  of  the  probability 
fluid  whereby  the  quantum  phenomena  are  explained  by  the  principles  of  the 
conservation  of  probability  fluid  density,  mechanical  energies,  including  external 
potential,  mean  flow  kinetic  energy  and  quantum  potential  energy.  For  example,  in  the 
recent  analysis  based  on  QFD  formalism,  Chiu  [24],  identified  that  the  quantized  energy 
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of  hydrogen  atom  is  equal  to  the  kinetic  energy  of  the  asymptotic  diffusion  velocity  of  the 
electron  probability  fluid,  and  that  the  radial  dynamic  equilibrium  of  the  electron  is 
maintained  by  the  balance  between  the  quantum  potential  forces,  induced  by  the  gradient 
of  dilatation  and  diffusion  kinetic  energy,  with  the  central  Coulomb  force.  The  kinetic 
energy  of  the  vortex  induced  flow  is  provided  by  the  quantum  dilatation  energy 
associated  with  the  diffusion  in  polar  angular  direction,  ft  is  also  reported  that  the  zero- 
point  energy  in  a  simple  harmonic  oscillator  is  found  to  be  the  dilatation  energy 
associated  with  the  asymptotic  diffusion  of  the  quantum  probability  fluid  at  the  zero-point 
state.  The  advantages  of  QFD  are  not  limited  to  the  provision  of  ontological  interpretation 
of  the  quantum  phenomena,  as  described  by  few  examples  illustrated  above,  but  also  the 
access  to  well-established  fluid  mechanical  theorems  and  principles,  together  with  the 
direct  simulation  of  complex  quantum  systems  by  computational  fluid  dynamics. 
Needless  to  say  that  the  quantum  fluid  dynamics  characterization  must  be  in  compliance 
with  established  quantum  uncertainty  principles  associated  with  measurements  and  the 
inherent  limitation  set-forth  by  quantum  theory. 

The  basic  issue  of  the  QM-QFD  reformulation  is  the  quantum  mechanical  equivalence 
between  the  “linear”  wave  field  description  by  QM  and  the  “non-linear”  fluid  flow  fields 
described  by  QFD.  In  these  regards  we  ask:  How  the  equivalence  between  linear  wave 
field,  QM,  and  the  non-linear  field,  QFD,  is  ensured?  Specifically,  what  are  the  necessary 
and  sufficient  conditions  for  the  preservation  of  quantum  mechanical  equivalence?  These 
are  the  basic  motivations  of  the  present  study.  Some  fundamental  aspects  concerning  the 
equivalence  of  two  formalisms  are  described  below. 

Quantum  fluid  dynamic  formulation  is  built  on  the  basis  of  the  Hamilton-Jacobi,  H-J, 
transformation,  which  relates  the  Schrodinger’s  wave  function  with  the  density  and  a 
velocity  potential  of  the  quantum  flow  field.  The  flow  variables  obey  the  H-J  equations  of 
the  continuity  of  probability  density  and  Bernoulli’s  equations.  The  quantum  mechanical 
equivalence  between  the  quantum  mechanics  and  fluid  dynamical  descriptions,  are 
preserved  if  (1)  the  flow  variables  obey  the  H-J  equations  of  the  continuity  of  probability 
density  and  Bernoulli’s  equations  and,  (2)  The  quantum  flow  field  satisfy  the 
quantization  conditions,  as  in  quantum  theory. [  17,  18,  21,22]. 

The  objectives  of  this  study  are  to  examine  the  necessary  and  sufficient  conditions  for 
the  QM-QFD  equivalence  by  adopting  the  view  that  a  quantum  system,  characterized  by 
a  linear  superposition  of  eigenfunctions,  is  equivalent  to  the  mixture  of  quantum  flowlets, 
termed,  for  simplicity,  “quamlets”.  The  quamlets  are  the  basic  flowlets  constructed  by  the 
bilinear  products  of  the  wave  function  and  its  complex  conjugate.  In  this  approach,  the 
question  of  the  QM-QFD  equivalence  is  reduced  to  seek  the  necessary  and  sufficient 
conditions  that  the  sum  of  the  quamlets  jointly  creates  the  non-linear  field  that  is 
equivalent  to  the  linear  wave  field. 

Section  2  presents  the  description  of  quantum  diffusive  fluid  dynamics,  (QDFD),  to 
introduce  basic  theoretical  and  physical  terms,  which  will  lead  to  the  physical  concepts 
behind  the  quamlets:  Section  3  introduces  the  flow  variables  of  quamlets  and  their 
mixture,  together  with  the  relevant  algebraic  rules  relating  to  the  laws  of  superposition  of 
elemental  and  interactive  quamlets.  Section  4  describes  the  laws  of  the  conservation  of 
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quamlets,  general  properties  and  structure  of  quamlets.  Finally,  section  5  presents  the 
necessary  and  sufficient  conditions  for  the  conditions  of  the  QM-QFD  equivalence.  An 
example  is  given  to  illustrate  the  application  of  quamlets  in  the  analysis  of  quantum  flows 
in  the  selected  simple  systems. 

II.  QUANTUM  DIFFUSIVE  FLUID  DYNAMICS  OF  QUANTUM  SYSTEMS 

A  wave  function  of  a  quantum  system  in  the  presence  of  external  potential  is  governed  by 
the  Schrodinger’s  equation,  given  by 

=-T^vV+V.r  (2.1) 

at  2\m^ 


where  ^is  the  wave  function,  me  is  the  mass  of  the  particle,  h\s,  the  Planck’s  constant 
divided  by  27i  and  V  e  is  the  potential.  The  QFD  conservation  equations  of  the  probability 
density  and  the  mechanical  energy  are  formulated  by  the  Hamilton-Jacobi  transformation. 
Following  Bohm  [3]  we  first  express  the  wave  function  in  polar  form, 

y/  =  exp(im^^/h)  where  p  is  the  probability  density,  and  is  a  scalar  potential.  By 

inserting  the  latter  expression  for  {/dnto  Eq.  (2.1)  and  separating  the  result  into  imaginary 
and  real  parts  we  obtain  the  equations  for  the  flow  variables  p  and  .  The  imaginary 
and  real  part  gives  the  conservation  equation  of  probability  density  and  energy, 
respectively.  The  continuity  equation  probability  density  is  given  by, 

^  +  V.J  =  0  (2.2) 


where  the  probability  density  p  and  current  density  /  or  Lagrangian  density  are, 


P  =  ¥¥* 


(2.3) 


3  =  p\i  = 


2im 


e  J 


[y/^V  y/  -yN 


(2.4) 


The  energy  equation  or  alternatively,  the  quantum  Bernoulli’s  equation,  is  given  by 


dt 


f-1 

v2y 


V 

u»u+^-n 

m„ 


Inp 


(2.5) 


where  is  the  external  potential,  E  is  the  Bernoulli’s  constant,  and  the  fourth  term  is 
the  quantum  potential  energy 

The  velocity  u  is  defined  by  the  ratio  of  the  probability  current  flux  J  with  the  density  p 
as 
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(2.6) 


u  = 


r  ^ 


V^imey 


In^ 

\  V'  J 


The  velocity  field  is  irrotational  except  at  the  nodal  singularities  where  the  wave 
function  vanishes  and  the  vortices  are  present.  Under  certain  regularity  conditions  on  the 
initial  wave  function  and  the  external  potential,  the  particle  has  zero  probability  of 
entering  the  nodal  surface  from  outside  region,  Carlen  [25],  Nelson  [7],  and  Bemdl,  et  al. 
[26].  It  is  added  here  that  since  the  determination  of  the  current  density,  J,  requires  the 
simultaneous  measurements  of  position  and  velocity  of  particle,  which  is  not  in 
compliance  with  the  Heisenberg  uncertainty  principle,  hence  J  cannot  be  treated  as  the 
averaged  measured  particle  flux  at  the  point  in  space  and  time  as  interpreted  in  classical 
fluid  dynamics. 

The  quantum  potential  flinp ,  which  appears  in  the  right  hand  side  of  Eq.  (2.5)  is 
expressed  by. 


V^(ln/?)  +  filv(ln/?).V(lnp) 


(2.7) 


It  can  be  shown  that  the  following  two  alternative  forms  of  quantum  potentials,  11^  and 
,  are  equivalent  to  TIi^^ .  They  are  given  by, 


V  p  1 
P  ^2 


n  = 


Vp 


1/2 


P 


1/2 


(2.8) 

(2.9) 


The  last  form  appears  most  frequently  in  the  traditional  Madelung  formulation 

Grosh  &  Deb  [27].  However,  the  quantum  potential  expressed  in  the  form  Tli^^  appears 

ontologically  and  methodologically  superior  to  other  forms  in  treating  the  quamlet 
analysis,  which  will  be  performed  in  section  5.  We  will  adopt  the  Ilj^^form  of  the 

quantum  potential  and  for  notational  convenience  will  denoted  it  by  just  n . 


The  physical  significance  of  the  quantum  potential  has  been  illustrated  by 

decomposing  the  quantum  potential  into  the  sum  of  the  dilatation  and  the  kinetic  energy 
due  to  the  diffusion  of  the  probability  fluid  as  described  below.  The  gradient  induced 
diffusion  velocity  V  is  defined  by  Chiu  [24] 


^  n  ^ 


V2w,y 


V  In  yC>  =  -D 


Yr 

p 


(2.10) 
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where  9-  is  the  diffusion  flux  of  the  probability  fluid  and  the  quantity  h  /  2m^ ,  which  has 

the  same  physical  dimension  as  the  mass  diffusivity,  is  the  quantum  diffusivity,  D,  in 
analogy  with  conventional  terminology.  The  diffusion  velocity  is  the  negative  of  the 
osmotic  velocity,  Neson  [7],  which  is  the  mean  of  the  difference  between  the  forward  and 
backward  velocities  of  Brownian  motion  without  friction. 


In  addition  to  the  flow  velocity,  u  ,  Eq.  (2.6)  and  diffusion  velocity  V,  Eq.  (2.11),  we 
define  a  complex  velocity  w  as  follows, 


w  =  ^  (V  -  m)  =  -Z)V  In  y/ 


w*  =  ^  (V  +  m)  =  -DV  In  y/ 


(2.11) 

(2.12) 


The  physical  significance  of  the  first  term  of  the  quantum  potential,  )V^  (In  p) , 

in  Eq.  (2.8)  is  obtained  from  definition  of  the  diffusion  velocity  V.  By  taking  the 
divergence  of  Eq.  (2.1 1)  we  have, 


M  = 


V\\np)  = 


P 


=  -DS/>V 


(2.13) 


Hence,  the  first  term  in  the  quantum  potential  per  unit  mass  represents  the  specific 
dilatation  energy  of  diffusion  field,  M  =-D'V*V  due  to  the  divergence  of  the  quantum 
diffusion  velocity. 


Next,  based  on  the  Eq.  (2.1 1)  we  show  that  the  term  (ti^  /  4mJ(l/2)V(lnp)»V(lnp) ,  in 
Eq.  (2.8)  is  the  kinetic  energy  associated  with  the  quantum  diffusive  velocity. 


n  = 


Uwt/J 

.2, 

V(ln/7)»V(ln/?)  =  |E.E 


(2.14) 


Thus,  quantum  potential  is  the  sum  of  the  dilatation  energy  M  and  diffusion  kinetic 
energy,  Q 


=  M 


+  Q  = 


rn 

-DV»V  + 

v*v. 

At  steady  state,  the  Bernoulli’s  equation  is  expressed  by 

f4Vu= + vq + £>v  •  V = [(A)  _  (A)] 


(2.15) 


(2.16) 


The  above  Bernoulli’s  equation,  as  we  shall  show  later,  can  be  rewritten  as  follows. 
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(2.17) 


Thus  a  wave  function  at  a  specific  quantum  state  can  be  represented  by  the  quanum 
flow  field,  whose  flow  variables:  density  and  velocities  are  obtained  from  the 
transformation  relations  and  obeys  the  continuity  equation  (2.2)  and  the  Bernoulli’s 
equation  (2.5). 

III.  QUAMLETS  AND  ALGEBRAIC  RULES  OF  QUAMLET  MIXTURE 

A  quantum  system  in  a  stationary  state  is,  in  general,  described  by  a  linear  superposition 
of  eigenfunctions  y/  a(E,  x ),  corresponding  to  the  same  degenerate  energy  level. 


N  -iEt/ 

'F(E,x,t)  =  E^„(E,x)e  a  =  l,2,3....N  (3.1) 

The  wave  function  function  is  expressed  by  \[/  „  =  Aa  (])  a,  where  Aa  is  a  constant  and 
(j)  a  is  a-th  eigenfunction  of  a  degenerate  energy  level. 

A.  Elemental  and  interactive  quamlet 

We  first  investigate  the  flow  field  corresponding  to  a  stationary  state  described  by  a 
single  wave  function,  y/a.  The  product  of  single  wave  function,  y/a,  with  its  complex 
conjugate  is  defined  as  an  “elemental  flowlef’,  which  is  termed,  in  short,  a  aa*-ih 
“quamlef’.  The  elemental  flow  field  is  characterized  by  a  set  of  quantum  flow  variables: 
the  density,  paa*  and  velocity  Uaa*  of  aq  elemental  quamlet  flow  field,  given  by. 


^aa*  T  a 


h 


2im„ 


-Vln 


v^g*  y 


-h 


Wa  =(V2)(V„„,-u„„.)  =  ^V(lnyJ 

fl 

w„*=(l/2)(V„„,-iu„„,)  =  ^V(lnv„,) 


(3.2) 

(3.3) 

(3.4) 

(3.5) 


Note  that  the  complex  velocity  Wa  does  not  depend  on  Iny/a*,  hence  the  second  index 
a^does  not  appear  in  the  complex  velocity.  Being  a  physical  entity  describing  a  fluid 
flow  field,  each  elemental  quamlet  satisfies  the  continuity  and  Bernoulli’s  equations,  i.e., 

%^  +  V'(p„u„)  =  0  (3.6) 
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Now  we  turn  our  attention  to  the  quamlets  representation  of  the  complete  flow  field, 
corresponding  to  the  wave  field  consists  of  the  linear  superposition  of  the  degenerate 
wave  functions.  Like  quantum  wave  formalism  wherein  a  superposition  of  the  degenerate 
eigenfunction  gives  a  complete  field,  one  is  led  to  conjecture  that  a  complex  quantum 
system  can  be  represented  by  the  superposition  of  “quamlets”  :  elemental  quamlets, ,  and 
the  interactive  quamlets,  which  are  created  by  the  bi-linear  interaction  between  a  pair  of 
eignefunctions:  aP*  where  the  latter  quamlet  is  the  complex  conjugate  of  the 

former.  Detailed  mathematical  characterization  of  the  flow  filed  by  quamlets  is  discussed 
in  following. 


B.  Quamlets  flow  variables 

/?  =  ¥pS\  =  V p 

^  P  ^  a  ^  P 


(3.8) 


In  the  above  expression  the  double  summation,  Za  (2p* )  is  carried  out  first  with  respect 
to  the  second  index,  |3*,  followed  by  the  first  index,  a,  and  S(p  ^  «)  represents  the  sum  of 
P  for  the  entire  set  of  degenerate  states,  except  a-th  eigenstate.  There  are  two  types  of 
densities:  real  valued  density,  paa*  and  complex  densities,  pap*  and  pap*  as  follows, 

P  ^a*  Paa*  ^p*  Pap*  ^a*  ^p  Pa*p 

(3.9) 

a*  p  *  a  5*  P 


where,  the  quantity  paa*  =  Pa*a  =  Va  Va* ,  is  the  density  of  the  elemental  quamlets  and 
Pa*p  is  the  complex  density  of  an  interactive  quamlet  field,  created  by  the  product  of  two 
different  eigenfunctions and  (pa*p)=y/a*Wp  =  (Pap*)*-^^^  elemental 
quamlets  obey  the  linear  superposition,  whereas,  the  interactive  quamlets  comply  with 
special  algebraic  rules  involving  the  summation  of  binomial  products.  The  addition  rule 
involving  the  linear  superposition  of  elemental  quamlets  and  summation  of  binomial 
products,  for  in  interactive  quamlets  is  the  consequence  of  the  non-linearity  of  fluid 
dynamics. 


C.  Algebraic  rules  of  quamlet  mixture 

Following  the  definition  of  the  density  of  quamlet  mixture,  we  define  the  quamlets 
concentration  of  aP*,  by  Zap,  expressed  by. 


PaB* 

V-  " - 

P 


(3.10) 


By  virtue  of  (3.9),  the  sum  of  the  concentrations  of  all  the  quamlets,  is  unity. 


(3.11) 


Pap* 

P  ) 


^ap*  )  1 


Alternatively,  we  write 

^ap*  )  ^a*  (^aa*  )  ^p*  ^ap*  ^a*  ^p  ^a*p  ^ 

P*7sa*  P?ia 


(3.12) 


Next,  by  using  the  definition  of  the  quamlet  concentration,  we  show  the  following  rules 
of  the  multiplication  of  the  concentrations  Zap*  and  Za*y, 

^ ap*^ a*Y  ~  I P  ~Paa*Pp*ylP  ~^aa*^p*y  (3.13) 

^a*p^aT*  ~  (Va*Vp)(H^aH^Y*)/P  ~  Paa*PpT*/P  ~  ^aa*^pY*  (3-14) 

Next,  by  (3.13)  and  (3.14),  we  obtain  the  rules  of  the  sum  of  the  product  of  the 
concentrations,  with  respect  to  /?  and  y  *,  as  follows. 


2jp*(2jy  ^aP*^a*Y 


z  z 

y*  a*p  ay* 


)  ^aa* 

)  =  Z 


Sb*(S..z 

a*a 


Ib(I.,.Z 


)  Z'oifx* 

^  =  z 

py*  /  a*a 


(3.15) 

(3.16) 


These  multiplication  rules  of  quamlets  concentration,  (3.13)  and  (3.14)  are  used  in  the 
analysis,  which  involves  the  algebraic  operations,  addition  and  multiplication  of  the 
quamlet  concentrations,  as  described  in  the  latter  part  of  the  analysis. 

D.  Quamlet  fluxes  and  velocities 

The  flux  of  probability  fluid,  /  =  /?«,  of  the  global  flow  field  is  expressed  by  the  sum 
of  the  flux  of  elemental,  J a  a*,  and  that  of  interactive  quamlets,  Jap  follows. 


j  =  y  j 

^  aa^ 


a/3^  /3  ^  a^/3 


(3.17) 


where,  fluxes  of  elemental  and  interactive  quamlets  are  given  by 

Jaa*  =  Paa*'^a  =  iP^aa  ( )  (3.1 8) 

Jap.  =  ipz„p.  (1/2)  [(V„„,  -  Vpp, )  -  i(u„„,  +  Upp, )]  =  ipz„p,  (w„  -  Wp  J  (3.1 9) 
Ja.p  =  -ipZa.p  (1/2)  [(V„„*  -  VppJ  +  i(u„„,  +Uppj]  =  -ipz„,p(w„,  -  Wp)  (3.20) 
j=ip2:„2:„.2:z„„.(w„-w„.) 


+  ip 


V  V  7 

^p*  ^ap* 


P*^a* 


(w„-wp.)-i:„.2:p  z„.p(w„,-wp) 

P^a 


(3.21) 
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By  dividing  the  /-flux  by  the  fluid  density,  we  obtain  the  velocity  of  the  complete  flow 
in  terms  of  the  quamlet  concentration  and  complex  velocity  as  follows 

u = j/p = i  ) 

+  i  Ip.  z„p,  ( w„  -  Wp, )  - 1„.  Ip  z„.p  ( -  Wp )] 

We  define  the  mean  velocities  of  quamlets  as  follows, 

Uaa*  ~  iZ  aa*  -Wa*)  Uap*  ~  iZ  ap*  -Wp*)  lla*p  ~  iZ  c(*p  (Wj^*  "Wp)  (3.23) 


One  important  aspect  of  the  equivalence  between  the  quantum  flow  description  and 
quantum  mechanics  is  the  quantization  condition  must  be  met  in  the  former  flow  field, 
i.e.. 


u»dr  =  2n:m 


c 


which  is  rewritten  by 

^ a*  ) 

C 


(3.24) 


where,  m  is  an  integer  and  the  subscript  c  is  any  closed  counter. 
The  flux  of  diffusion  of  probability  fluid  3  =  pV  ,  is  given  by 


i  =  y  y  1  -i-y  y  i  -i-y  y  i 

^  ^a*  ^aa*  ^  ^ p*  '^ap*  ^  ^a*  ^p  ^a*p 


a 


where,  fluxes  of  elemental  and  interactive  quamlet  are  given  by 

'^aa*  ~  Paa*Xx  ~  PZaa*(''^a  "^a*) 

=  PZap-(^)[(Vaa*  +  Vpp.)  "  i(u„„  +  Upp.)]  =  pZ„p.(w„  +  Wp.) 

3a*p  =  Pz„*p(|)[(V„„.  +  Vpp.)-i(u„„.  +Upp.)]  =  pz„.p(w„.  +  Wp) 

3  =  plz„„.(w„  w„.)  +  p[I„  Ip.  z„p.(w„  +  Wp.) 

+  Za*Zp  Z„.p(W„.+Wp)] 


The  diffusion  velocity  of  the  complete  flow  field  is  given  by 


(3.25) 


(3.26) 

(3.27) 

(3.28) 

(3.29) 

(3.30) 


Inm 
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(3.31) 


V  =  -  =  (w„  +  w„. )  +  [I„  Ip.  z„p.  (w„  +  Wp.  ) 

P  a%p* 

+  I„.Ip  Z„.p(w„.+Wp)] 

We  define  the  diffusion  velocities  of  quamlets  as  follows, 

Vaa‘=Zaa‘(Wa  +  Wa*),  Vap‘ =  Z  ap‘ (Wa  +  Wp*),  Va*p  =  Z a*p  (Wa*  +  Wp)  (3.32) 
Similarly,  the  probability  flux  based  on  the  complex  velocity,  w,  is  expressed  by 


X  =  (i)(3-iJ)  (3.33) 

Zaa*  =  Paa*W„  =  pZ„„,W„  (3.34) 

5CaP*  ~  P^ap*''^a  Xa*p  ~  P^a*p''^a*  (3.35) 

X  =  pI„I„.z„„*w„+p[Z„Sp*  z„p*w„+Z„*ZpZ„*pW„*]  (3.36) 

a%p'  “’^P 

The  complex  velocity  w  of  the  complete  field  is. 


w  =  -  ( V  - iu)  =  z„„*w„  +  p[Z„  Ip.  ^  z„p*w„ 

V  J  a%P* 

^a*  ^p  ^a*p  ^a*  ] 

The  complex  velocities  of  quamlets  are  given  by, 

Waa*  ~  Z  aa*  Wq[  W^p*  ~  Z  ap*  Wa*p  ~  Z  ^.p  Wq[*  (3.38) 

IV.  PRINCIPLE  OF  EQUIVALENCE 

The  equivalence  between  the  quantum  flow  field,  described  by  (3.9),  (3.22)  and  the 
quantum  mechanical  characterization  by  wave  function  (3.1)  requires  that  the  overall 
flow  field,  represented  by  the  sum  of  all  the  qaumlets,  must  obeys  the  conservation 
equations  (2.2)  and  (2.5),  and  moreover  the  global  flow  meet  the  quantization  condition 
as  described  by  Eq.  (3.24) 


A.  Continuity  equation 

We  shall  first  examine  the  necessary  condition  under  which  ensures  the  quamlet  mixture 
obeys  the  law  of  continuity  of  probability  density. 

The  application  of  a  divergence  operator  on  the  mass  flux,  Eq.  (3.17),  gives 
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V  •  pu  =  V  .  ( W„  -  ) + 

i  2:p.V.p„p.(w„-Wp.)-Z„.  ZpV.p„.p(w„.+Wp) 
aVp*  a^p  J 

At  stationary  state,  the  continuity  requires  that  the  above  expression  must  vanishes,  i.e., 
V-pu  =0. 

To  prove  that  the  above  expression  (4.1)  is  zero,  we  first  examine  the  divergence  of  the 
flux  of  an  elemental  quamlet  of  pa  a, 

V  •Pa„*«aa«  =  1  ^  V  •  p„„.  (w„  -  w„.)  =  AD  Z  V(\^„\^„.)(ln  \j/„  -  In 

=  D{(X1/„V  V„*  -  VX|/a*VX|/„)-(¥a*'^  V„  -  J  =  D(\1/„V  V„.  -  Va») 

(4.2) 

The  last  expression,  in  the  above  equation,  can  be  further  deduced  as  follows.  First,  we 
multiply  the  Schrodinger’s  equation  (2.1)  of  the  degenerate  eigenfunction  y/y*  by  y/y,By 
carrying  out  the  same  step  as  above  except  that  we  interchange  the  role  of  y/yhy  y/y*. 
Next,  by  subtracting  the  two  resulting  equations  from  each  other,  we  have 

-^(¥a^'¥a*-¥a*^^¥a)  =  -V,){¥ a¥ a* -¥  a¥ a*)  =  0  (4-3) 

2m^  V  > 

Thus,  the  elemental  quamlet  aa*  satisfies  the  law  of  continuity,  as  expected,  i.e., 

^•[Paa*('^a-W„.)]  =  0  (4.4) 

Next,  we  calculate  the  divergence  of  the  interactive  quamlet,  afi*,  represented  by  the 
second  term  of  Eq.  (4.1), 

•  Pap*  (w„  -  Wp. )  =  iDV  •  (y „i^p. )(V  In  I^p.  -  V  In  y  J 

=  iD{(Va'^Vp*  -V'fp.Vl^J-(Vp.VVa  -V'fp.V^a)  =  iD(l^„VVp* -Vp*VVa) 

(4.5) 

By  the  similar  procedure  used  in  the  above  analysis,  we  can  show  that  the  last  term  of  the 
above  equation  vanishes,  i.e., 

~ ¥py^¥a )  =  (^ -  VX¥a¥p*  - ¥p*¥a)  =  0  (4-6) 

Thus,  the  quamlet,  a  ^=1^ ,  satifies  the  continuity  of  the  probability  fluid  as  follows. 
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iV*Pap*(Wa-Wp.)  =  0 


(4.7) 


By  decomposing  the  above  equation  into  real  and  imaginary  parts,  we  get, 


n  f  n 

T  -  V.[pq„p.(u„„,-Uppj]  =  0 

I  \  I 


where,  i^ap  and  rjap  are  the  real  and  imaginary  part  of  given  by, 

PCap=(Op+®a®p) 

P^ap=(Va-<t>„®p) 


in  which  and  cox  are  the  real  and  imaginary  parts  of  xj/x  respectively. 
Similarly  the  continuity  of  a*P  is  given  by. 


(4.10) 


iV*Pa*p(Wa.-Wp)  =  0 


(l/2)/v{/?4'„^(v„„,  -v^^J]-(1/2)V{/?;7„^(u„„,  +u^^j]  =  0 


(4.11) 

(4.12) 


Thus  by  comparing  two  sets  of  continuity  equations  for  (4.8)  and  a*P,  (4.12)  we 
find  that  the  probability  of  the  two  quamlets  are  satisfied  provided. 


"^•\_PCapi^aa*-^ppS\  =  ^ 


(4.13) 

(4.14) 


We  conclude  that  the  necessary  condition  for  the  global  flow  to  obey  the  continuity  is 
that  each  quamlet  satisfies  independently  the  continuity  equation  of  each  quamlet  flow. 

C.  Bernoulli’s  equation  of  quamlets 

We  next  investigate  the  necessary  condition  required  for  the  quamlet  mixture  flow  to 
obey  the  Bernoulli’s  equation. 

The  dilatation  term,  D  V  »v  ,  represented  by  the  second  term  on  the  right  hand  side  of 
the  above  equation,  can  be  reformulated  in  the  following  form. 


D\  ]  \  D\ 

DV‘V  =  DV*  —  =  —  V»3  +  T)3V  —  =  —  +  (4.15) 

\P)  \P)  \P)  \P) 


Substitution  of  (4.15)  into  (2.16),  gives, 
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(4.16) 


i')(B=+V’)+-V.3  =  [(— )-(i^)] 


m„  m„ 


By  using  Eqs.(3.4)  (4-15)  we  re-express  the  above  equation  in  the  following  form, 
2(^ww*)- 


V. 

\P) 

_ 

Z  «  Z  Z  a*  Z /?Z3„.^ 

p^a  p^a 


=  [(— )-(^)] 


(4.18) 


The  right  hand  side  of  the  above  expression  can  be  recast  in  the  form  representing  the 
net  energy  (E-V)/me  contributed  by  the  sum  of  all  the  quamlets.  This  is  made  by  the 
following  two  steps.  First,  we  multiply  the  latter  sum  of  energies  on  the  right  hand  side  of 
the  equation.by  the  following  expression,  Eq.  (3.12  ). 

Z„  Z„.  )  +  z„  Up,  {z^p, )  +  Z„.  Z^  {z^,p  )  =  1  (4.19) 


Next  by  substituting  (3.37)  and  (3.26)  into  Eq.  (4.18)  and  by  expanding  each  term 
appears  in  the  above  expression  through  the  use  of  (3.13)  to  (3.16),  we  have 


(d 


2  {z„.  Tp  [z^  )  z,.  (zp^, ) 

+ z„  Zp,  y^lpap*  (w„ + Wp, )] + z„.  Zp  (w„. + )] 


-w^ }  +  J  {z„  Z„.  ( w„  +  )]} 


E  V 

[(  )  (  )]  I  'Zif^  Z(j..[!  'Z^aa*  Z^(x  'Zip,  -|-  'Zif^,  'Zip  Z^^p  j" 

m„  m„  6^  K  H  H ) 


(4.20) 


The  objective  of  the  remaining  part  of  the  analysis  is  to  prove  the  equality  of  the  left 
hand  side  of  Eq.(4.20)  with  that  of  the  right  handside.  The  detail  analysis  of  this  assertion 
is  given  below. 

Following  the  quamlet  analysis  for  the  continuity  equation,  we  examine  the 
conservation  of  the  energy  for  three  types  of  quamlets:  aa*,  ap*  ,and  a*p. 

First,  by  the  application  of  quamlet  multiplication  rules  (3.13)  and  (3.14),  we  have 


By  substituting  the  above  expressions  in  the  kinetic  energy  of  mean  motion  appear  in  Eq. 
(4.20)  we  have, 
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(4.22) 


(w„.w^,)  +  2{Z„,  z„,^  [Z^  Z,,(z^,,)](w„,.w^)} 


Secondly,  we  write  the  following  summation  rules  based  on  (3.1 1) , 

Z,,Z,(z,,,)=l  ,  Z,.Z,(z,*^)=l  ,  Z,Z,.(z,,*)=l  (4.23) 

Thus,  the  total  kinetic  energy  is  reduced  to  the  sum  of  the  energy  of  elemental 
quamlets  and  interactive  quamlets  as  follows, 

2  Z„  Z„*  (z„„* )(w„  )  +  2  Z„  Z^*  (z„^* )(w„ . )  +  2  Z„.  Z(z^  )(w„. •  )  (4.24) 


By  substituting  above  expression  into  (4.22)  and  collecting  the  terms  associated  with 
each  quamlet  group  as  follows  we  have  the  following  global  Bernoulli’s  equation, 


2Z„Z„.(z,„.)(w„.w„.)+  —  {Z„Z„.V{/?„„.(w„+w„.)]} 


\Pj 


+2Z„Z^.(z„^.)(w„*w^.)+  —  {Z,Z^.V.[/7„^.(w„+w^.)_ 


kPj 


+2Z„.Z^.(z^)(w„.*w^)+  —  {Z,.Z^V.[p„.^(w„.+w^) 


\Pj 


E  y. 


[(— )  -  (—)]{!«  Z„«.  +  I/?,  ^ap*  +  ^p  ^a*p  } 


(4.25) 


m„  m 


We  shall  examine  the  energy  balance  of  the  quamlet  mixture  as  follows. 


First  of  all,  by  collecting  all  the  terms  associated  with  ora terms  on  the  left  hand  side  of 
(4.25), 


2Z„Z„.(z„„,)(w„*w„.)+  -  {Z„Z„.V{/7„„.(w„+w„.)]} 

\P  J 


(4.26) 


By  using  (3.4)  and  (3.5),  we  re-express  the  dilatation  energy  given  by  the  second  term  in 
the  above  equation,  as  follows. 


u. 


Z„  Z„.  ( w„  +  )  =  -(^) z„  Z„.  [v {\ny/^+  In )]}  (4.27) 


By  expanding  the  gradients  of  logarithmic  terms  in  (4.27)  and  by  the  applications  of 
divergence  operator,  we  have. 
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-D" 


r((^„.VV„)  +  ((^„VV„.)]}-2{Z„  Z„.  (w„.w„.)}  (4.28) 


The  terms  in  the  first  curly  bracket  can  be  further  reduced,  as  follows.  Multiply  the 
Schrodinger  equation,  governing  y/a*  by  y/a  and  by  repeating  the  same  procedure,  except 
by  interchanging  the  role  of  y/a  with  y/a*  ■  Next,  by  adding  the  resulting  equations  we 
obtain, 


4m. 


^  1  ^ 


v2mey 


v2niey 


VV„ 


¥a¥aA 

-(7  \ 

\^aa^  ) 

\  P  J 

^ej 

(4.29) 


Thus,  the  dilatation  energy  of  all  the  elemental  quamlets,  aa*,  is 


\P/ 


( W„  +  )  =  z„  Z„.  )  [(— )  -  (-^)]  -  2  z„  Z„.  ( w„  .w„. ) 


(4.30) 

The  sum  of  the  kinetic  energy  and  dilatational  energy  of  elemental  quamlets  is 
obtained  by  substituting  (4.30)  into  (4.26)  as  follows, 

(  Yy\ 

2Z„Z„.(z„„.)(w„.w„.)+  -  {Z„Z„.V.[p„„.(w„+w„.)]} 

\P  J 

=  2  Z„  Z„.  (z„„* ) ( w„  .w„. )  +  Z„  Z„.[(— )  -  (-^)]  z„„.-  2 Z„  Z„.  (z„„. ) ( w„  .w„, ) 


(4.31) 


Following  the  similar  algebraic  procedure,  we  reduce  the  dilatation  energy  of  the 
interactive  quamlets  a/3*  and  a*^  appear  in  (4.25)  as  follows, 


kP/ 


[p 


z„  z^.  )] = z„  2:^.[(— )  -  i-^)]  -  2  z„  z^.  z„^.  (w„  ) 

(4.32) 

Z„*  Z^  ^{pa*p  (  )]  =  Z„.  Z^  [(— )  -  (-^)]  2  Z„.  Z^  T^a*p  a*  p  ) 

W-e 

(4.33) 
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The  dilatation  of  interactive  quamlets  consists  of  two  parts:  the  quantized  energy  of  the 
quamlets,  represented  by  the  first  term  appears  on  the  right  hand  side  of  the  above 
expressions,  and  the  negative  of  the  kinetic  energy  of  the  quamlets,  expressed  by  the 
second  term. 


Thus  by  adding  (4.30)  (4.32)  and  (4.33)  we  obtain  the  dilatation  energy  as  follows. 


(  d\ 

+  — 

\p) 


+ 


yp) 

E  V 

’[(  )~(  ’^a*p\ 

m„  m„  ^  ^  A- ; 


-2  ( w„  )  -  2  )  -  2  ) 


(4.34) 


Finally  by  substituting  all  the  dilation  energy  obtained  in  (4.36)  into  the  left  side  of  the 
global  Bernoulli’s  equation  (4.25),  we  show  that  the  left  hand  side  of  the  equation  is 
indeed  equal  to  the  right  hand  side,  i.e.  [(E/me )  -  (Vg/me)] 


)  +  2  z„  Z^,  (z„^, ) (w„  )  +  2  Z„,  Z^  (z„,^ ) ( ) 

.  I  Z„  Z„,  (  Z^,  (  T-aP*  )  +  Z„,  Z^  (  Z„,^  )|  [( )  -  (— )] 

-2 Z„  Z„,  (z„„, ) ( w„ )  -  2 Z„  Up.  (z^p. ) (w„ ) 

-2Z„.  Z^(z„.^)(w„..w^) 

=  {Z„  Z„*  (z„„. )  +  Z„  Zp.  (z^p. )  +  Z„*  Zp  (z^.p )}  [(— )  -  (-^)]  (4.35) 

/fie 

F  V 
=[(— )-(^)] 

We  have  proved  that  the  quamlet  mixture  satisfies  the  global  Bernoulli’s  equation. 
This  is  one  of  the  basic  criteria  for  the  QM-QFD  equivalence.  We  shall  show  that  the 
necessary  condition  to  meet  the  criterion  is  that  each  quamlet  must  be  the  solution  of  the 
Bernoulli’s  equation  of  individual  quamlet,  as  described  in  the  following  subsections. 


Energy  balance  of  elemental  quamlets:  aa* 

By  collecting  the  terms  associated  with  terms  in  the  left  hand  side  of  Eq.  (4.35),  we 
have, 
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2  )  +  f— 1  {l„  ^{paa*  ( W„  +  W„. )]} 

\P  J 

E  V 

=  2  z„  Z„.  (z„„. ) ( w„  .W„, )  +  z„  (z„„* )[(— )  -  i-^)] 

-2E„  (z«a.)('^a*W„.)  =  E„  {^aa^\i— )-{—)]  (4-36) 

m„ 

e  e 

The  degenerate  wave  functions  are  linearly  independent  solution  of  the  Schrodinger’s 
equation,  hence  if  the  wave  functions  are  continuous  up  to  the  first  derivatives,  then  the 
necessary  condition  for  the  validity  of  the  Eq.  (4.36)  is  that  each  elemental  quamlet, 
aa*  satisfies  the  following  the  “elemental  quamlet  Bemoullis  equation”  . 

2  )  +  (—]  (w„  +  )]  =  [(— )  -  (-^)]Z„„,  (4.37) 

VpJ 

Balance  of  interactive  quamlets:  af*  and  a*P 

Following  the  same  mathematical  steps,  we  find  that  the  sum  of  interactive  quamlets 
obeys  the  following  equation, 

+2 1,  Up,  [z^p, ) (w„  {S„  Tp,  (w„  +  )]} 

E  V 

=  2  ) ( w„  fLp,  (z^p^ ) [(— )  -  (^)]  -2  (z^p, ) ( W„  ) 

Me 

=  Z„Z,.(z„,.)[(— )-(ti)]  (4.38) 

e  e 

If  the  degenerate  wave  functions  for  Grand  ,  which  are  linearly  independent 
solutions  of  the  wave  equation  of  degenerate  energy  level,  are  continuous  up  to  the  first 
order  derivative,  then  the  above  equality  hold  when  each  interactive  quamlet, 
satisfies  the  interactive  quamlet  Bernoulli’s  equation. 

( r 

\PJ 

E  V 

=  2  (z„^. )  ( w„  )  [(— )  -  (^)]  -  2  (z„^, )  ( w„  ) 

T7l^  Tfl^ 

=  (v^.)[(— )-(—)]  (4.39) 

'  m  m 
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Thus,  we  conclude  that  the  quamlet,  aP*,  satisfies  the  Bernoulli’s  equation  of  . 


)(".  +  "/>• )]  =  )[  A  (4.40) 

By  the  similar  method  described  above,  we  show  that  satisfy  the  following  equation. 
a*P  quamlets 

°  j  (w„.  +  )]  =  )[(^)  -  (^)]  (4.41) 

By  expressing  the  above  equations  in  terms  of  the  mean  velocity,  and  the  diffusion 
velocity,  i.e.,  real  valued  velocities  we  arrive  at  the  following  Bernoulli’s  equations  for 
three  types  of  quamlets 

Bernoulli ’s  equation  of  aa*  elemental  quamlets 


Wz  ( 

2  act*  V 

y 


U  *11  +\  •\ 

act*  aa*  aa*  aa- 


)+f-l  V-(,£>„.v„. )  =  z„.[(— )  -  (i^)]  (4.42) 

VpJ 


The  real  part  of  the  Bernoulli’s  equation  ap*  is  given  by, 


[P^ap)[ 


^ aa*^ PP*  ^aa 


^pp*)  +  PVap[ 


^aa*^ PP*  ~^^aa*^ pp* 


^  D  ^ 

-  —  PP*)^  PVap{^aa*~'^  PP*)\ 

\P  J 


t/»C,.)[(— )-(—)] 


(4.43) 


The  imaginary  parts  is 

^  aa*^  PP*  ^aa*^PP*^  P^ap  {^aa*^  PP*  ^aa*^y0y0*)J 

^  D  ^ 

'^a^p*'^\p^ap*{^^  aa*'^^  PP*)~  P^ap{}^aa*~^pp*^ 

\P  J 

=  Z.Z,.(/jC^.)[(— )-(—)]  (4.44) 

^  '  m  m 

e  e 

It  is  concluded  that  the  necessary  and  sufficient  conditions  for  the  equivalence  between 
the  quantum  flow  field  with  the  wave  field  are  that  each  component  of  quamlet  satisfies 
the  equations  of  continuity  of  probability,  Eqs.(4.4)  (4.7),  (4.11),  the  Bernoulli’s 
equations  (4.40)  (  4.43)  and  (4.44)  and  the  global  flow  field  satisfies  the  quantization 
condition,  (3.24). 
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E.  Quamlet  energy  balance 


We  consider  electrons  moving  with  energy  E  in  two-dimensional  space  bounded  by 
parallel  plates  with  the  separation  between  two  plates  being  equal  to  b.  The  plates  extend 
from  X  =  0  to  infinity. 

Since  the  electrons  are  taken  as  free  particles  of  energy  E,  we  assume  their  wave 
functions  are  stationary  state  plane  waves  and  are  of  the  form  T  =  y/{x,  y)  exp(iEt  /  K) . 
The  time-dependent  Schrddinger  equation  of  motion  governing  the  electrons  is, 

(4.45) 

dt  2m^ 


where  h  is  the  Planck’s  constant  divided  by  2k,  and  is  the  rest  mass  of  an  electron. 
The  probability  density,  [].  of  the  electron  is  then, 

P  =  WW\  (4.46) 


where  \f/*  is  the  complex  conjugate  of  the  wave  function  y/ 

Time  independent  equation,  governing  the  stationary  wave  function,  y/,  is  obtained  by 
substituting  T  =  y/{x,y)  Qxp{iEtl  /t)into  Eq.  (4.45),  as  follows, 

VV(x,y)  +  A^V(^,T)  =  0,  (4.47) 


where,  V^=  d^/dx^+  d^/dy  ^ ,  x  =  nx  /h,  y  =  Ky_ /h  and  =  (2me/lf)( b/nf  E,  in  which 
b  is  the  plate  separation,  and  N  is  given  by. 


2m^E  b^ 

*2  ~T  ’ 
n  K 


(4.48) 


which  may  also  be  written  in  terms  of  the  deBroglie  wavelength,  /ijg  /  2m^E , 


\^dBj 


(4.49) 


where  =  b  /  tt  .  Alternatively,  N  can  be  view  as  a  quantum  Reynolds  number  that  is 

defined  by  the  ratio  of  typical  inertial  force  to  quantum  viscous  force.  Thus 
N  =  uj}l{fi7r  I  m^)  =  uJdIVq  ,  where  nie  The  solution  of  the  particle  field  is  constructed  by 

the  linear  superposition  of  the  wave  functions,  which  satisfies  Eq.  (4.47),  as  follows. 
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y/{E,x)  =  Y.A^y/^{E,x)  a  =  1,2,3,....  (4.50) 

The  wave  funetion  of  a-ih  degenerate  state  is  given  by, 

y^a  -  \  exp[-/ay-i(A^^  (4.51) 

y/p  =  exp[-i;^y-/(7V^  -y^")^x]  (4.52) 

The  above  wave  function  for  cr-th  state  creates  an  elemental  quamlet  a-a*  whose  density 
and  velocities  are  given  by, 

Paa  expHay-/(fV"  -  «" )^ x]  exp[/ay  + /(fV"  -a")^x]  =  4"  (4.53) 

Pap'  =y^a¥p-  cxp {-/(«- ;5)y - /(fV"  - {N^  - x]]  (4.54) 

=  -(hTT /  2m^b)V([-iay - i(N^  -a^)^x]  =  -{hTr / 2m^b)[-i{aey  +(N^  -a^)^e^]  (4.55) 
w^.  =  -(hTT  /  2m^b)V{[ij3y  +  i(N^  -  ^'  )>^  x]  =  -(Htt  /  2m^b)[i{pty  +  {N^  -  )^  e  J  (4.56) 

w„-w^.  =-(/?;r/2m^6)(H(a  +  y9)e^-/[(iV'-a')^+(iV'-y5')^eJ  (4.57) 

u  ,  =(h/  2imJp)V  ln(^^  /  y/J)  =  {h/  m^b)V[-ay  -  (N^  -  a^)^x] 

(4.58) 

=  (h  /  mfi)[aey  +  (N^  -  a' )/'  e  J 

(In  y/^  +lny/J)  =  0  (4.59) 

We  shall  prove  that  each  quamlet  satisfies  the  conservation  equation. 

The  elemental  quamlets,  satisfies  the  continuity  equation,  (4.4),  as  follows, 

=  A\h/m^b)V*{exp[-iay-i(N^  -a")^x][ae^  +  (A^"  -a")^ej} 


=  A/{h/m^b){a{-i){N^ -a^/^]  +  {N^  -a^/^(-ia)}exp[-iay-iiN^  -a")^x]  =  0 

This  proves  the  compliance  of  the  continuity  law  of  crcir^^th  elemental  quamlets,  aa*  and 
Next,  the  interactive  quamlet,  a/3* ,  satisfies  continuity  (4.7),  as  shown  below, 

iV*p^p.(yv^-w^,)  =  0 

=  -{hnl2mJb)A^A^.N  •  {expH(a  -  P)y  -  i{N^  -a")^  -  -  /?')^x]  (4.61) 

x{-/(a  +  ;^)e^ 

After  some  algebraic  manupilation  we  have, 
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•  Paff  -  w ». )  =  -ihn / 2m^b)A^A  -f)-{N^-a^)  +  (N^  - ^")] 

,/  ,/  ('*■«* 

exp {-i{a  +  P)y  - i[{N^  -ay^ +{N^  -  py^x'\}  =  0 


Thus,  the  afi* ,  satisfies  continuity  law.  By  the  similar  algebraic  procedure,  we  can  prove 
that  a*/3  satisfies  the  continuity  equation  (4.1 1) . 

Next  we  illustrate  that  each  type  of  quamlet  satisfies  the  Bernoulli’s  equation  as  follows. 

Based  on  (4.37),  the  elemental  quamlet  obeys  Benoulli’s  equation  , 

2(  V  )  +  {Dlp)V*  [p^^.  (w„  +  W^. )]  =  (^)z^^.  (4.63) 


Following  the  procedure  used  above,  we  substitute  all  the  terms  appear  on  the  left  hand 
side  of  the  above  equation  to  establishe  the  equality,  as  follows. 


2Z) V •  [p^^,  (in +  In )]  =  ) 


(4.64) 


By  substituting  the  density  and  velocities  given  by  (4.54)  in  the  above  equation  and 
expanding  the  left  hand  side  of  above  equation  we  have. 


D^AA.  exp{-iia-/3)y-iiy-a^/^  -(A^^ 


?2\)21 
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D-^Ay.  Qxp{-i{a- P)y-i[y  -a^)X  -y  -,5^)>^]}(2yV^) 

-  D^pyiN^)  =  {n/4my^/bfp^^.  =  p^^.  (E/mJ 


(4.65) 


We  have  shown  that  the  sum  of  the  left  hand  side  of  the  above  equation  agrees  with  the 
right  hand  side  of  Eq.  (4.64),  i.e.. 


2Z)  V  )  .  V  ln(^„  )  -  Z)  V  . 


Pap[^y^a+^¥p) 


=  (— Xv) 

m„  ^ 


(4.66) 


Thus,  we  prove  that  the  quamlet  aP*  satisfies  the  Bernoulli’s  equation  of  aP* 

The  proof  of  the  compliance  of  crcr^^can  be  made  by  replacing  /?«^by  in  the  above 
analysis,  whereas  that  of  a*l3  can  be  demonstrated  by  exchange  the  role  of  a  and  /?  form 
that  of  the  quamlet  aP*. 
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CONCLUSIONS 


The  quantum  fluid  dynamics  (QFD)  formalism  based  on  the  Hamilton-Jacobi 
transformation,  incorporated  with  the  quantum  diffusion,  is  adapted  to  examine  the 
equivalence  between  QFD  non-linear  field  with  quantum  mechanical  description  based 
on  Schrodinger’s  equation  of  linear  wave  field.  The  latter  systems,  which  are 
characterized  by  a  linear  superposition  of  the  wave  functions  of  same  degeneracy  energy 
level,  envisioned  as  the  mixture  of  both  elemental  and  interactive  quantum  flowlets, 
termed  quamlets.  Each  elemental  quamlet  has  the  probability  density  equal  to  the  product 
of  a  degenerate  eigenfunction  and  its  complex  conjugate,  whereas  the  density  of  an 
interactive  quamlet  is  the  product  of  a  degenerate  wave  function  and  the  complex 
conjugate  of  other  degenerate  wave  function.  The  necessary  and  sufficient  conditions  for 
the  equivalence  between  QM  and  QFD  are  that  each  quamlet  obeys  the  laws  of 
conservation:  the  continuity  of  probability,  and  mechanical  energy  of  each  quamlet,  and 
that  the  global  QFD  field  satisfies  the  quantization  condition.  The  quamlets  formalism 
provides  the  theoretical  guides  and  criteria  required  for  the  analytical  and  numerical 
simulation  of  complex  quantum  systems  by  fluid  dynamic  methods,  in  particular,  the 
quantum  computational  fluid  dynamics,  (QCFD).  The  validity  of  the  later  approach  is 
ensured  through  the  implementation  of  the  quamlet  formalism,  which  provides  a  physical 
and  mathematical  equivalence  between  the  quantum  fluid  dynamic  and  Schrodinger’s 
wave  field  descriptions.  Specifically,  the  formalism  ensures  the  compliance  of  the 
quantization,  uncertainty  principles,  and  exclusion  principles,  which  are  embodied  in  the 
wave  formalism. 
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1  Nanojet ,  sine  qua  non  tool  of  quantum 
nanoprocessing 

J  Quantum  nanoJet  theory  of  dual-slit  electron 
interferometer 

j  Structural  and  dynamic  complexities 

-  Quantum  transition:  quantum-classical  transition 

-  Quantum  branching:  multi-branch  configurations 

-  Quantum  clustering:  dynamic  clustering 

-  Interference  fringe  complexities 

J  Conclusions 


I  Definition: 

Quantum  nanojet  refers  to  the  streams  or  arrays  of  particles  of 
atoms,  sub-atoms,  molecules,  neutrons,  ions  or  photons,  injected  or 
emitted  from  a  particle  source,  directly  or  through  a  injector  of 
nanoscale,  to  a  designated  target  space  in  the  presence  or  absence 
of  an  external  potential,  whereby  exhibiting  structural  and  dynamic 
quantum  phenomena. 

1  Quantum  wave  modes; 

Quantum  waves  in  a  jet  exhibit  propagation,  tunneling,  interference, 
reflection  and  refraction,  and  create  hierarchical  jet  structures. 

J  Nanojet:  sine  qua  non  tooi  of  functionai  quantum  processes 

-  Transport  of  particles  and  creation  of  dynamic  field 

-  Formation  of  particle  density  and  flux  distribution  patterns, 

-  Excitations  of  particle-particle  interaction, 

-  Photon-particle  interaction, 

-  Transition  of  particle  states 

-  State  transition  and  formation  of  interference  fringes 
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Figure  B.  This  schemalic  shows  a  modulated  molecular-beam  mass  spectrometer 
bieing  used  to  study  the  effects  of  water  vapor  on  lithium  hydride.  The  beam  of 
water  vapor  direct^  at  the  lithium  hydride  is  turned  on  and  off  at  brief,  regular 
intervals,  as  shown  in  Ihe  square  pattern  of  the  incoming  beam.  The  outgoing 
bieam  of  molecular  hydrogen  has  a  different  shape,  indicating  a  brief  residence 
time  in  the  lithium  hydride. 
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1  Nanojet  Type 

-  Single  injector 

-  Double-slit,  or  Multi-slit  injectors 

-  Free  jet,  confined  jet 

-  External  potential 
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1  Quantum  nanojet  control  param^ers: 

-  Injection  conditions:  particle  Kinetic  energy, 
injector  geometry. 

-  Quantum  Reynolds  Number  N 


1  Physical  significance 

-  Inertial  force/quantum  force 

-  Injector  width/de-Bloglie  wave  length 


Double-slit  quantum  jet  -  Electron  interfeometer 


Quantum  mechanical  analysis  of  two-siit  quantum  jet 
Schrodinger’s  equation  of  Double-slit  jet 
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Density-Velocity  Structures  1 
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density  and  velocity  distributions  viewed  from  toward  upstream 


Density-Velocity  Structures  2 
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1  Processes  characteristic  time 

>  Quantum  diffusion  time:  Td=R/vd=R/pD  [ 

>  Convection  time:  L/u  |  ^ 

Quantum  diffusion  was  T 

introduced  by  Chiu  as  foiiows  2b  I 

vd=-pVlnp  D=h/2m 

:i  Nanojet  behaviors  I 

Quantum  behavior:  Td«  xc  R/v  «L/u,  N«L/R 
Classical  convection:  xd»  xc  R/v  »L/u,  N»L/R 


Djf^u^n  velocity  Vd 


Convection  velocity  u 
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1  Quantum-quasi-classical  transition 

Critical  state:  xd=xc  R/v  =L/u,  N=Ncr=L/R=2b/R 
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''Functional  quantum  processes”  and  quantum  devices  applications 
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I  Two  nanojets,  1  and  2,  with  the  same  slit  geometry  and  inlet 
profiles,  i.e.,  wave  function  and  its  derivative  in  x  direction,  are 
similar  if  the  following  condition  is  satisfied: 

Ni  =N2 

where,  N  is  the  quantum  Reynolds  number,  N  =  mUb/hp. 

Let  b  is  fixed,  then  N  oc  mU,  thus  the  similarity  condition  is 

U2/U1  =  mi/m2 

I  Example,  for  the  same  slit  geometry,  the  structure  of  hydrogen 
atom  Jet,  (ionized)  is  similar  to  that  of  electron. 

For  H  atom  and  electron  we  have  Ue  /UH  =  1380 

For  electron  N=100=mUb/hp,  m=9. 1095x10-31  kg  ,  b=200x10-9 

h  =  1.05  X  10-34  ,  p  =  3.14 

The  electron  velocity  is  Ue  =  100x3. 14x  1.05  x  10-34  /(9.1x10- 
31  )x2x  10-7  =  1.8  x104  m/sec 

The  ionized  hydrogen  with  velocity  of  UH  =  1400m/sec  will  have 
similar  Jet  structure. 


Quantum  nanojets  are  the  smallest  particle  Jets,  which  exhibit  the  fundamentai 
quantum  structurai  and  dynamic  compiexities,  and  offer  broad  appiications  in 
quantum  nanodevices. 

I.  Quantum  nanojets  structural  and  dynamic  complexities 

Quantum  behaviors: 

:  muiti-branches  and  two-branches  (duai-  slit 

nanoJets). 

:  Base  cluster  and  primary  cluster. 
Quantum-quasi-classical  modes  and  transition 

N  <  Ncr  =S/R 

Quasl-classical  Jet  mode  N  >  Ncr  =S/R 

■ 

■ 

Multi-branches  are  formed  by  the  intersecting  right-and 

left  running  quantum  waves  created  in  the  global  jet  field. 

■ 

■ 

Clusters  are  formed  by  the  interacting  right-and-left  quantum  waves  at 
the  axis  of  each  primary  Jet. 

Clusters  are  characterized  by  the  location,  cluster  sizes,  density  peak 
and  the  laws  of  growth  and  decay. 


II.  Nanojet  applications 

Functional  quantum  process  based  applications 
•Cluster  based  applications: 

Mass  spectrometer:  determination  of  atomic  weight 
(Concept  is  proposed,  need  experimental  validation) 

Imaging,  sensor,  lithography,  nanoelectronics. 

Multi-and  two-branch  based  applications: 

Mass  spectrometry:  determination  of  molecular  weights.  (Concept 
has  been  proposed,  need  experimental  validation). 

Interferometer  and  mass  spectrometer  (Concept  has  been 
proposed,  experimental  validation  is  needed. ) 

MLAreas  of  future  research 

Controlled  nanojet,  by  active  and  passive  controls,  for  optimum 
functional  quantum  processes,  at  steady  state  and/or  time 
dependent  is  for  broader  applications. 
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4  New  types  of  quantum  nanojet-based  diagnostic  devices;  mass 
spectrometer  and  energy  sensor,  each  of  which  determines  the 
mass  of  the  particie,  and  the  energy  of  the  particie  of  the 
nanojet  respectively. 

Theoretical  basis  of  quantum  nanojet  diagnostic  devices: 

A  cluster  parameter,  C,  given  as  function  of  the  quantum 
Reynolds  number  N,  by 


ji 


Mass  spectrometer : 


ffj  =^(rjii/2Eb)r(L  -^gj/KJCJ/y) 


Energy  sensor :  ■ 


£  ^  (rjmifw)f(L  ^G)/Kj(i/y) 


The  concepts  described  above  have  not  been  experimentaiiy 
vaiidated.  Further  studies  inciuding  the  experimental  studies 
for  the  vaiidations  and  the  test  of  accuracy  of  these  devices 
must  be  carried  out. 
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Abstract 

Nano-jets  of  atoms,  molecules  or  charged  particles  are 
anticipated  to  play  new  and  important  roles  in  the 
nanoengineering:  nano-  fabrication,  electronics,  and 
avionics;  molecular  assembling;  micro-  and  nano¬ 
propulsion  systems;  atomic  optics  and  interferometery; 
lithography;  drug  delivery  systems;  and  quantum 
sensors.  The  objective  of  this  paper  is  to  understand 
the  fundamentals  driving  the  quantum  dynamics  and 
structural  characteristics  of  nano-jets.  The  ultimate 
goal  is  to  establish  nano-jet  design  and  control 
guidelines  for  specific  applications.  We’ve  conducted 
extensive  numerical  simulations  based  on 
Schrodinger’s  equation,  which  simulate  a  double  slit 
particle  jet  under  various  jet  inlet  conditions.  The 
simulations  are  characterized  by  the  quantum  a 
Reynolds  number,  defined  as  the  ratio  of  the  inertia 
force  to  quantum  force,  N  =  mUb/h,  where  m  is  the 
mass  of  particle,  U  is  the  inlet  velocity,  b  is  the 
characteristic  jet  size,  and  h  is  the  Planck  constant. 
With  the  slit  separation  and  widths  fixed  at  of  200  nm 
and  20  nm,  respectively,  we  conducted  parametric 
studies  varying  the  range  of  N  from  10  to  1000, 
corresponding  to  the  electron  initial  velocity  of  10^  to 
10^  m/s,  to  investigate  the  nature  of  the  dynamic  and 
structural  variations  of  the  jet  streams.  At  lower 
Reynolds  number,  0  <  N  <  100,  or  0  <  U  <  lO^m/s, 
the  dual-slit  jet  reveals  multi-branching  configurations. 
For  example,  when  N=10,  Fig. la,  quantum  branching 
is  observed  at  a  distance  from  a  vertical  plane  passing 
through  the  jet  nozzles  that  is  comparable  to  the 
separation  width  of  the  slits.  Branching  is  a  result  of 
the  generation  of  higher  harmonic  probability  density 
waves  within  the  stream.  The  center,  or  forward 
propagating  branch,  is  a  result  of  the  interaction  of  the 
streams  effusing  from  each  slit.  Branching,  as 
expected,  has  reflection  symmetry  through  a 


horizontal  plane  dividing  the  two  slits.  As  N  increases 
the  number  of  branches  increases  and  their  angular 
displacement  with  respect  to  the  plane  of  symmetry 
decreases.  Moreover,  the  center  branch  is  generated 
at  greater  and  greater  distances  from  the  jet  nozzle(s). 
By  the  time  the  Reynolds  number  reaches  1000  the 
center  branch  occurs  at  inflnity  and  the  angular 
displacements  of  the  side  branches  have  converged  on 
each  other  producing  two  distinct  forward  propagating 
beams  of  particles  (See  Fig.  1).  A  uniquely  interesting 
phenomenon,  which  we  call  quantum  clustering,  is 
found  over  a  broad  range  of  Reynolds  numbers.  When 
N  =  100  the  probability  density  in  the  region  of  0.5  < 
X  <  1.7,  or  50  nm  <  x  <  170  nm,  increases  to  the 
maximum  value  of  0.9.  This  is  twice  its  value  found  at 
the  plane  of  the  jet  nozzles  and  is  an  example  of  the 
quantum  clustering  of  the  probability  current  density. 
The  detailed  structures  of  these  clusters  vary  as 
function  of  Reynolds  number.  Clustering  is 
characterized  by  an  initial  axial  growth  in  the 
probability  density  immediately  after  the  two  branches 
intersect  on  the  reflection  plane  of  symmetry,  and  thus 
constitutes  the  major  self-interference  of  a  jet.  The 
first  cluster,  or  base  cluster,  decays  rapidly  but  is 
followed  immediately  by  another  rapid  compression  to 
peak  density  to  form  a  primary  cluster,  which  has  a 
larger  amplitude  than  the  base  cluster.  The  primary 
cluster  decays  exponentially  due  to  a  quantum 
tunneling  process,  which  is  a  diffusion  transport 
process  associated  with  the  low  linear  momentum  of 
the  beam  particles  in  the  forward  direction.  The  region 
dominated  by  group  wave  propagation  with  long 
group  wavelength  as  shown  in  Fig.  1,  follows  the 
primary  cluster.  The  wavelength  of  the  wavelet  in  the 
base  and  primary  cluster  and  the  group  wave 
propagation  region  are  approximately  equal.  The 
structural  variation  of  a  single-slit  nanojet  with  respect 
to  the  quantum  Reynolds  number  is  also  presented  in 
Fig.2.  Numerical  simulations  for  the  low,  and 
intermediate  energy  ranges  reveal  profoundly  complex 
structural  variations  of  the  interference  fringe, 
quantum  potential,  and  velocity  distributions  as  shown 
in  Fig.  3.  Finally,  we  conclude  that  nanojets  in  a  broad 
range  of  Reynolds  number  lend  themselves  to  various 
nanodevices  applications,  depicted  in  Fig.4.  Here  N* 
is  the  critical  Reynolds  number  at  which  the  system 
transitions  from  a  quantum  to  a  classical  regime. 
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Fig  3.  Interference  fringes,  quantum  potentials  and  velocity  distributions  at  various  N 
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Fig  4.  Application  of  nanojets  in  nanodevices 
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